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ABSTRACT 

The theory of phase ordering dynamics - the growth of order through domain 
coarsening when a system is quenched from the homogeneous phase into a broken- 
symmetry phase - is reviewed, with the emphasis on recent developments. Interest 
will focus on the scaling regime that develops at long times after the quench. How 
can one determine the growth laws that describe the time-dependence of character- 
istic length scales, and what can be said about the form of the associated scaling 
functions? Particular attention will be paid to systems described by more compli- 
cated order parameters than the simple scalars usually considered, e.g. vector and 
tensor fields. The latter are needed, for example, to describe phase ordering in ne- 
matic liquid crystals, on which there have been a number of recent experiments. The 
study of topological defects (domain walls, vortices, strings, monopoles) provides a 
unifying framework for discussing coarsening in these different systems. 
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1 INTRODUCTION 



The theory of phase ordering kinetics, or 'domain coarsening' following a tempera- 
ture quench from a homogeneous phase into a two-phase region, has a history going 
back more than three decades to the pioneering work of Lifshitz and coworkers 
, 0, ^. Since that time, many excellent reviews have appeared, including those by 
unton et al. Binder [Q, Furukawa and Langer In the present article I 
will not, therefore, attempt to cover the complete history of the field. Rather, I will 
concentrate on some of the recent developments, over the last five years or so, which 
in my view are interesting and represent significant advances or new directions of 
research, e.g. the recent interest in systems with non-scalar order parameters. The 
fundamental concepts and background material necessary for the understanding and 
appreciation of these new developments will, nevertheless, be discussed in some de- 
tail. It follows that, while this article does not aim to be a complete or comprehensive 
account, it does aspire to be self-contained and comprehensible to non-experts. By 
adopting a fairly pedagogical approach, I hope that the article may also serve as a 
useful introduction to the field. 

In order to keep the length of the article within reasonable bounds, I will concen- 
trate primarily on theoretical developments, although important results from exper- 
iment and simulations will be cited as appropriate. For the same reason, I apologise 
in advance to all those authors whose work has not been explicitly discussed. 

Systems quenched from a disordered phase into an ordered phase do not order 
instantaneously. Instead, the length scale of ordered regions grows with time as 
the different broken symmetry phases compete to select the equilibrium state. To 
fix our ideas, it is helpful to consider the simplest, and most familiar, system: the 
ferromagnetic Ising model. Figure 1 shows the spontaneous magnetization as a 
function of temperature. The arrow indicates a temperature quench, at time t = 0, 
from an initial temperature Tj above the critical point Tq to a final temperature 
Tp below Tq- At Tp there are two equilibrium phases, with magnetization ±Mq. 
Immediately after the quench, however, the system is in an unstable disordered 
state corresponding to equilibrium at temperature T/. The theory of phase ordering 
kinetics is concerned with the dynamical evolution of the system from the initial 
disordered state to the final equilibrium state. 

Part of the fascination of the field, and the reason why it remains a challenge 
more than three decades after the first theoretical papers appeared, is that, in the 
thermodynamic limit, final equilibrium is never achieved! This is because the longest 
relaxation time diverges with the system size in the ordered phase, reflecting the 
broken ergodicity. Instead, a network of domains of the equilibrium phases develops, 
and the typical length scale associated with these domains increases with time t. 
This situation is illustrated in Figure 2, which shows a Monte Carlo simulation of a 
two-dimensional Ising model, quenched from Tj = oo to Tp = 0. Inspection of the 
time sequence may persuade the reader that domain growth is a scaling phenomenon 
- the domain patterns at later times look statistically similar to those at earlier 
times, apart from a global change of scale. This 'dynamic scaling hypothesis' will 
be formalized below. 

For pedagogical reasons, we have introduced domain-growth in the context of the 
Ising model, and will continue to use magnetic language for simplicity. A related 
phenomenon that has been studied for many decades, however, by metallurgists, 
is the spinodal decomposition of binary alloys, where the late stages of growth are 
known as 'Ostwald ripening'. Similar phenomena occur in the phase separation of 
fluids or binary liquids, although in these cases the phase separation is accelerated 
by the earth's gravitational field, which severely limits the temporal duration of 
the scaling regime. The gravitational effect can be moderated by using density- 
matched binary liquids and/or performing the experiments under microgravity. All 
of the above systems, however, contain an extra complication not present in the 
Ising ferromagnet. This is most simply seen by mapping an AB alloy onto an Ising 
model. If we represent an A atom by an up spin, and a B atom by a down spin, 
then the equilibrium properties of the alloy can be modelled very nicely by the Ising 
model. There is one important feature of the alloy, however, that is not captured by 
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the Ising model with conventional Monte-Carlo dynamics. Flipping a single spin in 
the Ising model corresponds to converting an A atom to a B atom (or vice versa), 
which is inadmissible. The dynamics must conserve the number of A and B atoms 
separately, i.e. the magnetization (or 'order parameter') of the Ising model should 
be conserved. This will influence the form of the coarse-grained equation of motion, 
as discussed in section ^ and lead to slower growth than for a non-conserved order 
parameter. 

In all the systems mentioned so far, the order parameter (e.g. the magnetization 
of the Ising model) is a scalar. In the last few years, however, there has been 
increasing interest in systems with more complex order parameters. Consider, for 
conceptual simplicity, a planar ferro magnet, in which the order parameter is a vector 
confined to a plane. After a quench into the ordered phase, the magnetization 
will point in different directions in different regions of space, and singular lines 
(vortex lines) will form at which the direction is not well defined. These 'topological 
defects' are the analog of domain walls for the scalar systems. We shall find that, 
quite generally, an understanding of the relevant topological defects in the system, 
combined with the scaling hypothesis, will take us a long way towards understanding 
the forms of the growth laws and scaling functions for phase ordering in a wide 
variety of systems. 

The article is organised as follows. The following section introduces most of the 
important concepts, presents dynamical models appropriate to the various phase- 
ordering systems, and analyses these models using simple physical arguments. Sec- 
tion ^ broadens the discussion to more general phase-ordering systems, with non- 
scalar order parameters, and introduces the key concept of 'topological defects' 
which, in later sections, will provide a unifying framework for analytic treatments. 
Section ^ involves a temporary excursion to the realm of exactly soluble models. 
These models, although of interest in their own right, unfortunately lack many of 
the physical features of more realistic models. In section |^, approximate analytical 
treatments are presented for the more physical models introduced in sections |] and 
R|. Finally, sections ^, |^ and S present some exact results for the short-distance 
behaviour and the growth laws for these systems, and make some observations con- 
cerning universality classes for the dynamics of phase ordering. 



2 DYNAMICAL MODELS 

It is convenient to set up a continuum description in terms of a coarse-grained order- 
parameter field (e.g. the 'magnetization density') (^(x, t), which we will initially take 
to be a scalar field. A suitable Landau free-energy functional to describe the ordered 
phase is 

m = J d'x (V0)^ + V^(0)) , (1) 

where the 'potential' V{(j)) has a double-well structure, e.g. V{(f)) = (1 — 0^)^. We 
will take the minima of V{(j)) to occur at = ±1, and adopt the convention that 
V{±1) = 0. The potential V{(f)) is sketched in Figure 3. The two minima of V 
correspond to the two equilibrium states, while the gradient-squared term in (|I]) 
associates an energy cost to an interface between the phases. 

In the case where the order parameter is not conserved, an appropriate equation 
for the time evolution of the field is 

d(f)/dt = -6F/6(f) 

= v2 0-n0), (2) 

where V'{(j)) = dV/dcf). A kinetic coefficient F, which conventionally multiplies 
the right-hand side of (^, has been absorbed ibto the timescale. Eq. (Q), a sim- 
ple 'reaction-diffusion' equation, corresponds to simple gradient descent, i.e. the 
rate of change of is proportional to the gradient of the free-energy functional in 
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function space. This equation provides a suitable coarse-grained description of the 
Ising model, as well as alloys that undergo an order- disorder transition on cooling 
through Tc, rather than phase separating. Such alloys form a two-sublattice struc- 
ture, with each sublattice occupied predominantly by atoms of one type. In Ising 
model language, this corresponds to antiferromagnetic ordering. The magnetization 
is no longer the order parameter, but a 'fast mode', whose conservation does not 
significantly impede the dynamics of the important 'slow modes'. 

When the order parameter is conserved, as in phase separation, a different dy- 
namics is required. In the alloy system, for example, it is clear physically that A and 
B atoms can only exchange locally (not over large distances), leading to diffusive 
transport of the order parameter, and an equation of motion of the form 

d(j)/dt = V^6F/6(j) 

= -V^ [V^ - m] , (3) 

which can be written in the form of a continuity equation, dt(f) = —V ■ j, with 
current j = —X'V{6F/6(f)). In (|^), we have absorbed the transport coefficient A into 
the timescale. 

Eqs. (0) and (0) are sometimes called the Time-Dependent-Ginzburg-Landau 
(TDGL) equation and the Cahn-Hilliard equation respectively. A more detailed 
discussion of them in the present context can be found in the lectures by Langer 
1^. The same equations with additional Langevin noise terms on the right-hand 
sides are familiar from the theory of critical dynamics, where they are 'model A' 
and 'model B' respectively in the classification of Hohenberg and Halperin 

The absence of thermal noise terms in (|]) and (|^) indicates that we are effectively 
working at T = 0. A schematic Renormalization Group (RG) flow diagram for T 
is given in Figure 4, showing the three RG fixed points at 0, Tq and oo, and the 
RG flows. Under coarse-graining, temperatures above Tq flow to infinity, while 
those below Tc flow to zero. We therefore expect the final temperature Tp to be 
an irrelevant variable (in the scaling regime) for quenches into the ordered phase 
This can be shown explicitly for systems with a conserved order parameter |10 



For this case the thermal fluctuations at Tp simply renormalize the bulk order 
parameter and the surface tension of the domain walls: when the characteristic 
scale of the domain pattern is large compared to the domain wall thickness (i.e. the 
bulk correlation length in equilibrium), the system behaves as if it were T = 0, with 
the temperature dependence entering through T-dependent model parameters. 

In a similar way, any short-range correlations present at Tj should be irrelevant in 
the scaling regime, i.e. all initial temperatures are equivalent to Tj = oo. Therefore 
we will take the initial conditions to represent a completely disordered state. For 
example, one could choose the 'white noise' form 

(0(x,O)0(x',O)) = A5(x-x') , (4) 

where (■ ■ ■) represents an average over an ensemble of initial conditions, and A 
controls the size of the initial fluctuations in 0. The above discussion, however, 
indicates that the precise form of the initial conditions should not be important, as 
long as only short-range spatial correlations are present. 

The challenge of understanding phase ordering dynamics, therefore, can be posed 
as finding the nature of the late-time solutions of deterministic differential equations 
like d^) and (^, subject to random initial conditions. A physical approach to this 
formal mathematical problem is based on studying the structure and dynamics of 
the topological defects in the field (p. This is approach that we will adopt. For scalar 
fields, the topological defects are just domain walls. 



2.1 THE SCALING HYPOTHESIS 

Although originally motivated by experimental and simulation results for the struc- 
ture factor and pair correlation function 0, 12, 13, 14] for ease of presentation it 
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is convenient to introduce the scaling hypothesis first, and then discuss its imphca- 
tions for growth laws and scaling functions. Briefly, the scaling hypothesis states 
that there exists, at late times, a single characteristic length scale L{t) such that 
the domain structure is (in a statistical sense) independent of time when lengths are 
scaled by L{t). It should be stressed that scaling has not been proved, except in 
some simple models such as the one-dimensional Glauber model |]T5[ and the n- vector 
model with n = oo |T^. However, the evidence in its favour is compelling (see, e.g.. 
Figure 5). 

We shall flnd, in section |7, that the scaling hypothesis, together with a result 
derived in section ^ for the tan of the structure factor, is sufficient to determine the 
form of L{t) for most cases of interest. 

Two commonly used probes of the domain structure are the equal-time pair 
correlation function 

C(r,t) = (0(x + r,t)0(x,t)), (5) 
and its Fourier transform, the equal-time structure factor. 



5(k,t) = (0k(t)0-k(t)) 



(6) 



Here angle brackets indicate an average over initial conditions. The structure factor 
can, of course, be measured in scattering experiments. The existence of a single 
characteristic length scale, according to the scaling hypothesis, implies that the pair 
correlation function and the structure factor have the scaling forms 



C(r,t) 



fir/L) , 



(7) 



where d is the spatial dimensionality, and g{y) is the Fourier transform of f{x). 
Note that /(O) = 1, since (at T = 0) there is perfect order within a domain. 

At general temperatures T < Tc, /(O) = M^, where M is the equilibrium value 
of the order parameter. (Note that the scaling limit is defined hy r ^ ^, L ^ ^, 
with r/L arbitrary, where ^ is the equilibrium correlation length). Alternatively, we 
can extract the factor explicitly by writing C(r, t) = f{r/L). The statement 
that T is irrelevant then amounts to asserting that any remaining temperature de- 
pendence can be absorbed into the domain scale L, such that the function f{x) is 
independent of T. 

The scaling forms (|^) are well supported by simulation data and experiment. As 
an example. Figure 5 shows the scaling plot for f{x) for the 2-D Ising model, with 
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For future reference, we note that the different-time correlation function, defined 
by C{r,t,t') = (0(x + r, t) J)(x., t')), can also be written in scaling form. A simple 
generalization of (0) gives |l8l 



C{r,t,t') = f{r/L,r/L') 



(8) 



where L, V stand for L{t) and Lit'). Especially interesting is the limit L ^ L', 
when (^) takes the form 



C{v,t,t')^{L'/Lfh{r/L) 



(9) 



where the exponent A, first introduced by Fisher and Huse in the context of non- 
equilibrium relaxation in spin glasses |]T9i, is a non-trivial exponent associated 



with phase ordering kinetics p^. It has recently been measured in an experi- 



ment on twisted nematic liquid crystal films pTj. The autocorrelation function 



A{t) = C{0,t,t') is therefore a function only of the ratio L'/L, with A{t) ~ (L' / L)^ 
for L > L'. 

In the following sections, we explore the forms of the scaling functions in more 
detail. For example, the linear behaviour of f{x), for small scaling variable x in 
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Figure 5, is a generic feature for scalar fields, both conserved and non-conserved. 
We shall see that it is a simple consequence of the existence of 'sharp' (in a sense 
to be clarified), well-defined domain walls in the system. A corollary that we shall 
demonstrate is that the structure factor scaling function g{y) exhibits a power-law 
tail, g{y) ~ y-(^+^') for y ^ 1, a result known as 'Porod's law' ||22|, In section |^ 
we shall show that this result, and its generalization to more complex fields, together 
with the scaling hypothesis, are sufficient to determine the growth law for L{t). 



2.2 DOMAIN WALLS 

It is instructive to first look at the properties of a fiat equilibrium domain wall. 
From d^) the wall profile is the solution of the equation 

d^Hdg' = , (10) 

with boundary conditions 0(±oo) = ±1, where g is a coordinate normal to the wall. 
We can fix the 'centre' of the wall (defined by = 0) to be at = by the extra 
condition 0(0) = 0. Integrating (|10J) once, and imposing the boundary conditions, 
gives {d(f)/dgY = 2V{(t>). This result can be used in (^ to give the energy per unit 
area of wall, i.e. the surface tension, as 

a = dg {dcP/dgf = £ d<P pV{ct>) . (11) 

Note that, for scalar fields, the two terms in (P contribute equally to the wall energy. 

The profile function (t){g) is sketched in Figure 6. For (7 — ^ ±oo, linearizing ([lOD 
around = ±1 gives 

1 T ~ exp(-[y"(±l)]i/2|^|) , ^7 ^ ±00 , (12) 

i.e. the order parameter saturates exponentially fast away from the walls. It follows 
that the excess energy is localized in the domain walls, and that the driving force for 
the domain growth is the wall curvature, since the system energy can only decrease 
through a reduction in the total wall area. The growth mechanism is rather different, 
however, for conserved and nonconserved fields. 



2.3 NONCONSERVED FIELDS: THE ALLEN-CAHN EQUA- 
TION 

The existence of a surface tension implies a force per unit area, proportional to 
the mean curvature, acting at each point on the wall. The calculation is similar 
to that of the excess pressure inside a bubble. Consider, for example, a spherical 
domain of radius in three dimensions. If the force per unit area is F, the work 
done by the force in decreasing the radius by dR is AnFR'^dR. Equating this to 
the decrease in surface energy, SnaRdR, gives F = 2a /R. For model A dynamics, 
this force will cause the walls to move, with a velocity proportional to the local 
curvature. If the friction constant for domain-wall motion is rj, then this argument 
gives rjdR/dt = —2a/R. For general dimension d, the factor '2' on the right is 
replaced by {d — 1). 

It is interesting to see how this result arises directly from the equation of motion 
(0). We consider a single spherical domain of (say) = — 1 immersed in a sea of 
+1. Exploiting the spherical symmetry, (^ reads 

dt dr'^ r dr 
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Provided the droplet radius R is much larger than the interface width ^, we expect 
a solution of the form 

<P{r,t) = f {r - R{t)) . (14) 

Inserting this in ([T3|) gives 

= r + [{d - l)/r + dR/dt]f - V'if) . (15) 

The function f{x) changes from -1 to 1 in a small region of width near x = 0. It's 
derivative is, therefore, sharply peaked near x = (i.e. near r = R{t)). Multiplying 
(P^) by /' and integrating through the interface gives 

= {d-l)/R + dR/dt , (16) 

where we have used /' = far from the interface, and V{f) has the same value on 
both sides of the interface (in the absence of a bulk driving force, i.e. a magnetic 
field). Integrating (|T^) gives R'^it) = _R^(0) — 2{d — l)t, i.e. the collapse time scales 
with the initial radius as t ~ i?^(0). Equation (O) is identical to our previous result 
obtained by considering the surface tension as the driving force, provided the surface 
tension a and friction constant r] are equal. This we show explicitly below. 



The result for general curved surfaces was derived by Allen and Cahn |Q, who 
noted that, close to a domain wall, one can write = {d(f)/dg)t g, where g is 
a unit vector normal to the wall (in the direction of increasing 0), and so = 
{d'^(j)/dg'^)t + {d(l)/dg)t V-g- Noting also the relation {d(j)/dt)g = —{d(p/dg)t {dg/dt)^, 
(0) can be recast as 

- {d<P/dg)t (dg/dt)^ = {d<P/dg)t V-g + {d^<i^/dg\ - V\<^) . (17) 

Assuming that, for gently curving walls, the wall profile is given by the equilibrium 
condition (lU), the final two terms in ( [T7| ) cancel. Noting also that [dg/dt)^ is just 



the wall velocity v (in the direction of increasing 0), ( p!7D simplifies to 

v = -V-g = -K , (18) 

the 'AUen-Cahn equation', where K = V ■ g is {d — 1) times the mean curvature. 
For brevity, we will call K simply the 'curvature'. An alternative derivation of ([THD 
follows the approach used for the spherical domain, i.e. we multiply Eq. ([T7|) by 
{d(j)/dg)t and integrate (with respect to g) through the interface. This gives the 
same result. 

Equation (0) is an important result, because it establishes that the motion of the 
domain walls is determined (for non-conserved fields) purely by the local curvature. 
In particular, the detailed shape of the potential is not important: the main role of 
the double- well potential V{(f)) is to establish (and maintain) well-defined domain 
walls. (Of course, the well depths must be equal, or there would be a volume driving 
force). We shall exploit this insensitivity to the potential, by choosing a particularly 
convenient form for V{(j)), in section ^ 

For a spherical domain, the curvature K is {d — l)/R, and ([TB|) reduces to (0). 
Our explicit treatment of the spherical domain verifies the Allen-Cahn result, and, 
in particular, the independence from the potential of the interface dynamics. 

A second feature of ( |TBD is that the surface tension a (which does depend on the 
potential) does not explicitly appear. How can this be, if the driving force on the 
walls contains a factor a? The reason, as we have already noted, is that one also 
needs to consider the friction constant per unit area of wall, t]. The equation of 
motion for the walls in this dissipative system is rjv = —aK. Consistency with ([ISD 
requires r] = a. In fact, r] can be calculated independently, as follows. Consider a 
plane wall moving uniformly (under the influence of some external driving force) at 
speed V. The rate of energy dissipation per unit area is 

dE/dt = / dg 

J —oo 



-oo UlfJ dt 



(19) 
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using (0). The wall profile has the form (f){g, t) = f{g—vt), where the profile function 
/ will, in general, depend on v. Putting this form into ( jlQ] ) gives 

dE/dt = -v"^ j dg {d(t)/dgf = -av^ , (20) 

where the definition ([TT| ) of the surface tension cr was used in the final step, and the 
profile function f{x) replaced by its f = form to lowest order in v. By definition, 
however, the rate of energy dissipation is the product of the frictional force t]v and 
the velocity, dE/dt = —riv"^. Comparison with (pO|) gives r] = a. We conclude 
that, notwithstanding some contrary suggestions in the literature, the Allen-Cahn 
equation is completely consistent with the idea that domain growth is driven by the 
surface tension of the walls. 



2.4 CONSERVED FIELDS 

For conserved fields the interfaces cannot move independently. At late times the 
dominant growth mechanism is the transport of the order parameter from interfaces 
of high curvature to regions of low curvature by diffusion through the intervening 
bulk phases. To see how this works, we first linearize (0) in one of the bulk phases, 

with say ~ 1. Putting = 1 + in (^), and linearizing in 0, gives 

d4>/dt = -V'0 + V"(1)V20 . (21) 

Since the characteristic length scales are large at late times, the term is negligible, 
and (pi]) reduces to the diffusion equation, with diffusion constant D = ^"(1)- The 
interfaces provide the boundary conditions, as we shall see. However, we can first 
make a further simplification. Due to the conservation law, the interfaces move little 

during the time it takes the diffusion field to relax. If the characteristic domain 
size is L, the diffusion field relaxes on a time scale tn ~ L^. We shall see below, 
however, that a typical interface velocity is of order 1/L^, so the interfaces only 
move a distance of order unity (i.e. much less than L) in the time t^- This means 
that the diffusion field relaxes quickly compared to the rate at which the interfaces 
move, and is essentially always in equilibrium with the interfaces. The upshot is 

that the diffusion equation can be replaced by Laplace's equation, V^0 = 0, in the 
bulk. 

To derive the boundary conditions at the interfaces, it is convenient to work, not 

with directly, but with the chemical potential fi = 6F/6(j). In terms of fi, (^) can 
be written as a continuity equation, 

d(f)/dt = -V-j (22) 
J = -V/i (23) 
l_i = V'(0) - V20 . (24) 

In the bulk, fi and are proportional to each other, because (^) can be linearized 

to give fi = V"{l)(j) — V^0, and the term is again negligible. Therefore fi also 
obeys Laplace's equation, 

VV = , (25) 

in the bulk. 

The boundary conditions are derived by analysing (E^ near an interface. As 
in the derivation of the Allen-Cahn equation, we consider surfaces of constant 
near the interface and introduce a Cartesian coordinate system at each point, with 
a coordinate g normal to the surface (and increasing with increasing 0). Then ( P^ 
becomes (compare Eq. (|T7p, 

/i = V'(0) - {d(P/dg)t K - {d^/dg\ (26) 
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near the interface, where K = 'V ■ g is the curvature. The value of jj, at the inte rfac e 
can be obtained (just as in our treatment of the spherical domain in section |2.3| ), 
by multiplying through by {d(f)/dg)t, which is sharply peaked at the interface, and 
integrating over g through the interface. Noting that /i and K vary smoothly through 
the interface, this gives the completely general relation 

/iA0 = AV-aK (27) 

at the interface, where A(f) is the change in (p across the interface, and AV is the 
difference in the minima of the potential for the two bulk phases. In deriving (^), 
we have used {d(f)/dg)t far from the interface, and made the identification 
/ dg{d(t)/dg)1 = a, as in ([TT|), with a the surface tension. Simplifying to the case 
where the minima have equal depth (we shall see that the general case introduces 
no new physics), and taking the minima to be at = ±1 as usual, gives AV = 
and A0 = 2. Then (^) becomes 

/i = -aK/2 . (28) 

This (or, more generally, Eq. (^7|)) is usually known as the Gibbs- Thomson boundary 
condition. Note that we have assumed that the order parameter takes its equilibrium 
value (±1) in both bulk phases. This is appropriate to the late stages of growth in 
which we are primarily interested. 

To summarize, (^) determines fx on the interfaces in terms of the curvature. 
Between the interfaces, fj, satisfies the Laplace equation (^5|) . The final step is to 
use (^) to determine the motion of the interfaces. An interface moves with a velocity 
given by the imbalance between the current flowing into and out of it: 

vA(f) = jout - Jin = -[dfi/dg] = -[g-Vfi] , (29) 

where v is the speed of the interface in the direction of increasing 0, g is the usual 
coordinate normal to interface, [■ ■ ■] indicates the discontinuity across the interface, 
and we have assumed as usual that ~ ±1 in the bulk phases. 

To illustrate how (pSj), (^) and (|29|) are used, we consider again the case of a 
single spherical domain of negative pEase (0 = —1) in an infinite sea of positive 
phas e (0 = +1). We restrict ourselves to d = 3 for simplicity. The definition of fi, 
Eq. (M), gives = at infinity. Let the domain have radius R{t). The solution of 
(p5|) mat obeys the boundary conditions yU = at infinity and (^) at r = R, and 
respects the spherical symmetry is (using K = 2/R for d = 3) fi = —a/r for r > R. 
Inside the domain, the 1/r term must be absent to avoid an unphysical singularity 
at r = 0. The solution of (^) in this region is therefore fi = const. The boundary 
condition (^Sj) gives /i = —ajR. 



To summarize. 



/X = —a/R , r < R 

= -a/r , r>R. (30) 



Using (29), with A0 = 2, then gives 



dR/dt = v = ~[dfi/dr]^tl = -a/2R^ , (31) 

and hence R^{t) = R^{0) — 3cxt/2. We conclude that a domain of initial radius R{0) 
evaporates in a time proportional to i?^(0). This contrasts with the i?^(0) result 
obtained for a non-conserved order parameter. In the non-conserved case, of course, 
the domain simply shrinks under the curvature forces, whereas for the conserved 
case it evaporates by the diffusion of material to infinity. 

We now briefly discuss the case where the potential minima have unequal depths, 
as sketched in Figure 7. Consider first a planar interface separating the two equi- 
librium phases, with order parameter values 0i and 02. Since no current flows. 
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j = — V/i = gives /i = constant. From the definition p^ ) of /i, and tlie fact that 
vanishes far from the interface, it follows that fi = V'{(j)i) = V'{(j)2)- On the 
other hand, the Gibbs-Thomson boundary condition (|27|) for a flat interface {K = 0) 
gives n = AV/A(j). Combining these two results gives 

= V'{<f)2) = AV/A<j) , (32) 

leading to the common tangent construction, shown in Figure 7, that determines 0i 
and 02 as the points where the common tangent touches the potential. If one now 
repeats the calculation for a spherical drop, with a domain with = 0i immersed 
in a sea with = 02, one obtains the equation of motion for the radius, dR/dt = 
— 2(T/(A0)^i?^, a simple generalisation of (|3l|) . Henceforth, we will consider only 
the case of degenerate minima. 



2.5 GROWTH LAWS 

The scaling hypothesis suggests a simple intuitive derivation of the 'growth laws' for 
L{t), which are really just generalizations of the calculations for isolated spherical 
domains. For model A, we can estimate both sides of the Allen-Cahn equation 
(p!8|) as follows. If there is a single characteristic scale L, then the wall velocity 
V ~ dL/dt, and the curvature ~ 1/L. Equating and integrating gives L{t) ~ t^/^ 
for non-conserved scalar fields. 

For conserved fields (model B), the argument is slightly more subtle. We shall 
follow the approach of Huse |25|. From (pSD, the chemical potential has a typical 
value /i ~ (t/L on interfaces, and varies over a length scale of order L. The current, 
and therefore the interface velocity v, scale as V/i ~ cr/L"^, giving dL/dt ~ a/L"^ 
and L{t) ~ (aty^^. A more compelling argument for this result will be given in 
section |^. We note, however, that the result is supported by evidence from com- 
puter simulations [B^, ^ (which usually require , however, some extrapolation into 
the asymptotic scaling regime) as well as a Renormalization Group (RG) treatment 
0, |T0[ . In the limit that one phase occupies an infinitesimal volume fraction, the orig- 
inal Lifshitz-Slyozov- Wagner theory convincingly demonstrates a t^^^ growth. This 
calculation, whose physical mechanism is the evaporation of material (or magneti- 
zation) from small droplets and condensation onto larger droplets, will be discussed 
briefly in the following subsection. 

It is interesting that these growth laws can also be obtained using naive ar- 
guments based on the results for single spherical domains |^. For nonconserved 
dynamics, we know that a domain of radius R collapses in a time of order R^. 
Therefore, crudely speaking, after time t there will be no domains smaller than t^/^, 
so the characteristic domain size is L(t) ~ t^/^. Of course, this is an oversimplifica- 
tion, but it captures the essential physics. For conserved dynamics, the same line of 
argument leads to t^^^ growth. In fact, this approach can be used rather generally. 



for a variety of systems , and gives results which agree, in nearly all cases, with 
the exact growth laws that will be derived in section |^. 



2.6 THE LIFSHITZ-SLYOZOV- WAGNER THEORY 

In their seminal work, Lifshitz and Slyozov, and independently Wagner, derived 
some exact results in the limit that the minority phase occupies a negligible volume 
fraction [0, |^. in particular, they showed that the characteristic size of the minority 

phase droplets increases like t^^^. 

We begin by considering again a single spherical droplet of minority phase (0 = 
— 1), of radius R, immersed in a sea of majority phase, but now we let the majority 
phase have order parameter = 0o < 1 at infinity, i.e. the majority phase is 
'supersaturated' with the dissolved minority species. If the minority droplet is large 
enough it will grow by absorbing material from the majority phase. Otherwise it will 
shrink by evaporation as before. A 'critical radius' Rc separates these two regimes. 
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With the convention that V{±1) = 0, the boundary condition ( p?] ) at r = i? 
becomes (1 + 0o)/^ = ^(0o) ~ 2a /R, while the boundary condition at infinity is 
II = V'{(j)o). Solving the Laplace equation for /i with these boundary conditions 
gives 

^ - n.„).(fM_n.o))^^i, r>_R^ (33) 
. ^-^i r<fi. (34) 

1 + 00 1 + 00 -K 

Eq. ( p9|) gives the interface velocity, 

dR ^ ( l^(0o) _ n0o) \ ^ _ 2a l_ 
dt V(l + 0o)2 1 + 00/ R (l + 0o)2 R'' ■ 



(35) 



Now consider the limit of small supersaturation, 0o = 1 — e with e -C 1. To 
leading non-trivial order in e, the velocity is 

dR cr / 1 1 , , , 

(36) 



dt 2R \Rc R 

where Rc = a/V"{l)e is the critical radius. 

In the Lifshitz-Slyozov- Wagner (LSW) theory, an assembly of drops is considered. 
Growth proceeds by evaporation of drops with R < Rc and condensation on to 
drops with R > Rc. The key idea is to use the time-dependent supersaturation e{t) 
as a kind of mean- field, related to the time-dependent critical radius via Rc{t) = 
a/V"{l)e{t), and to use ( p6| ) with the time- dependent Rc for the dynamics of a given 
drop. 

So far the discussion has been restricted to spatial dimension d = 3. However, 
a result of the form (^) can be derived (with a d-dependent numerical constant 
multipying the right-hand side) for general d > 2. The next step is to write down a 
scahng distribution of droplet radii, 

obeying the continuity equation 

W + 4 = ' (38) 

where v{R) is just the velocity dR/dt. 

Suppose the spatial average of the order parameter is 1 — eo- At late times the 
supersaturation e{t) tends to zero, giving the constraint 

/■OO /"OO 

eo = 2Vd / dRR'^n{R,t) = 2Vd / dxx'^f{x) , (39) 
JO Jo 



where Vd is the volume of the d-dimensional unit sphere. Eq. (|39|) fixes the normal- 
isation of f{x). A linear equation for f {x) can be derived by inserting the scaling 
form (|37D into the continuity equation (^) and using (|36D for v{R), which we write 
in the form (valid for general d > 2) 



R \Rr R 
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where = a/2. This procedure gives 



where Rc = dRc/dt and /' = df/dx. A consistent solution requires that the Rc- 

dependence drop out from this equation. This means that R^Rc = c^dl, a constant, 
giving 

Rcit) = iSadit)'/' , (42) 

Integrating (^l]) then gives 

,^ .^d,2-y-^(d^ (43) 

It is clear that f{x) cannot be non-zero for arbitrarily large x, or one would have 
the asymptotic behaviour /(x) ~ x~'^'''^^\ and the normalisation integral (|39|) would 
not exist. Therefore f{x) must vanish for x greater than some value Xmax-, which 
must be the first pole of the integrand on the positive real axis. The existence of 
such a pole requires 7 < 4/27. Lifshitz and Slyozov argue that the only physically 
acceptable solution is 7 = 70 = 4/27, corresponding to a double pole at x^ax = 3/2. 
The argument is as follows. In terms of the scaled radius x = R/Rc, Eqs. (^0]) and 
m imply 



dx _ 1/1 1 

dt 3'jt \x X 



2 



^.(x), (44) 

the last equality defining the function g{x). The form of g{x) is sketched in Figure 
8, where the arrows indicate the fiow of x under the dynamics (?3). From Figure 



8(a) it is clear that for 7 < 70 all drops with x > Xi will asymptotically approach 
the size X2Rc(t), which tends to 00 with t. Therefore it will not be possible to 
satisfy the condition (PSp which imposes the conservation of the order parameter. 
Similarly, Figure 8(c) sEows that for 7 > 70, all points move to the origin and the 
conservation condition again cannot be satisfied. The only possibility is that 7 tends 
to 7o asymptotically from above (it cannot reach 70 in finite time, otherwise all drops 
with X > 3/2 would eventually arrive at x = 3/2 and become stuck, much like the 
case 7 < 7o). 

Evaluating the integral (^Hf ) with 7 = 70 = 4/27 gives the scaling function for 
the droplet size distribution: 

f{x) = const, x^ (3 + x)'^"'^'''^^ ^- — x^ exp ^— — (45) 

for X < 3/2, and f{x) = for x > 3/2. The constant prefactor is fixed by the 
normalisation integral (P5|). 

Lifshitz and Slyozov nave shown that the above scaling solution is obtained for 
generic initial conditions in the limit of small volume fraction v of the minority 
phase. The general-ci form ( ^3[ ) for f{x) has been derived by Yao et al. PS|. Note 



that the method described above only works for d > 2: it is easy to show that 
the constant ad in (|40|) , which sets the time scale for the growth via (j^), vanishes 
linearly with {d — 2) for d —>■ 2, refiecting the singular nature of the Green's function 
for the Laplacian in d = 2. The general expression is ad = (d— l){d — 2)a/4. 



Working in the limit of strictly vanishing v, Rogers and Desai found the same 



scaling form (HH) for d = 2, but with R^. ~ (t/lnt)^/^. For small, but non-zero v, 
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Yao et al. [^] find Rc ~ (t/Jhit'|)^/^. The two results correspond to taking the limit 
f — s> before |^ or after the limit t — oo. 

Many groups have attempted, with varying degrees of success, to extend the 
LSW treatment to non-zero f , either by expanding in v (actually, a/u) [0, 0, or by 



the use of physically motivated approximation schemes and/or numerical methods 



J5|, |3^, p^- When v is of order unity, however, such that both 

pEases are continuous, dnlerent techniques are required, discussion of which will be 
postponed to section 



2.7 BINARY LIQUIDS 

The phase separation of binary liquids is a phenomenon of considerable experimental 
interest. Model B is inappropriate for this system, since it takes no account of the 
transport of the order parameter by hydrodynamic flow. Here we briefly review the 
modifications to model B needed to describe binary liquids. 

The principal new ingredient is 'advection' of the order parameter by the fluid. 
The appropriate modification of (|^) is 

d(l)/dt + V ■ V0 = AVV , (46) 

where v is the (local) fluid velocity, and we have reinstated the transport coefficient 
A. The velocity obeys the Navier-Stokes equation which, with the simplification 
that the fluid is incompressible, reads 



P + (v ■ V)vj = r^VV - Vp - , (47) 

where p is the pressure, Tj the viscosity, and the density p is constant. The final term 
in ( ^7]) arises from the free energy change per unit volume (j)5p that accompanies 
the transport of a fluid region with order parameter over a distance for which the 
change in the chemical potential is 5/i : chemical potential gradients act as a driving 
force on the fluid. 

In the overdamped limit appropriate to most experimental systems, the left side 
of (P7|) can be set to zero. The velocity is then 'slaved to the order parameter' 
p8|, |39|| . The resulting linear equation for v can be solved in Fourier space: 

vM = ^ {-^Kpiy) + Fa(k)) , (48) 

where F = — 0V/i. The incompressibility condition k ■ v(k) = determines the 
pressure. Putting the result into (|48| ) gives, with the summation convention for 
repeated indices, 

i;,(k) = T„^(k)F^(k) 

r..(k) ^ ^ (4.-^1^). (49) 

where T is the 'Oseen' tensor. In real space (for = 3) 



To.,{v) = -^(5^, + '^] . (50) 



Putting everything together gives the equation of motion in real space, 

d(t)/dt = AVV - / dr' [V0(r) ■ r(r - r') ■ V'0(r )] /i(r') , (51) 
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where we recall that /i(r) = 6F/6(j){r). In deriving the final form (0), integration 
by parts (exploiting the transverse property, Taf3(k)k/^ = 0, of the Oseen tensor) was 
used to convert 0(r')V'/i(r') to — /i(r')V'0(r') inside the integral. 

It should be emphasized that, as usual, thermal fiuctuations have been neglected 
in (|5T|). We have previously argued, on rather general grounds, that these are 
negligible at late times, because the coarsening is controlled by a strong coupling 
renormalisation group fixed point (see, in particular, the discussion in section |3). 
For binary liquids, however, a rather subtle situation can arise when the nominally 
dominant coarsening mechanism does not operate. Then thermal fiuctuations do 
contribute. We will enlarge on this below. 

We can use dimensional arguments to estimate the sizes of the two terms on the 
right-hand side of (|5TD. Using /x ~ cr/L, T ~ ^/v^y ^ind V0 ~ 1/L, gives Xa/L^ 
for the first ('diffusive ) term and a/r]L for the second ('advective') term. Advective 
transport of the order parameter therefore dominates over diffusion for L ^ (Ar/)^/^. 
To determine L{t) in this regime, we use the expression for the fluid velocity, 

v(r) = J dr [T(r - r') • V0(r )] /x(r') . (52) 



Using the same dimensional arguments as before, and also v ~ L/t, gives L{t) 
at/r], a result flrst derived by Siggia This result has been conflrmed by experi- 



ments I^TJ and by numerical simulations ^ Because the inertial terms are 



negligible compared to the viscous force here, we will call this the 'viscous hydrody- 
namic' (or just 'viscous') regime. 

Under what conditions is it correct to ignore the 'inertial' terms on the left-hand 
side of (^7|) ? Using dimensional arguments again, we see that these terms are of order 
pL/t^. Comparing this to the driving term 0V/i ~ a/L"^ on the right, (the viscous 
term rjV'^v is of the same order in the viscous regime), and using the result derived 
above, t rjL/a, for this regime, shows that the inertial terms are negligible when 
L -C rf/ap. At sufficiently late times, when this inequality is violated, the inertial 
terms will therefore be important. In this 'inertial' regime, L{t) is determined by 
equating the inertial terms, which scale as pL/t^, to the driving term 0V/i, which 
scales as cr/L^ (and the viscous term is negligible) to give L [aip' / pY^^ . The t^/^ 
growth in the inertial regime was flrst predicted by Furukawa ||45|| . 

To summarise, there are in principle three growth regimes for phase separation 
in binary liquids, after a deep quench, with the growth laws 

Lit) ~ [XatY'"^ , L < {Xr]Y'^ , ('diffusive') , (53) 
~ <yt/rj , {XviY^"^ ^ L ^ rf / pa , ('viscous hydro dynamic') (54) 
~ [at'^/pY^^, L^rj^/pa, ('inertial hydrodynamic') . (55) 



These results basically follow from dimensional analysis. The 'inertial hydrody- 
namic' regime has not, to my knowledge, been observed experimentally and we will 
not consider it further. However, a t^/^ regime has been observed at late times in 
simulations of two-dimensional binary liquids |^4|, ^ 



Siggia fi^ has discussed the physical origin of the linear growth in the 'viscous 
hydrodynamic' regime. He argues that the essential mechanism is the hydrodynamic 
transport of fluid along the interface driven by the surface tension. This mechanism, 
however, can only operate if both phases are continuous. If, by contrast, the minor- 
ity phase consists of independent droplets (which occurs for volume fractions less 
than about 15%), this mechanism tends to make the droplets spherical but does not 
lead to any coarsening (it is easy to show, for example, that the hydrodynamic term 
in (|5T|) vanishes for a single spherical droplet). In the absence of thermal fluctu- 
ations, therefore, the Lifshitz-Slyozov evaporation-condensation mechanism deter- 
mines the growth even beyond the nominal crossover length given above. Thermal 
fluctuations, however, facilitate a second coarsening mechanism, namely droplet co- 
alescence driven by Brownian motion of the droplets. Again, Siggia has given the 
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essential argument. The mobility /i of a droplet of size L is of order l/rjL, so the 
diffusion constant is given by the Einstein relation as Z) = fiksT ~ ksT/riL, where 
is Boltzmann's constant. The time for the droplet to diffuse a distance of or- 
der L (and coalesce with another droplet) is t ~ L'^/D ~ rjL^/kBT, which gives 

L ~ (kBTt/r^y/^ 

The presence of two different mechanisms leading to the same growth exponent 
suggests a 'marginal operator' in the theory, and the Renormalisation Group treat- 
ment of section ^ lends support to this view. The RG approach also shows how, 
in this nominally irrelevant variable (temperature) can affect the late-stage 

growth in a non-trivial way. From a physical point of view, it seems plausible that 
the presence of competing mechanisms will lead to a late-stage morphology that de- 
pends on the ratio of the amplitudes derived for the two mechanisms, i.e. that scal- 
ing functions will depend continuously on this ratio. This could tested by numerical 
simulations, where the transport coefficient A and viscosity rj can be independently 
varied. In real binary liquids, however, these coefficients are related, and the ratio 
of amplitudes for the two mechanisms depends only on the volume fraction v (see 
pO| and the discussion in section |8.4| ). Note that, even without hydrodynamics, 
the scaling functions will depend on v, since the morphology does. The role of the 
Lifshitz-Slyozov mechanism can be enhanced by going to small v, since the growth 
rate due to evaporation-condensation is independent of v for small v. By contrast, 
the coalescence rate increases with v. To see this, we refine the argument given in 
the previous paragraph . 

Let Rhe a typical droplet radius. Then the droplet number density is n ~ v/E?. 
The time for a droplet to diffuse a distance of order its radius is ~ R^/D. The 
volume swept out by the drop in time t (for t > tn) is of order R?t/tR ~ RDt. In a 
'coalescence time' tc the expected number of drops in this volume is of order unity, 
i.e. nRDtc ~ 1, giving tc ~ R^/vD ~ rfR^/vksT, where we used D ~ kBT/rjR in 
the last step. This implies that R grows with time as i? ~ {vkBTt/r]Y^^, a result 
first given by Siggia 



3 TOPOLOGICAL DEFECTS 

The domain walls discussed in the previous section are the simplest form of 'topolog- 
ical defect', and occur in systems described by scalar fields ||^. They are surfaces, 
on which the order parameter vanishes, separating domains of the two equilibrium 
phases. A domain wall is topologically stable: local changes in the order parameter 
can move the wall, but cannot destroy it. For an isolated flat wall, the wall profile 
function is given by the solution of (p!0[) , with the appropriate boundary conditions, 
as discussed in section ^]2| (and sketched in Figure 6). For the curved walls present 
in the phase ordering process, this will still be an approximate solution locally, pro- 
vided the typical radius of curvature L is large compared to the intrinsic width (or 

'core size'), ^, of the walls. (This could be defined from ( [I2D as ^ = [l^"(l)]~^/^, 
say). The same condition, L ^, ensures that typical waTTseparations are large 
compared to their width. 

Let us now generalize the discussion to vector fields. The 'O(n) model' is de- 

scribed by an n-component vector field 0(x, t), with a free energy functional 

that is invariant under global rotations of (j). A suitable generalization of (|I|) is 

F[4 = J d'^x (V0)^ + \/(0)) , (56) 

where (V(/>)^ means Z^iLi Z]a=i('9i0")^ (i-e. a scalar product over both spatial and 
'internal' coordinates), and V{(f)) is 'mexican hat' (or 'wine bottle') potential, such 
as (1 — 0^)^, whose general form is sketched in Figure 9. It is clear that F[(j)] is 
invariant under global rotations of (f) (a continuous symmetry), rather than just the 



16 



inversion symmetry (0 —(f), a discrete symmetry) of the scalar theory. We will 

adopt the convention that V has its minimum for 0^ = 1. 

For non-conserved fields, the simplest dynamics (model A) is a straightforward 
generalization of namely 

d$/dt = VV - dV/d$ . (57) 

For conserved fields (model B), we simply add another (— V^) in front of the right- 
hand side. 

Stable topological defects for vector fields can be generated, in analogy to the 
scalar case, by seeking stationary solutions of (|57D with appropriate boundary con- 
ditions. For the 0{n) theory in d-dimensional space, the requirement that all n 

components of vanish at the defect core defines a surface of dimension d — n (e.g. 
a domain wall is a surface of dimension d — 1: the scalar theory corresponds to 
n = 1). The existence of such defects therefore requires n < d. For n = 2 these 
defects are points ('vortices') for d = 2 or lines ('strings', or 'vortex lines') for d = 3. 
For n = 3, d = 3 they are points ('hedgehogs', or 'monopoles'). The field configura- 
tions for these defects are sketched in Figures 10(a)-(d). Note that the forms shown 
are radially symmetric with respect to the defect core: any configuration obtained 

by a global rotation is also acceptable. For n < d, the field only varies in the n 
dimensions 'orthogonal' to the defect core, and is uniform in the remaining d — n 
dimensions 'parallel' to the core. 

For n < d, the defects are spatially extended. Coarsening occurs by a 'straighten- 
ing out' (or reduction in typical radius of curvature) as sharp features are removed, 
and by the shrinking and disappearance of small domain bubbles or vortex loops. 
These processes reduce the total area of domain walls, or length of vortex line, in 
the system. For point defects {n = d), coarsening occurs by the mutual annihilation 
of defect-antidefect pairs. The antidefect for a vortex ('antivortex') is sketched in 
Figure 10(e). Note that the antivortex in not obtained by simply reversing the di- 
rections of the arrows in 10(b): this would correspond to a global rotation through 
TT. Rather, the vortex and antivortex have different 'topological charges': the fields 
rotates by 27r or — 27r respectively on encircling the defect. By contrast, an anti- 
mo nopole is generated by reversing the arrows in 10(d): the reversed configuration 
cannot be generated by a simple rotation in this case. 

For the radially symmetric defects illustrated in 10(b)-(d), the field has the 

form 0(r) = -r f{r), where f is a unit vector in the radial direction, and /(r) is the 



profile function. Inserting this form into (|57|) , with the time derivative set to zero, 
gives the equation 

^^(^^_(^ (58) 

with boundary conditions /(O) = 0, f{oo) = 1. Of special interest is the approach 
to saturation at large r. Putting /(r) = 1 — e(r) in (|58|), and expanding to first 
order in e, yields 

r^oo. (59) 

This should be contrasted with the exponential approach to saturation (|T2|) for scalar 
fields. A convenient definition of the 'core size' is through / ~ 1 — ^^/r^ for large 
r. This gives ^ = [{n - l)/1/"(l)]i/2 for n > I. 

3.1 DEFECT ENERGETICS 

Consider an isolated, equilibrium defect of the 0{n) model in (i-dimensional space 
(with, of course, n < d). For a radially symmetric defect, 0(r) = /(r) r, the energy 
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per unit 'core volume' (e.g. per unit area for a wall, per unit length for a line, or per 
defect for a point) is, from ( |56D 



E = Sj drr--' f + - {VfY + V{f) , (60) 



i n -I) 2,1 
^ 2r2 ^ 2 

where Sn = 27r"/2/r(n/2) is the surface area of an n-dimensional sphere. 

For scalar fields (n = 1), we have seen (section |2.2| ) that the terms in (V/)^ 
and V{f) contribute equally to the wall energy. For n > 2, the first term in (|60]) 
dominates the other two because, from (^), the three terms in the integrand fall on 
with distance as r"^, and V{f) ~ V"{1){1 — f)"^ ~ respectively as r — »• cxo. 
For n > 2, therefore, the first term gives a divergent integral which has to be cut 
off as the system size Lgyg, i.e. E ~ hi{Lsys/^) for n = 2 and E ~ L^ys^ for ^ > 2. 
Actually, the second and third terms give divergent integrals for n > 6 and n > 4 
respectively, but these are always subdominant compared to the first term. 

The above discussion concerns an isolated defect. In the phase ordering system 
the natural cut-off is not Lgyg but L{t), the characteristic scale beyond which the 
field of a single defect will be screened by the other defects. Of particular interest are 
the dynamics of defect structures much smaller than L{t). These are the analogues 
of the small domains of the scalar system. For d = n = 2, these are vortex-antivortex 
pairs, for d = 3, n = 2 they are vortex rings, while for d = 3 = n they are monopole- 
antimonopole pairs. For such a structure, the pair separation r (for point defects) 
or ring radius r (for a vortex loop) provide the natural cut-off. Including the factor 
^d-n ^Yie volume of defect core, the energy of such a structure is 

E ~ r'^-^ ln(r/0 , d>n = 2, 

~ , d>n>2. (61) 

The derivative with respect to r of this energy provides the driving force, —dE/dr, 
for the collapse of the structure. Dividing by r'^~"' gives the force F acting on a unit 
volume of core (i.e. per unit length for strings, per point for points, etc.): 



F{r) ~ -r"^ , d = n = 2 , 

~ -r"-3 ln(r/^) , d>n = 2 



„n— 3 



d>n>2. (62) 



In order to calculat e th e collapse time we need the analogue of the 'friction 
constant' t] (see section p.3|) for vector fields. This we calculate in the next section. 
Before doing so, we compute the total energy density e for vector fields. This can 
be obtained by putting r ~ L{t) in ([UT|), and dividing by a characteristic volume 
L[tY (since there will typically be of order one defect stucture, with size of order 
L{t), per scale volume L{tY), 

e ~ L{t)-^ In (L(t)/0 , d>n = 2 , 

~ L{t)-^ , d>n>2. (63) 

For scalar systems, of course, e ~ L{t)~^. 

As a caveat to the above discussion, we note that we have explicitly assumed 
that the individual defects possess an approximate radial symmetry on scales small 



compared to L{t). It has been known for some time [^, however, that an isolated 
point defect for d > 3 can lower its energy by having the field uniform (pointing 
'left', say) over most of space , with a narrow 'flux tube' of field in the opposite 
direction (i.e. pointing 'right'). The energy is then linear in the size of the system, 
E ~ Lsys^'^~^, which is smaller than the energy, ~ L'^ys, of the spherically symmetric 
defect, for (i > 3. A defect-antidefect pair with separation r, connected by such a 
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flux tube, has an energy E ~ r^*^"^, which implies an r- independent force for all 
d > 3, in contrast to (§2]). 

How relevant are these considerations in the context of phase-ordering dynamics? 
These single-defect and defect-pair calculations treat the field as completely relaxed 
with respect to the defect cores. If this were true we could estimate the energy 
density for typical defect spacing L{t) as ^'^~^L(t)^~'^ for d > 3. However, the 
smooth variation ('spin waves') of the field between the defects gives a contribution 
to the energy density of (V0)^ ~ L(t)~^, which dominates over the putative defect 
contribution for d > 3. Under these circumstances, we would not expect a strong 
driving force for point defects to adopt the 'fiux tube' configuration, since the energy 
is dominated by spin waves. Rather, our tentative picture is of the point defects 
'riding' on the evolving spin wave structure for d > 3, although this clearly requires 
further work. Note, however, that these concerns are only relevant for d > 3: Eq. 
(|6^) is certainly correct for the physically relevant cases d < 3. 



3.2 DEFECT DYNAMICS 

Here we will consider only nonconserved fields. Using the methods developed in 
section |^, however, it is possible to generalise the results to conserved fields . 
The caveats for d > 3 discussed in the previous subsection also apply here. 

The calculation of the friction constant r) proceeds as in section |2.3| . Consider an 
isolated equilibrium defect, i.e. a vortex for d = n = 2, a monopole for n = d = 3, 
a straight vortex line for n = 2, d = 3 etc. Set up a Cartesian coordinate system 
xi, . . . ,Xd. For extended defects, let the defect occupy the (hyper)-plane defined by 
the last d — n Cartesian coordinates, and move with speed v in the xi direction. 

Then only depends on coordinates Xi, . . . , a;„, and the rate of change of the system 
energy per unit volume of defect core is 

dE/dt = I dxi... dxn iSF/6$) ■ d$/dt 

dxi... dxn {d$/dtf . (64) 

The defect profile has the form 0(xi, . . . ,x„) = /(xi — vt,X2, ■ ■ ■ ,x„), where the 
function / depends on v in general. Putting this into (|64D gives 



dE/dt = J dxi . . .dxn{d(j)/dxif' 

= -{v^/n) f d"r{V(pf = -r]v^ , (65) 



where the function / has been replace by its f = form to lowest order in w, and 
7] is the friction constant per unit core volume. The final expression follows from 
symmetry. It follows that r] is (up to constants) equal to the defect energy per unit 
core volume. In particular, it diverges with the system size for n > 2. For a small 
defect structure of size r, we expect the divergence to be effectively cut off at r. 
This gives a scale-dependent friction constant, 

r/(r) ~ r"-2 ln(r/0 , d>n = 2 , 

~ r"-2 , d>n>2. (66) 

Invoking the scaling hypothesis, we can now determine the growth laws for non- 
conserved vector systems. Eqs. (B3) and (|66D give the typical force and friction 
constant per unit core volume as F(L) ancrrf(L). Then a typical velocity is f ~ 
dL/dt ~ F{L)/t]{L), which can be integrated to give, asymptotically, 

L{t) ~ {t/lnty/\ d = n = 2, 

~ t^^"^ , otherwise . (67) 
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The result for n 
by Yurke et al. 



= d = 2 was derived by Pargellis et al. [5^, and checked numerically 
The method used here follows their approach |^2|]. The key 
concept of a scale-dependent friction constant has been discussed by a number of 
authors 
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A detailed analysis of monopole-antimonopole annihilation, in the 
context oTnematic liquid crystals, has been given by Pismen and Rubinstein p4 |. 

A more general and powerful method to derive growth laws, valid for both con- 
served and nonconserved systems, is the subject of section |^. The results agree with 
the intuitive arguments presented so far, with the possible exception of the case 
d = 2, for which evidence suggestive of scaling violations will be presented. 



n 



3.3 POROD'S LAW 

The presence of topological defects, seeded by the initial conditions, in the system 
undergoing phase ordering has an important effect on the 'short-distance' form of the 
pair correlation function C(r, t), and therefore on the 'large-momentum' form of the 
structure factor ^(k, t). To see why this is so, we note that, according to the scahng 

hypothesis, we would expect a typical field gradient to be of order |V0| ~ 1/L. At a 
distance r from a defect core, however, with ^ ^ r ^ L, the field gradient is much 

larger, of order 1/r (for a vector field), because = f implies (V*/))^ = (n — l)/r^. 
Note that we require r ^ for the field to be saturated, and r <^ L for the defect 
field to be largely unaffected by other defects fwhich are typically a distance L 
away). This gives a meaning to 'short' distances <^ r <^ L), and 'large momenta' 
<^ k -C ^~^). The large field gradients near defects leads to a non-analytic 
behaviour at x = of the scaling function f{x) for pair correlations. 

We start by considering scalar fields. Consider two points x and x + r, with 
C -C r <C L. The product 0(x) 0(x -|- r) will be —1 if a wall passes between them, 
and +1 if there is no wall. Since r <^ L, the probability to find more than one 
wall can be neglected. The calculation amounts to finding the probability that a 
randomly placed rod of length r cuts a domain wall. This probability is of order 
r/L, so we estimate 

C(r,t) ~ (-1) X (r/L) + (+l) X (1-r/L) 

= 1 - 2r/L , r < L . (68) 

The factor 2 in this result should not be taken seriously. 

The important result is that ( |68D is non- analytic in r at r = 0, since it is linear 
in r = |r|. Technically, of course, this form breaks down inside the core region, when 
r < ^. We are interested, however, in the scaling limit defined by r ^ ^, L ^ ^, 
with X = r/L arbitrary. The nonanalyticity is really in the scaling variable x. 

The nonanalytic form (|68D implies a power-law tail in the structure factor, which 
can be obtained from (|6^) py simple power-counting: 

^(k,t)~-^^, (69) 

a result known universally as 'Porod's law'. It was first written down in the general 
context of scattering from two-phase media [^. Again, one requires fc^ <C 1 for 
the scaling regime. Although the fc-dependence of ( |69D is what is usually referred to 
as Porod's law, the L-dependence is equally interesting. The factor 1/L is simply 
(up to constants) the total area of domain wall per unit volume, a fact appreciated 
by Porod, who proposed structure factor measurements as a technique to determine 



the area of interface in a two-phase medium |]22[. On reflection, the factor 1/L 
is not so surprising. For kL ^ 1, the scattering function is probing structure on 
scales much shorter than the typical interwall spacing or radius of curvature. In 
this regime we would expect the structure factor to scale as the total wall area, 
since each element of wall with linear dimension large compared to 1/k contributes 
essentially independently to the structure factor. 
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This observation provides the clue to how to generalize (|69| ) to vector (and other) 
fields p5| , The idea is that, for kL ^ 1, the structure factor should scale as the 
total volume of defect core. Since the dimension of the defects is d — n, the amount 
of defect per unit volume scales as L~". Extracting this factor from the general 
scaling form (^ yields 



S{k,t) 



1 



kL > 1 



(70) 



for the 0{n) theory, a 'generalized Porod's law'. 

Equation (^) was first derived from approximate treatments of the equation of 



motion 



for nonconserved fields 
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In these derivations, however, 

the key role of topological defects was far from transparent. The above heuristic 
derivation suggests that the result is in fact very general (e.g., it should hold equally 
well for conserved fields), with extensions beyond simple 0{n) models. The appro- 
priate techniques, which also enable the amplitude of the tail to be determined, were 
developed by Bray and Humayun [pH], and will be discussed in detail in section 10. 



3.4 NEMATIC LIQUID CRYSTALS 

Liquid crystals have been a fertile area for recent experimental work on the kinetics of 



phase ordering, largely due to the efforts of Yurke and coworkers ||21|, ^ 0, ^ . 
Here we will concentrate on the simplest liquid crystal phase, the nematic. In a 
simple picture which captures the orientational degrees of freedom of the nematic. 



the liquid crystal can be thought of as consisting of rod-like molecules [£4[ which 
have a preferential alignment with the 'director' n in the ordered phase. Because 
the molecules have a head-tail symmetry, however, the free energy is invariant under 
the local transformation n — n. As a result, the nematic is usually described by 
a tensor order parameter Q that is invariant under this local transformation, and 
has in the ordered phase the representation 

Qab = S {naUb - i 6ab) ■ (71) 

Note that Q is a traceless, symmetric tensor. The simplest free energy functional is 



one that is invariant under global rotations of the field n(r). It has the form ||64 



m=J d'x (1 Tr {VQr + "-TTQ'-^TrQ' + l (Tr Q^)^) , (72) 

where we have retained only terms up to order Q^, and the lowest order term 
involving spatial gradients. (Note that TrQ^ = (1/2) (TrQ^)^ for a 3 x 3 traceless, 
symmetric tensor, so we don't include this term separately). The presence of the 
cubic term, allowed by symmetry, leads to a first-order phase transition in mean-field 
theory p4|. In experimental systems the transition is weakly first order. 

The iiFee energy functional (^) is an idealisation of real nematics, in the spirit of 
the Lebwohl-Lasher lattice Hamiltonian = —JJ2<i,j>{^i ■ ^jY^ where iij is the 
local director at site i. Both models are invariant under global rotations of n, as well 
local inversions, n — n. The gradient terms in ( [72| ) can be written, using ([H]), 
as (ignoring constants) Y.i.aidn'^ / dxiY , i.e. as (Vn), which is a continuum version 
of Hll. The overarrow on n here indicates that, in these models, the spatial and 
'internal' spaces can be considered as distinct (much as in the 0{n) model, where 
n and d can be different). In real nematics, however, these spaces are coupled. An 
appropriate gradient energy density is the Frank energy 

Ef = Ki{V ■ n)2 + (n • V X n)2 + [n x (V x n)]^ , (73) 

where the Frank constants Ki, K2 and are associated with 'splay', 'twist' and 
'bend' of the director p|[. The isotropic models discussed above correspond to the 
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case Ki = K2 = K3, the much-used 'equal-constant approximation'. We will limit 
our considerations exclusively to this case. 

In a similar spirit, we adopt the simplest possible dynamics, namely the purely 
relaxational dynamics of model A. This captures correctly the non-conserved nature 
of the dynamics, but ignores possible complications due to hydrodynamic inter- 
actions. Recent work, comparing experimental results with simulations based on 
relaxational dynamics, provide some justification for this approach [Q. The equa- 
tion of motion is dQ/dt = -6G[Q]/6Q, where G[Q] = F[Q] - J d^x X{x.)Q{x.) and A 
is a Lagrange multiplier introduced to maintain the condition Tt Q = 0. Imposing 
the constraint to eliminate A gives 

T-r2/^ „^ , r rp,, ^2 



V'Q-rQ + w[Q' --ITiQ'j-uQTiQ', (74) 
whe re / is the unit tensor. This equation will be discussed in more detail in section 



5.4 



Due to the extra local symmetry (compared to 0(3) models) under n — >■ — n, 
nematic liquid crystals support a number of defect types |4^. In the present context, 
the most important are string defects, or 'disclinations', in which the director rotates 
through ±7r on encircling the string, as sketched in Figure 11. These ±l/2-strings 
are topologically stable, in contrast to ±1 string configurations (in nematics and the 
0(3) model) which can be relaxed by smoothly canting the order parameter towards 
the string axis ('escape in the third dimension'). The presence of string defects 
(which have been observed by Yurke's group ||61|), makes the nematic more akin to 



the 0(2) model than the 0(3) model as far as its ordering kinetics are concerned. In 
particular, if we make the natural assumption that the total string length decreases 
as L ^, a Porod law of the form (|70D with n = 2 is predicted. Scattering data were 
originally interpreted as being consistent with n = 3-like behaviour, i.e. an effective 
exponent d + n = 6.0 ± 0.3 in ( |70D was observed but it is not clear whether the 

16611 . Numerical simulations 



appropriate region of the structure factor tail was ntted 
of a 'soft-spin' version of the Lebwohl-Lasher model ||67|| are f ully consistent with a 
tail exponent of 5 [BH]. We will return to this point in section IST 



4 EXACTLY SOLUBLE MODELS 

There are few exactly solved models of phase ordering dynamics and, unfortunately, 
these models are quite far from describing systems of physical interest. However, 
the models are not without interest, as some qualitative features survive in more 
physically relevant models. In particular, such models are the only cases in which 
the hypothesized scaling property has been explicitly established. 

We begin by discussing phase ordering of a vector field in the limit that the 
number of vector components of the field, n, tends to infinity. This limit has been 
studied, mostly for nonconserved fields, by a large number of authors [T^, [H^, 
[701 , [ni . In principle, the solution is the starting point for a systematic treatment m 
powers of 1/n. In practice, the calculation of the 0(1/?!) terms is technically difficult 
|7T| . Moreover, some important physics is lost in this limit. In particular, there 



are no topological defects, since clearly n > d + I for any d as n ^ 00. As a 
consequence, Porod's law (|70|) , for example, is not found. It turns out, however, 
that similar techniques can be applied for any n after a preliminary transformation 

from the physical order parameter field to a suitably chosen 'auxiliary field' rh. 
This is discussed in section ^. The topological defects are incorporated through the 

functional dependence of on rh, and Porod's law is recovered. 

4.1 THE LARGE-n LIMIT: NONCONSERVED FIELDS 

Although not strictly necessary, it is convenient to choose in (|57D the familiar '0^' 



potential, in the form V{(f)) = (n — 0^)^/4n, where the explicit n-dependence is for 
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later convenience in taking the limit n oo. With this potential, (0) becomes 

= vV+0- -(0^)0 • (75) 
n 

The simplest way to take the limit is to recognize that, for n oo, can be 

replaced by its average, to give 

d(j)/dt = V'^(f) + a{t)(j) (76) 
ait) = 1 - (02) , (77) 

where now stands for (any) one of the components of 0. Eq. ( [76|) can alternatively 
be derived by standard diagrammatic techniques |2^. Eq. (^) can be solved exactly 
for arbitrary time t after the quench. However, we are mainly interested in late times 
(i.e. the scaling regime), when the solution simplifies. After Fourier transformation, 
the formal solution of (|76D is 

0k(t) = 0k(O) exp[-A;2t + 6(t)] , (78) 

b{t) = f'dt'a{t'), (79) 



givmg 

a{t) = db/dt = 1 - A ^ exi)[-2kH + 2b{t)] , (80) 

k 

where (^) has been used to eliminate the initial condition. Since we shall find a 
posteriori that a(t) -C 1 at late times, the left side of (pO| ) is negligible for t — >• oo. 
Using J2kexp{-2kH) = [Mty^'^ gives b{t) (rf/4) ln(t/to), where 

to = A2/787r . (81) 
Therefore, a{t) d/At for t oo, and the solution of (|78D, valid at late times, is 

0k(t) = 0k(O) (t/to)"/' exp(-Pt) . (82) 

Using (1) once more, we obtain the structure factor, and its Fourier transform, the 
pair correlation function as, 

S{k,t) = {8nty/^ exp{-2kH) , (83) 
C{r,t) = exp(-rV8t) . (84) 

These results exhibit the expected scaling forms (|^, with length scale L{t) oc t^/^. 
Note that the structure factor has a gaussian tail, in contrast to the power-law tail 
(froj) found in systems with n < d. It might be hoped, however, that the large-n 
forms (^) and (|8^ ) would be qualitatively correct in systems with no topological 
defects, i.e. ioi n > d + 1. These cases will be discussed in section ||. 



4.2 TWO-TIME CORRELATIONS 



Within the large-n solution, we can also calculate two-time correlations to test the 
scaling form (|g). It turns out (although this becomes apparent only at 0(l/n)) 
pO| , [7l| that there is a new, non-trivial, exponent associated with the limit when 
the two times are well separated [[72l . 
From (^2]) it follows immediately that 



5(k,t,t') 
Cir,t,t') 



(0k(t)0-k(t')) = 
(0(r,t)0(O,t')) 



in 



{tt 



Att' 



r\l/2]d/2 



exp[-k^{t + t')] 



d/A 



exp 



4(t + 1' 



(85) 
(86) 
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Eq. ( p6D indeed has the expected form (||). In the hmit t ^ t', ( p6D becomes 



C(r,t,t') = (4t7t)'^/%xp(-rV4t) , (87) 
= {L'/Lfh{r/L), (88) 

where the last equation defines the exponent A through the dependence on the later 
time t. Clearly, A = d/2 for n = oo. When the Oil/n) correction is included, 



however, an entirely non-trivial result is obtained P0| , ffl 

It is interesting to consider the special case where the earlier time t' is zero. Then 
C(r,t,0) is just the correlation with the initial condition. This quantity is often 
studied in numerical simulations as a convenient way to determine the exponent A. 



Within the large-n solution, Eqs. (0) and (|8lD give, in Fourier and real space 



S(k,t,0) = [87i{ttoY/Y^^ exp{-kH) , (89) 
C(r,t,0) = (4to/t)'^/%xp(-rV4t) . (90) 

This is just what one gets by replacing t' by to in ( pS] ) and ( ^6] ) (with to playing the 
role of a short-time cut-off), and then neglecting tg comparea to t. 

A related function is the response to the initial condition, defined by 

G{k,t) = {dMt)/dMO)) ■ (91) 
Within the large-n solution, (^2p gives immediately 

G(k,t) = {t/toY/* exp{-kH) . (92) 

Comparing (|89[) and (^), and using (|8lD once more gives the relation. 

5(k,t,0) = AG(k,t) . (93) 

In fact, this is an exact result, valid beyond the large-n limit, as may be proved 
easily using integration by parts on the gaussian distribution for {0k(O)}. The 
general scaling form for (^(k, t), 

G{k,t) = gnikL) , (94) 

defines a new exponent A, equal to d/2 for n = oo. Since, however, the correlation 

with the initial condition has the scaling form C(r,t,0) = f{r/L), the identity 
(0) gives immediately ||7^ 

X = d-X. (95) 

(The symbol A is also used for the transport coefficient in systems with conserved 
dynamics. This should not be a source of confusion, as the meaning will be clear 
from the context). 

Before leaving this section, it is interesting to consider to what extent the results 
depend on the specific form (^) chosen for the correlator of the initial conditions. 
Let us replace the right-hand side of (^ by a function A(|x — x'|), with Fourier 
transform A{k). Then A{k) will appear inside the sum over k in Eq. (|80|) . The 

dominant k values in the sum, however, are of order t~^/^, so for late times we can 
replace A (A;) by A(0), provided the latter exists. This means that universal results 
are obtained when only short-range spatial correlations are present at t = 0. For 
sufficiently long-range correlations however, such that A{k) diverges for — ^ 0, new 
universality classes are obtained. We shall return to the role of initial conditions, 
from a more general perspective, in section |^. 
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4.3 THE LARGE-n LIMIT: CONSERVED FIELDS 

For conserved fields, the calculation proceeds as before, but with an extra (— V^) on 
the right-hand side of the equation of motion. Making as before, the replacement 

(fp' /n — > {(jp) for n —>■ oo, where (j) is (any) one component of (p, one obtains 

d(j)/dt = -VV - ci{t) , (96) 
with a{t) still given by (ffTl) . Transforming to Fourier space, the solution is 

Mt) = 0k(O) exp[-kH + eb{t)] , (97) 
The function b(t), defined as in ([79|), satisfies the equation 

a{t) = db/dt = 1 - A ^ exp[-2kH + 2k%{t)] . (98) 

k 



This equation was solved by Coniglio and Zannetti fig] , by first expressing the sum 
over k as a parabolic cylinder function, then taking the large-t limit. Here we will 
take the large-t limit from the outset, and recognize that the sum can then be 
evaluated using steepest descents. Just as for the nonconserved case, we can show a 
posteriori that db/dt <^ 1 at late times, so that this term can be dropped from (0). 
After the change of variable k = [6(t)/t]^/^x, we obtain 

1 = /\Cd{b/tYI^ / x'^-^dx exp[2/3(x2 - x% , (99) 
Jo 

where Cd is an uninteresting constant, and 

/?(t) = b'^{t)/t . (100) 

Provided j3{t) oo for t ^ oo (which can be verified a posteriori), the integral 
on the right of (|99|) can be evaluated by steepest descents. Including the gaussian 
fluctuations around the maximum of the integrand at x = 1/ ^J2 gives 

1 = const A p-^/^ {P/tf^ exp(/3/2) , (101) 
with asymptotic solution 

/3~(d/2)lnt, t-*oo, (102) 

justifying the use of the steepest descents method for large t. Putting this result 
into (|97|) gives the final result for the structure factor [|l6l 

5(k,t) - fid/mk/km) (103) 

(\){x) = l-{l-x^f. (105) 

Here km{t) is the position of the maximum in S{\<.,t). A slightly more careful 

treatment (retaining the leading subdominant term in ( |102|) ), gives an additional 

logarithmic prefactor, of order (lnt)^^~°')/^, in (103), such that (asymptotically in 
time) Ek5Xk,t) = 1. 

Eq. (|l03D is interesting because, in contrast to the nonconserved result (H), 
it does not have the conventional scaling form. Rather it exhibits 'multiscaling' 
T6| . In particular there are two, logarithmically different, length scales, k^ and 



L = t^/^. For simple scaling, these two scales would be the same. Furthermore, 
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for fixed scaling variable, which can be written as k/km, the structure factor would 
vary as L{tY, with a j)re/actor depending on the scaling variable. In the multiscaling 
form (|103|) , for fixed scaling variable, S{h,t) ~ L'^'t'C^/^^) ^ i.e. the ex^^onent depends 
continuously on the scaling variable. 

After the discovery of multiscaling in the n — > oo limit, some effort was devoted 
to looking for similar phenomena at finite n, notably for scalar systems |7^, but also 



for n = 2 [7^, |76| and n = 3 [f7|. However, no evidence was found for any departure 
from simple scaling for any finite n. At the same time. Bray and Humayun showed, 
within the context of an approximate calculation based on an idea of Mazenko, 
that simple scaling is recovered asymptotically for any finite n [|T9[. This result is 
discussed in detail in section ^. 

4.4 THE ONE-DIMENSIONAL ISING MODEL 

An exceptionally simple system that can be solved exactly [O is the Ising model 
in one dimension with Glauber dynamics. It is defined by the Lllauber equation for 
the spin probability weight: 

+ 1:P(S......-S......S„;«) fi±™?5.)(106) 



where /3 = 1/T, hi = J{Si-i + Si+i) is the local field at site i, and periodic boundary 
conditions , Sj+ N = Si, have been adopted. 

From ( |106| ), it is straightforward to derive the equation of motion for the pair 
correlation function, Cjj(t) = {Si(t)Sj{t)) , where the brackets indicate an average 
over the distribution P. After averaging also over the initial conditions, Cij depends 
only on the difference r = |i — j| if the ensemble of initial conditions is invariant 
under translations. Then one obtains 



{d/dt)C{r,t) = C{r + l,t) -2C{r,t) + C{r -l,t) , r^O. (107) 

For r = one has trivially C(0, t) = 1 for all t . To solve for C in the scaling limit, it 
is simplest to take the continuum limit, when ( p.07| ) reduces to the diffusion equation, 
dC/dt = d'^C/dr'^, with constraint C(0, t) = 1. A scaling solution obviously requires 
L(t) = t^/^. Inserting C{r,t) = /(r/t^/^) in the diffusion equation gives /" = 
— (x/2)/', which can be integrated with boundary conditions /(O) = 1, /(oo) = 
to give f{x) = erfc {x/2), where erfc is the complementary error function. Thus the 
scaling solution is 

C(r,t) = erfc (r/2t^/2) . (108) 

In particular, the solution exhibits the expected Porod regime, C = 1 — r/(7ity^'^ + 
0{r^/t^^'^) at short distance. A more complete discussion can be found in 

The scalar TDGL equation (^ is also soluble in one dimension, in the sense that 
the scaling functions can be exactly calculated [^. When L{t) ^ ^, neighboring 
domain walls interact only weakly, with a force of order exp(—L/^), leading to a 
logarithmic growth law, L ~ ^Int. Moreover, in the limit L/^ — > oo, the closest 
pair of domain walls interact strongly compared to other pairs, so that the other 
walls can be treated as stationary while the closest pair annihilate. This leads to 



a simple recursion for the domain size distribution, with a scaling solution It 
is interesting that the fraction of the line which has never been traversed by a wall 
decays with a non-trivial power of the mean domain size |81|. A similar phenomenon 
(but with a different power) occurs for Glauber dynamics 
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4.5 THE ONE-DIMENSIONAL XY MODEL 

As our final example of a soluble model, we consider the case d = 1, n = 2 with 



non-conserved order parameter. The solution, first given by Newman et al is 

interesting for the 'anomalous' growth law obtained, L(t) ~ t^^^. Here we shall give 
a more detailed discussion than appears in [^, emphasizing the scaling violations 
exhibited by, in particular, the two-time correlation function. In section 0, we shall 
present a general technique, developed with A. D. Rutenberg ||83|], for determining 
growth laws for phase ordering systems. The scaling form (j^) plays an important 
role in the derivation. For the d = 1, n = 2 model, however, our method fails to 
predict the correct t^^^ growth. The reason is precisely the unconventional form (i.e. 
different from (|^)) of the two-time correlation function for this system. 

It is simplest to work with 'fixed length' fields, i.e. 0^ = 1, with Hamiltonian 

F = (1/2) / dx{d(j)/ dxY . The constraint can be eliminated by the representation 

= (cos^, sin6'), where 6 is the phase angle, to give F = (1/2) / dx{dO/dxY. The 
'model A' equation of motion, 89 /dt = —6F/S6, becomes 

de/dt = d^e/dx^ , (109) 

i.e. a simple diffusion equation for the phase. In general dimensions, it is difficult 
to include vortices, which are singularities in the phase field, in a simple way. Such 
singularities, however, are absent for d = 1. 

Eq. ( |1U9D has to be supplemented by suitable initial conditions. It is convenient 
to choose the probability distribution for 9{r, 0) to be gaussian: in Fourier space 

Pmm) oc exp[- E(/5fc/2)^fc(0)^^-fc(0)] . (110) 

k 

Then the real-space correlation function at t = is readily evaluated using the 
gaussian property of the {6'fc(0)}: 

C{r,Q) = (cos[0(r,O) -^(0,0)]) 

= exp{-(l/2)([%,0)-e(0,0)n} 

= exp{-^(l-cosA;r)//3fc} . (Ill) 

k 

The choice Pk = (^/2)A;^ yields C(r, 0) = exp(— |r|/^), appropriate to a quench from 
a disordered state with correlation length ^. 

The general two-time correlation function can be calculated in t he s ame way 
|7|. Using (3k = (^/2)fc2, and the solution Okit) = OkiO) exp{-kH) of (|109D, gives 



C(r, ti,t2) = (cos[^(r,ti)-e(0,t2)]) 

= exp{-il/2)mr,t,)-9{0,t,)r)} 

= exp {[exp{-k'^ti) - exp{-kH2)]'^ 

+2(1 -cosA;r)exp{-A;2(ti + t2)}}) • (112) 

Since we shall find that r is scaled by (^1+^2)^'^^? we can take fcr <^ 1 in the summand 
for the r-values of interest, i.e. we can replace (1 — cosfcr) by {krY/2. Evaluation 
of the sums then gives 



C{r,ti,t2) = exp 



2{h + hy/^ 



+ 2(ti + t2)'/'-(2ti)^/'-(2t2)'/' 



(113) 
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For the special case ti = ^2 = t, ( |113| ) reduces to C(r, t, t) = exp(— r^/2^(27rt)^/^), 
which has the standard scaling form (|^), with growth law L(t) ~ t^^^. This growth 
law is unusual: we shall show in section ^ that the generic form for nonconserved 
fields is L{t) ~ t^/^ just as in the large-n result (|8^) . Another, related, feature of 
( | 113|) is the explicit appearance of ^, the correlatioiTlength for the initial condition. 
I'he large- solution, fo r ex ample, becomes independent of ^ for L{t) ^ ^. The 
most striking feature of ( |113|) , however, is the breakdown of the scaling form (|8[) for 
the two-time correlations! It is this feature that invalidates the derivation of the 
result L(t) ~ t^/^ given in section |^. It is possible that a similar anomalous scaling 
is present in the conserved d = 1 AY model, for which simulation results [^, 77[ 

suggest L{t) ~ t^^^^ instead of the t^^^ growth derived in section |^ assuming simple 
scaling for two-time correlations. Unfortunately, an exact solution for the conserved 

case is non-trivial. 

The explicit dependence of ( |113| ) on C, suggests an unusual sensitivity to the initial 
conditions in this system. A striking manifestation of this is obtained by choosing 
initial conditions with a non-exponential decay of co rrelations. For example, choos- 
ing Pk ~ 1^1" for small \k\ in ( JllOD gives, via ( 111| ), C{r,0) ~ exp(— const |r| 



a- 



for large |r|, provided 1 < a < 3. The calculation of the pair correlation function 
is again straightforward. For example, the equal-time function has a scaling form 
given by C{r,t) = exp(— const r^/t^^"")/^), implying a growth law L{t) ~ 
but the two-time correlation still does not scale properly. 

An especially interesting case is a = 1, which generates power-law spatial corre- 
lations in the initial condition. Thus we choose = 1^1/7 for 1^1 ^ ^5 and /^^ = oo 
for \k\ > A, where A is an ultraviolet cut-off. Then the initial-condition correlator 

has the form C(r, 0) ~ {Ar)~"'^'^ for Ar 3> 1. The general two-time correlation 
function now has the conventional scaling form (^), with L{t) ~ t^/^. Its form is 

C{r,h,h) = f (—-^] (jT^^Y^'^ , (114) 



where f{x) is the equal-time correlation function. In particular, for t2 ^ ti this gives 

C{r,ti,t2) ~ f{r/V^2) (4ti/t2)^/^'', so the exponent A defined by (|) is -f/2n for this 
model. Also the large distance behaviour of the equal-time correlation function is 
f{x) ~ x"'^/'^, exhibiting the same power-law decay as the initial condition. These 
two results are in complete agreement with the general treatment MM of initial 



conditions with power-law correlations given in section |8.3.1 



5 APPROXIMATE THEORIES FOR SCALING 
FUNCTIONS 

While the determination of growth laws (i.e . th e form of L{t)) has proved possible 
using fairly simple arguments (as in section p.5|) , which can be made precise by the 
use of exact relations between correlation functions (section |^), or Renormalization 
Group methods (section ^), the calculation of scaling functions, e.g. the pair cor- 
relation scaling function f{x) (see Eq. (|^)), has been a long-standing challenge. In 
the previous section we have shown that this function can be calculated exactly in a 
number of soluble models. With the exception of the 1 — d Glauber model, however, 
these models lack the topological defects that play such an important role in realistic 



models. In particular, these defects are responsible for the power-law tail ([70|) in 
the structure factor. 

In this section we will review some of the approximate theories that have been 
put forward for the scaling function f{x) of the pair correlation function. The most 
successful by far are theories for nonconserved fields. Even these, however, are not 
quite as good as has been believed, as we shall show. We shall propose a new 
approach which can in principle lead to systematically improvable calculations of 
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scaling functions for nonconserved fields. For conserved fields the theory is in a less 
satisfactory state. We shall try to give some indication of why this is so. Finally we 
emphasize that the discussion is limited throughout to the late-stage scaling regime. 



5.1 NONCONSERVED FIELDS 

A number of approximate scaling functions have been proposed for non-conserved 
fields, but in my view none of them is completely satisfactory. The most physically 
appealing approach for scalar fields is that of Ohta, Jasnow and Kawasaki (OJK) 
which starts from the Allen-Cahn equation (ITB|) for the interfaces. Below we 



wm review the OJK method, as well as an earlier approach by Kawasaki, Yalabik 
and Gunton (KYG) and more recent work by Mazenko |S^, Finally we 
discuss in detail a new approach which has the virtue that it can, in principle, 
be systematically improved. 



5.1.1 The OJK Theory 

A common theme, introduced by Ohta, Jasnow, and Kawasaki (OJK), in the 
approximate theories of scaling functions is the replacement of the physical field 
0(x, t), which is ±1 everywhere except at domain walls, where it varies rapidly, by 
an auxiliary field m(x, t), which varies smoothly in space. This is achieved by using 
a non- linear function (j){m) with a 'sigmoid' shape (such as tanhm). In the OJK 
theory, the dynamics of the domain walls themselves, defined by the zeros of m, are 
considered. The normal velocity of a point on the interface is given by the Allen- 
Cahn equation (9), w = —K = — V-n, where K is the curvature, and n = Vm/| Vm| 
is a unit vector normal to the wall. This gives 

f = {— V^m + nanfoVaVfom}/! Vm| . (US) 

In a frame of reference comoving with the interface, 

dm/dt = = dm/dt + v ■ Vm . (116) 

But since v is parallel to Vm (and defined in the same direction), v.Vm = f |Vm| 
so 

1 dm 

v = -T=^^ . 117 
|Vm| dt ^ ' 

Eliminating v between ( [115| ) and ( |117[ ) gives the OJK equation 

dm/dt = V^m — naUbV a^bf^ ■ (118) 

Since n = Vm/|Vm|, this equation is non-linear. To make further progress, OJK 
made the simplifying approximation of replacing UaUb by its spherical average Sab/d, 
obtaining the simple diffusion equation 

dm/dt = DV'^m , (119) 

with diffusion constant D = {d — l)/d. 

Providing there are no long-range correlations present, we do not expect the form 
of the random initial conditions to play an important role in the late-stage scaling. 
A convenient choice is a gaussian distribution for the field m(x, 0), with mean zero 
and correlator 

(m(x, 0)m(x', 0)) = A(5(x - x') . (120) 
Then the linearity of ( 119 ) ensures that the field m(x, t) has a gaussian distribution 



at all times. Solving ( |119| ), and averaging over initial conditions using ( |120| ) gives 
the equal-time correlation function 



^2 



(m(l)m(2)) = -TT^r exp , (121) 
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where '1' and '2' represent space points separated by r. Of special relevance in what 
follows is the normalized correlator 



7(12) 



(m(l)m(2)) 



(m(l)2)i/2(m(2)2)i/2 



exp 



(122) 



8Dt J 

and will be given explic- 



The generalization to different times is straightforward 
itly below. 

To calculate the pair correlation function of the original field (p, we need to know 
the joint probability distribution for m(l) and m(2). For a gaussian field this can 
be expressed in terms of the second moments of m: 



P(m(l),m(2)) = exp 



1 



2(1-7^ 



m(l)2 m(2)2 m(l)m(2) 

+ ^ 27- 



Sail] 



So{2) 



5o(l)5o(2), 



where 7 = 7(12), and 
Soil) = (m(l)2), 5o(2) = (m(2)2), 



(123) 



N = (27r)-i[(l - 7')So(l)So(2)]-i/2 . (124) 



Note that (|123|) is a general expression for the joint probability distribution of a 
gaussian field, with 7 defined by the first part of (|l22|) . Now '1' and '2' represent ar- 
bitrary space-time points. For the special case where m obeys the diffusion equation 
(HI), 7 is given by 



7 



d/i 



.{tl+t2)\ 



exp 




(125) 



a simple generalization of 

The pair correlation function is given by C(r,t) = (0 (m(l)) (m(2))). In the 
scaling regime, one can replace the function 0(m) by sgn(m), because the walls 
occupy a negligible volume fraction. In a compact notation. 



C(12) = (sgnm(l)sgnm(2)) = (2/77) sin-^(7) 



(126) 



The g auss ian av erag e over the field m required in ( |126|) is standard (see, e.g. 
Eqs. ( |122| ) and (|126|) define the 'OJK scaling function for equal-time pair corre 
tions. Note that (apart from the trivial dependence through D) it is independent 
of the spatial dimension d. We will present arguments that it becomes exact in the 
large-c? limit. The OJK function fits experiment and simulation data very well. As 
an example, we show the function f{x) for the d = 2 scalar theory in Figure 12. 

The general two-time correlation function is especially interesting in the limit 
ti ^ t2 that defines (see, e.g., (P5|)) the exponent A. Since 7 -C 1 in this limit, ( |126D 

can be linearised in 7 to give C{r,ti,t2) ~ (4^1/^2)'^'^'^ exp(— r2/4Dt2), i-e. A = d/2 
within the OJK approximation. 



5.1.2 The KYG Method 

An earlier approach, due to Kawasaki, Yalabik and Gunton (KYG) |S^, building 
on still earlier work of Suzuki [^, was based on an approximate resummation of 
the direct perturbation series in the non-linearity, for the quartic potential V{(j)) = 
(1/4) (1 — 0^)^. The equation of motion (0) for this potential is 

d(j)/dt = + - , (127) 

with g = 1. The basic idea is treat g as small, expand in powers of g, extract the 
leading asymptotic (in t) behaviour of each term in the series, and set g = 1 at the 
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end. However, an uncontrolled approximation is made in simplifying the momentum 
dependence of each term (the expansion is performed in Fourier space) . After setting 
g = 1, the final result can be expressed in terms of the mapping 

0(m) = m/(l + m2)i/2 . (128) 

It is found that m obeys the equation 

dm/dt = V'^m + m , (129) 

instead of (|119|) , which gives an exponential growth superimposed on the diffusion. 
After the rep lacement 0(m) — sgn(m), how ever , this drops out: the OJK scaling 
function (|126|) is recovered, with 7 given by ( |122| ) (but with D = I). 



The nature of the approximation involved can be clarified by putting ( |128|) into 



|127|) (with g = 1) to derive the exact equation satisfied by m: 



dm „9 m(Vm)^ 

= V2m + m-3-^ L. 130 

ot 1 + 



In contrast to a claim made in |]86[ , there is no reason to neglect the final term. On a 
physical level, the fact that this approach gives the correct growth law, L{t) ~ t^/^, 
seems to be fortuitous (see the discussion in section 0). In particular, the crucial 
role of the interfacial curvature in driving the growth is not readily apparent in this 
method. By contrast the OJK approach, while giving the same final result, clearly 
contains the correct physics. 

Despite its shortcomings, the KYG method has the virtue that it can be readily 
extended to vector fields |5^, . Eq. (|129|) is again obtained, but with m replaced 

by a vector auxiliary field rh, with = m/(l + m^)^/^. At late times, — *• m, a 
unit vector, almost everywhere and C(12) = (m(l) • m(2)). Taking gaussian initial 
conditions for m, the resulting scaling function is with 7 again given by 
(but with D = l), 



C(12) = 

ZTT 



B 



n + 1 1 



1 n + 2 



f Ij.j;— ;7^l . (131 



where B{x, y) is the beta function and F{a, b; c; z) the hypergeometric function 2F1. 
The same scaling function was obtained independently by Toyoki We will call 
it the 'BPT scaling function'. The result ( |122| ) for 7 implies L{t) ~ t^/^ for all n 
within th is ap proxi matio n. 

Both ( |119| ) and ( [L29| ) suffer from the weakness that (for scalar fields) the width 
of the interface changes systematically with time. Since 0(m) is linear in m for 
small m, and |V0| is fixed (by the interface profile function) in the interface, we 
expect ((Vm)^) = const. Eqs. ( |119| ) and (|129|) , however, give ~ t-('^+2)/2 and 
~ exp(2t)/t'^'^"'"^^/^ for this quantity, corresponding to increasing and decreasing in- 
terface widths respectively. Oono and Puri |9^ showed t hat t his unphysical feature 
can be eliminated by introducing an extra term h{t) m in (|119|) . Since this term van- 
ishes at the interfaces, where m = 0, it's inclusion does not change the underlying 
physics. Fixing h{t) by the requirement ((Vm)^) = const, gives h{t) ~ ((i + 2)/4t at 
late times. The scaling function ( |126|) , however, is unaffected by the presence of the 
extra term. In section ^.2| we shall find that the Oono-Puri result arises naturally 
within a systematic treatment of the problem. 



5.1.3 Mazenko's Method 

In an interesting series of papers, Mazenko [p^,j88|, p9| has introduced a new approach 
that deals with the interface in a natural way. Tms approach combines a clever choice 
for the function 0(m) with the minimal assumption that the field m is gaussian. 
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Specifically 0(m) is chosen to be the equilibrium interface profile function, defined 
by (compare Eq. (|T0|)) 

0"(m) = V'(0) , (132) 

with boundary conditions 0(±cxd) = ±1, 0(0) = 0. The field m then has a phys- 
ical interpretation, near walls, as a coordinate normal to the wall. Note that this 
mapping transforms a problem with two length scales, the domain scale L{t) and 
the interface width ^, into one with only a single length scale, namely L{t) (see Fig. 
13). With the choice (|132|) for (f){m), the TDGL equation (H) becomes 

dt(j) = - 0"(m) . (133) 

Multiplying by at a different space point and averaging over initial conditions 
gives 

(l/2)9iC(12) = V2C(12) - (0"(m(l)) 0(m(2))) . (134) 

So far this is exact. In order to simplify the final term in (|134|) , Mazenko assumes 
that m can be treated as a gaussian field. Then the final term can be expressed in 
terms of C(12) itself as follows, exploiting the Fourier decomposition of 0(m) and 
the gaussian property of m 



(0"(m(l))0(m(2))) = Y.^k,Ki-kl){eM^kim{l) + ik2mi2)]) 



ki,k2 



2 9C(12)/95o(l) . (135) 



where 5*0(1), 5'o(2), are given by ( |12^ ) and Co(12) = (m(l)m(2)). The deri vativ e 
in ( |135| ) is taken holding 5*0(2) and 6o(12) fixed. Since, from the definition ( |122|) , 

7(12) = Co(12)/^5o(l)5o(2), the general result ( p.26| ) for gaussian fields implies 

dC{V2)_ ^ dC{l2) ^7(12) 
dSo{l) rf7(12) 95o(l) 

where 

a{t) = l/5o(l) = {m{lf)-' . (137) 
Putting it all together, and suppressing the arguments, the final equation for C is 

{1/2) dtC = V^C + a{t)-fdC/dj . (138) 

Using (|12|) for ^(7) gives -fdC/d-f = (2/7r) tan[(7r/2)(:7]. Then (|13|) becomes 
a closed non-linear equation fo r C . For a scaling solution, one requires L(t) ~ t^/^ 
and a{t) = X/2t for large t in ( |138| ), so that each of the terms scales as 1/t times a 
function of the scaling variable r/t^/^. Setting C(r,t) = /(r/t^/^) gives the equation 

for the scaling function f{x). The constant A is fixed by the requirement that the 
large-distance behaviour of C be physically reasonable [Q. Linearization of ( |139| ) 
(valid for large x) leads to two linearly independent large-x solutions with gaussian 
and power-law tails. The constant A is chosen to ehminate the 'unphysicaF power- 
law term. 
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It is straightforward to adapt this approach to nonconserved vector fields [5^, 
50| . A significant simplification is that for gaussian fields, the joint probability 
distributi on fo r m(l) and m(2) factors into a product of separate distributions of 

]138| ) for any 
:!26|). Again, 



ing function 



the form ( |123l ) for each component. This resu lts is an equation of form 
n, but with the function C{'y) given by ( |131| ) for general n instead of 
a(t) = X/2t, with A chosen to eliminate the power-law tail in the sea 
f {x). The values A for various n and d are given in table 1. 

It is interesting that the 'unphysical' power-law tails in real space become phys- 
ical when sufficiently long-range spatial correlations are present in the initial state. 
This will be shown using Renormalisation Group arguments |^ in section |8.3.1| . It 
also emerges within the Mazenko treatment 

The general two-time correlation function C7(r, ti, ^2) can al so be evaluated within 
It is given by a simple generalisation of ( |138|) , namely 



this scheme 59, 60 



dt^C = V^C + a{ti)-fdC/d-f 



(140) 



with a{ti) = X/2ti. This equation simplifies for ti ^ t2, because C is then small and 
the linear relation between C and 7 for small C (see Eq. ( |131|) ) implies 'ydC/d'j = C, 
i.e. 



dt,C = V'C + {X/2ti)C , ti>t2 



(141) 



This linear equation can be solved by spatial Fourier transform. Choosing an initial 
condition at ti = at2, with a ^ 1 to justify the use of ( |141| ) for all ti > at2, gives 



S(k,ti,t2 



at2 



A/2 



exp{-A;^(ti - at2)} S(k,at2,t2) 



(142) 



Imposing the scaling form ^(k, 0^2,^2) = ^2^^5'(^^^2), with g{0) = constant, and 
Fourier transforming back to real space gives, for ti ^ at2, 



C(12) = constant ( — 



{d-X)/2 



(143) 



The constant cannot be determined from the linear equation alone: it is, of course, 
independent of a. 

Comparison of ( |143| ) with the general form (B^), shows that X = d — X, i.e. the 
parameter A of the Mazenko theory is precisely tne exponent A associated with the 
response function G(k,t) (Eq. (p^)), related to A by (pSj). This connection was 
first pointed out by Liu and Mazenko [Q. The values of A obtained (table 1) are 



in reasonable agreement with those extracted from simulations |T9|, |70|, ^ p^]. For 
example, for the scalar theory in d = 2 simulations \^^, Wl\, p6| give A ~ 0.75 (argued 
to be 3/4 exactly in |]19|), compared to 0.711 from table 1. 



n=l 



n= 



0.301 0.378 0.414 



2 0.711 0.829 0.883 0.912 



3 1.327 1.382 1.413 1.432 



Table 1 Exponent A within the Mazenko theory. 
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To summarise, the virtues of Mazenko's approach are (i) only the assumption 
that the field m is gaussian is required, (ii) the scaling function has a non-trivial 
dependence o n d (w here as, ap art fr om the trivial dependence through the diffusion 
constant D, ( |122| ), ( [126| ) and ( |131D are independent of d), and (iii) the non-trivial 



behaviour of different-time co rrel ation functions [M emer ges i n a natural way M . 
In addition, the OJK result ( |126| ), and i ts g eneralisation (|131| ), are reproduced for 
d ^ oo, while the exact scaling function flTOB) of the 1—d Glauber model is recovered 



d oo, w hile the exact scaling function ( |10^ ) of the 1—d Glauber model is recovered 
from ( |139| ) in the limit d ^ 1 |99 1 . In practice, however, for d > 2 the shape o f the 



scaling function f(x) differs verjTlittle from that of the OJK function given by 
and ( |122| ), or its generalization ( |131D for vector fields p^. All these functions are 



in good agreement with numerical simulations. The Mazenko function for n = 1, 
d = 2 is included in Figure 12, while the BPT results for vector fields are compared 
with simulations in Figures 14 and 15. The Mazenko approach can also be used. 



with some modifications, for conserved scalar |^ and vector [jTOf fields. 



To conclude this section we note that the crucial gaussian approximation, used 
in all of these theories, has recently been critically discussed by Yeung et al. ||101| . 



By explicit simulation they find that the distribution P(m) for the field m at a single 
point is flatter than a gaussian at small m. In section y we shall show that the joint 
distribution P(m(l),m(2)) can be calculated analytically when |m(l)|, |m(2)| and 
r are all small com pare d to L{t). The result is non-gaussian, but is consistent with 
the gaussian form ( |123| ) in the limit ci — > oo. Below, we present evidence that the 
gaussian approximation becomes exact for ^ oo. Finally we note that very recent 
work by Mazenko presents a first attempt to go beyond the gaussian approximation 

Hi. 



5.2 A SYSTEMATIC APPROACH 

All of the treatments discussed above suffer from the disadvantage that they invoke 
an uncontrolled approximation at some stage. Very recently, however, a new ap- 
proach has been developed which recovers the OJK and BPT scaling functions 



in leading order, but has the advantage that it can, in principle, be systematically 
improved. 

5.2.1 Scalar Fields 

For simplicity of presentation, we will begin with scalar fields. The TDGL equation 
for a non-conserved scalar field 0(x,t) is given by Eq. (H). We recall that, accord- 
ing to the Allen-Cahn equation (|18|), the interface motion is determined solely by 
the local curvature. It follows that the detailed form of the potential V{(f)) is not 
important, a fact that we can usefully exploit: the principal role of the double- well 
potential is to establish and maintain well-defined interfaces. 

11 w- 
TWe 



Following Mazenko 
conditions (l)(±oo) = ± 



we define the function (j){m) by Eq. ( |13^ ) with boundary 
e have noted that 0(m) is just the equilibrium domain- 
wall profile function, with m playing the role of the distance from the wall. Therefore, 
the spatial variation of m near a domain wall is completely smooth (in fact, linear). 
The additional condition 0(0) = locates the center of the wall at m = 0. Figure 
13 illustrates the difference between and m for a cut through the system. Note 
that, while saturates in the interior of domains, m i s typ ically of order L{t), the 
domain scale. Rewriting (H) in terms of m, and using (|132|) to eliminate V', gives 



dm = V'm - (1 - (Vmf) . (144) 

0'(m) 



For general potentials V{(j)), Eq. ( | 144|) is a complicated non-linear equation, not 
obviously simpler than the original TDGL equation (0). For reasons discussed in 
section p.3| , however, we expect the scaling function f{x) to be independent both of 
the detailed form of the potential and of the particular choice for the distribution 
of initial conditions. Physically, the motion of the interfaces is determined by their 
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curvature. The potential V{(f)) determines the domain wall profile, which is irrelevant 
to the large-scale structure. 

Similarly, the initial conditions determine the early-time locations of the walls, 
which should again be irrelevant for late-stage scaling properties. For example, in 
Mazenko's approximate theory, both the potential and the initial conditions drop 
out from the equation for the scaling function f{x). 

The key step in the present approach is to exploit the notion that the scaling 
function should be independent of the potential (or, equiv alent ly, independent of the 
wall profile) by choosing a particular ^(0) such that Eq. ( |144| ) takes a much simpler 
form (Eq. ( |148D ) . Specifically we choose the domain-wall profile function 0(m) to 
satisfy 

0"(m) = -m0'(m) . (145) 



This is equivalent, via (|132| ), to a particular choice of potential, as discussed below. 
First we observe that ( |145| ) can be integrated, with boundary conditions 0(±oo) = 
±1 and 0(0) = to give tlie wall profile function 

rm 

0(m) = (2/7r)^/2 / exp(-xV2) = erf {m/^y2) , (146) 
Jo 



where erf (x) is the error function. Also, (|132|) can be integrated once, with the zero 



of potential defined by V^(±l) = 0, to give 

\/(0) = (1/2) (00' = (1/^) exp(-m2) = (I/tt) exp(-2[erf-i (0)]^) , (147) 

where erf^^ (x) is the inverse function of erf (x). In particular, V{(f)) ~ l/vr — 0^/2 
for 02 < 1, while V(4)) - (1/4)(1- ln(l - 0^)1 for (1 - 0^) < 1 M . 



With the choice (|145|) , Eq. (|146|) reduces to the much simpler equation 

dtm = V^m + {1- {Vmf)m . (148) 

This equation, though still non-linear, represents a significant simplification of the 
original TDGL equation. It is clear, however, on the basis of the physical arguments 
discussed above, that it retains all the ingredients necessary to describe the universal 
scaling properties. 

We now proceed to show that the usual OJK result is recovered by simp ly re- 
placing (Vm)2 by its average (over the ensemble of initial conditions) in (|148|) , and 
choosing a gaussian distribution for the initial conditions. In order to make this 
replacement in a controlled way, however, and to facilitate the eventual computa- 
tion of corrections to the leading order results, we systematize the treatment by 
attaching t o th e field m an internal 'colour' index a which runs from 1 to A^, and 
generalize ( p.48|) to 

N 

dtrria = V^m^ + (1 - iV"^ ^ (Vm^)^) m„ . (149) 

/3=1 

Eq. (|148|) is the case = 1. The OJK result is obtained, however, by taking the 
limit A^ — > oo, when A^~^ Z]^=i(Vm/3)2 may be replaced by its average. In this limit 
( |148|) becomes (where m now stands for one of the rria) 

dtm = V'^m + a{t)m (150) 
a{t) = 1 - ((Vm)2) , (151) 

a self-consistent /mear equation for m{'x,t). 

It is interesting that the replacement of (Vm)^ by its average in (|148|) is also 
justified in the limit d ^ oo, because (Vm)^ = YA=i{dm/ dxiY . If m is a gau s- 
sian random field (and the self-consistency of this assumption follows from ( |150| ) - 
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see below) then the different derivatives dm/dxi at a given point x are indepen- 
dent random variables, and the central limit theorem gives, for d — *• cx), (Vm)^ — *• 
d{{dm/dxiY) = ((Vm)^), with fluctuations of relative order l/y/d. While this ap- 
proach is not so simple to systematize as that adopted above, it seems clear that 
the leading order results become exact for large d. 

As discussed above, we will take the initial conditions for m to be gaussian, with 
mean zero and correlator (in Fourier space) 

(mk(0)m_k'(0)) = A5k,k' , (152) 

representing short-range spatial correlations at t = 0. Then m is a gaussian field at 
all times. The solution of (|150|) is 

mk(t) = mk(0)exp(-A;2t + 6(t)) , (153) 

b{t) = f dt'ait') . (154) 

Inserting this into (|151|) yields 







a{t) = db/dt = 1 - A ^ A;^ exp{-2kh + 2b) . (155) 

k 

After evaluating the sum one obtains, for large t (where the db/dt term can be 
neglected), exp (26) ~ {At/ Ad) (Svrt)'^/^, and hence a{t) ~ (ci + 2)/4t. This form for 
a{t) in ( |150| ), arising completely naturally in this scheme, reproduces exactly the 
Oono-Puri modification of the OJK theory [p2|, des igned to keep the wall- width 



5.1.2 



finite as t — oo, which was discussed in section 

The explicit result for mk(t), valid for large t, is 

mk(t) = mk(0) (4t/Arf)i/2 (^gj^tfi exp{-kH) , (156) 

from which the equal-time two-point correlation functions in Fourier and real space 
follow immediately: 

(mk(t)m_k(t)) = (4t/rf) (87rt)'^/2 exp(-2A;2t) , (157) 
(m(l)m(2)) = (4t/rf) exp(-rV8t) , (158) 

where '1', '2', are the usual shorthand for space-time points (ri,t), (r2,t), and 
r = |ri - r2|. 

We turn now to th e evaluation of the correlation function of the original fields 
Since, from ( 158|) , m is typically of order y/t a.t late times it follows from 



( |146|) that the field is saturat ed (i .e. = ±1) almost everywhere at late times. 
As a consequence, the relation ( |146| ) between and m may, as usual, be simpli- 
fied to = sgn (m) as far as the late-time scaling behavior is concerned. Thus 
C(12) = (sgn (m(l)) sgn (m(2))). The calculation of this average for a gaussian field 
m pr oceeds just as in the OJK calculation. The OJK result, given by (|126|) and 
( |122|) , (with D = 1) is recovered. The present approach, however, makes possible a 
systematic treatment in powers of 1/N. The work involved in calculating the next 
term is comparable to that required to obtain the 0{l/n) correction to the n = oo 
result for the 0{n) model pO|, 7T|. 



5.2.2 Vector Fields 

For vector fields, the TDGL equation is given by (|57|) , where V{(f)) is the usual 
'mexican hat' potential with ground-state manifold 0^ = 1. This time we introduce 
a vector field m(x, t), related to by the vector analog of (|132|) , namely ^| 



V^0 = dV/d(f) , (159) 
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where means YTa=\^'^ I ^"^a n-component field. We look for a radially 

symmetric solution of ( I 159| ), 0(m) = ihglp), with boundary conditions g{0) = 0, 
g{oo) = 1, where p = \m\ and m = m/ p. Then the function g{p) is the defect profile 
function for a topological defect in the n-component field, with p representing the 
distance from the defect core In terms of m, the TDGL equation for a 

vector field reads 

Just as in the scalar theory, we can attach an additional 'colour' index a (= 1, . . . , A^) 
to the vector field m, such that the theory in the limit — > oo is equivalent to 



replacing Vm-b ■ Vmc by its mean, ((Vmfc) ) (5b e in ( |16CI|) . Noting also that ((Vmb) ) 




is independent of h from global isotropy, (|160| ) simplifies in this limit to 

where mi is any component of ffi. Finally, this equation can be reduced to the linear 
form ( |150|) , with m replaced by m, through the choice V^0a = —J2b{9(pa/dmb) mi, 

or, more compactly, = — (m ■ Vm) 0, to determine the function 0(m). Substi- 
tuting the radially symmetric form = rhg{p) gives the equation 

-^^ = 0, (162) 
P 

a generalization of (|145|) , for the profile function g{p), with boundary conditions 
g{0) = 0, g{oo) = 1. The solution is linear in p for p 0, while (^(p) ~ 1 — (n— l)/2p^ 

for p —» oo . The potential V^(0) corresponding to this profile function can be deduced 
from ( p.59| ), though we have been unable to derive a closed form expression for it. 
Note that we are making here the natural assumption that scaling functions are 
independent of the details of the potential for vector fields, as well as for scalar 
fields. 

For the vector theory, Eqs. ( |15CI| ) and (|151|) hold separately f or ea ch component 
of the field. Taking gaussian initial conditions, with correlator ( |120|) , yields a{t) ~ 
{d + '^)/4:t again, giving ( |158|) for each component. The final step, the evaluation 

of the two-point function C(12) = (0(1) • 0(2)), proceeds exactly as in the KYG 

treatment of section 5.1.2| : since \m\ scales as ^/t, we can replace the function 0(m) 



by m at late times. Then C(12) = (m(l) ■m(2)) in the scaling regime. The required 
gaussian average over the fields m(l), m(2) yields the BPT scaling function ( |131j ). 
Again, it can be systematically improved by expanding in 1/A^. 

5.2.3 The Porod Tail 

It is easy to show [^, ^] that ( |131j ) contains the singular term of order 



(with an additional logarithm for even n) that generates the Porod tail ( ffOD in 
the structure factor. This feature was effectively built into the theory through the 

mapping 0(m). Specifically, the singular part of (|131D for 7 ^ 1 is [|104| 

/ 11. -t- X, I \ ^ 

B 



■"sing 



2tx 



Using 7 = exp(— r^/8t) = 1 — r^/8t + ■ ■ ■ for r <^ t^/^, and simplifying the beta and 
gamma functions, gives 



TT 



r(n/2) (4t)"/' 
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It will be interesting to compare this result with the exact short- distance singularity 
derived in section |[ 

In Figures 14 and 15, w e compare the BPT scaling function with numerical 
simulation results ||105| , |107|| , both for the pair correlation function C{r,t) and the 
structure factor S'(k, t). Since the defect density p scales natural choice 

for the scaling length L is Note that p can be measured independently in the 

simulation, so using rp^/" as scaling variable provides a direct, zero parameter test 
of the scaling hypothesis itself. For the scalar systems, the scaling variable r(l — cfP') 
was employed ||107|| : because 1 — cfp' is non-zero only near domain walls, (1 — 0^) is 
equal to p, up to a time- independent constant. 

The resulting scaling plots (Figure 14) provide very good evidence for scaling, 
except for d = 2 = n where clear scaling violations are apparent: the data drift 
to the right with increasing time, i.e. they are 'undercoUapsed'. In this case we 
can apparently make the data scale, however, by plotting against r/L(t) with L(t) 
chosen independently at each time t to provide the best data collapse. The collapse 
is then as good as for any of the other systems. 

The theoretical curves in Figure 14 are the BPT function ( |131|) , which reduces 
to the OJK scaling function for n = 1. In making the fits, 7 was replaced by 
exp(— Q;r^/L(t)^) with the scale factor a{n,d) adjusted to give the best fit by eye. 

The structure factor plots of Figure 15, on a log-log scale, confirm the existence 
of the Porod tail ( [70D in the data. On the logarithmic scale, the poor scaling of the 
d = 2 = n data against rp^/^ is reflected in a slight spreading of the data at small 
k in the corresponding structure factor plot. We do not show the structure factor 
plots for n = 1: the existence of the Porod tail for scalar systems is implicit in the 
linear short-distance regime in the real-space plots. 

It should be noted that the real-space correlation function (|131|) is independent 
of the space dimensionality d. The (i-dependence of the structure factor enters only 
through the process of Fourier transformation. Within the BPT theory, therefore, 
the Porod tail is obtained for any n and d. The same feature is present in the 
structure factor computed using Mazenko's method |5^, [6^]. In section ^, however, 
we shall show that the Porod tail is a direct consequence of the presence of stable 
topological defects in the system and, furthermore, that the amplitude of the tail 
can be evaluated exactly in terms of density of defect core. Since stable de fects are 
only possible for n < d, the Porod tail obtained from the BPT function ( p.31| ) for 
n > d is an artefact of the approximations invoked. This scenario is consistent with 
our claim that the BPT function actually represents an exact solution in the limit 
d — > 00. In this limit, of course, the condition n < li is always satisfied! 



5.2.4 External Fields, Thermal Noise, Quenched Disorder . . . 

Remarkably, the systematic approach can be readily extended to treat the situation 
where a general (space- and/or time-dependent) external field is present and/or the 
initial conditions contain a bias. This also allows the effects of thermal fluctuations 
to be incorporated to a limited extent. For simplicity we will only treat scalar fields. 
Consider the following equation of motion: 

d(j)/dt = V^(f) - V^{(j)) + /i(x, t)Vl{(f)) . (165) 

Here Vo(0) is the usual symmetric, double- well potential sketched in Figure 3, while 
Vi{4>) has the sigmoid form sketched in Figure 16(a). The full potential, V{(f); x, t) = 
Vo(0) — /i(x, )f:)Vi(0), has (for given x and t) the asymmetric double- well form shown 
in Figure 16(b), with the right-hand minimum lower for h > 0. 

As in our treatment of the case h = 0, we can exploit the insensitivity of the 
domain growth to specific details of the potential to choose especially convenient 
forms for Vq and Vi. This rests on the physical truth that the motion of an interface 
depends only on the local curvature K and the local field. To see this, consider again 
Eq. (p!?!) for the interface motion, this time for a general potential V{(j)). Multiplying 
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through by {d(l)/dg)t, integrating over g through the interface, and using (|TT|), gives 
the local velocity of the interface as 

v = -K + AV/a , (166) 

instead of (|18|), where AV is the potential difference across the interface. The 
essential point is that the interface motion depends on the external field only through 
AV. This gives us a great deal of flexibility in the choice of Vi{(f)), since all that 
matters is the potential difference between the minima of V{(j)). For example, we 
can choose the minima to remain at = ±1, as in Figure (16b). 

With these insights, we now change variables to the auxiliary field m, with 
(p = (j){m). Then ([165|) becomes 



0' dm/dt = ^'V^m + <i)"iymf - V^{<p) + /i(x, t)y/(0) , (167) 



where 0' = d(f)/dm etc. Simplifications analogous to those that led to (|148|) are 
achieved through the choices 

V^{(j)) = (f)" = -m(f)' (168) 

VM) = 0', (169) 

which give immediately 

= V^m+ (1 - (Vm)2)m + /i(x,t) , (170) 



a simple extension of (|148D . 



The right part of ( |168[ ) gives (with the appropriate boundary conditions) the 
usual error function profile (|146| ), while the left part leads to the previous form 
( p.471 ) for Vo(0). Integrating ( |i69D gives, with the boundary condition Vi(0) = 0, the 
result 



V^i(0) = - / dt exp(-t2) = — erf (m) = — erf [v^erf"^ (0)] . (171) 

TT Jo y/TT y/TT 

Again, this only defines Vi(0) for — 1 < < 1, but this is the only region we require 
for the T = dynamics. 

The difference Vi(l) — Vi{—1) is 2/y/7c, so the difference between the minima of 
the full potential, V = Vq — hVi, is —2h/y/7i, corresponding to an effective magnetic 
field heff = h/y/n as far as the interface dynamics are concerned. 
External Fields/Initial Bias 

As a simple application of ( p.7CI| ), consider the case /i(x, t) = h, representing a uni- 
form, time-independent magnetic field. In order to solve the equation, we take 
the same limit (d — *■ cxo, or number of 'colours', A^, large) as in section |5.2| , en- 
abling the replacement of (Vm)^ by its mean. Additionally, we allow for a bias, 
(m(0)) = mo(0), in the (gaussi an) initial conditions, while the other Fourier compo- 
nents (k 7^ 0) of m still satisfy (|152|) . Then the equations for the k 7^ components of 
m, and the self-consistency condition, are unchanged by the field: a{t) = 1 — (Vm)^, 
and b{t) = Jq dt' a{t') are the same as for = 0. The equation for the k = compo- 
nent is druo/dt + a(t)mo = h, with solution 

mo{t) = mo(0) exp{6(t)} + h dt' exp{6(t) - b{t')} . (172) 

JO 

Inserting the result exp(6) ^ {At / Ad)^^"^ (Srrt)'^^^ , valid for large t, from section 5^ 
gives, for large t, 

mo(t) = mo(0) (^— j iSntr^' + hj^ dt' (^-j , (173) 
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where to ~ (Afi)^/*^'^^^) is a short-time cut-off (to allow for the breakdown, at short 
times, of the form used for b(t)). 

Exploiting the gaussian property of m (which now has a non-zero mean given by 
( |173| )), we can calculate the expectation value of the original field 0: 

(0) = (sign (m)) = erf f ^tLj , (174) 

where {m'^)c = (m^) — (m)^ is the second cumulant of m. It is given by the same 
expression, Eq. (|158|) with 1 = 2, as for h = 0: {m?')c = At/d. So, 



(0) = erf 



The time-dependence of the mean order parameter (0) depends on d. Consider 
the argument of the error function. The initial bias mo(0) gives a contribution of 
order t'^/'* for any c?, but the contribution from the external field h scales as t^^^ for 
d < 2 (when times t' of order t dominate the integral in ( |175|) , as t^^'^\n{t/t(j) for 
d = 2, and as f^^^ for d > 2 (when the integral is dominated by times near the 
lower cut-off). Thus for t large the external field dominates over the initial bias for 
d < 2, whereas for d > 2 both terms are of the same order. For an arbitrary time- 
dependent field h{t), the final term in the argument of the error function is simply 
{d/8ty/^ Jl dt' h{t'){t/t'Y'^+^^/'^. This shows that the field becomes less important at 
late times and, for > 2, a constant field has all its effect at early times of order to. 
For d < 2, a. constant field continues to have an effect at late times. 



In fact one can make exact statements ||106|| about the 'initial growth' regime, 
where (0) <^ 1 . The main modification is that t'^^'^ (= L'^/^) gets replaced by 
L^, where A is the exponent in the scaling form (Q) for the response to the initial 
condition. This result is essentially obvious from tEe definition (^) of the response 
function. The crossover (in d) betwe en t he two regimes no lon ger occurs at d = 2, 



but at the dimension where A = 1 [p.06|| . The virtue of ( |175| ) is that it gives the 



complete time dependence, from the initial regime to final saturation ((0) = 1). 
Thermal Fluctuations 
Thermal fluctuations can be included, to some extent, within the present formalism 
by choosing /i(x, t) to be a gaussian white noise, with mean zero and correlator 
(/i(x, t)h{x.', t')) = 2D6{k — x.')6{t — t'). The original equation of motion ( |165| ) may 

be recast using ( [16^ ) and ([l69|) as d(j)/dt = VV - V^{(j)) + ^2Vo{^) /i(x,t). Recall 



that Vo(0) vanishes is the bulk phases, Vo(±l) = 0. The noise in the equation, 
therefore, also vanishes in the bulk phases, differing from zero only in the interfaces. 
Consequently, this noise will be effective in thermally roughening the interfaces, but 
will be incapable of nucleating bubbles of stable phase from a metastable state, or 
thermally exciting reversed regions within a domain. 
Quenched Disorder 

Quenched random fields are generated by a time independent, spatially random field 



/i(x). Again, in the original variable the field is multiplied by -y/2Vo(0), and so 

is active only at the interfaces. Since driving forces due to the field only act at the 
interfaces, this way of including a random field is perfectly adequate. Unfortunately, 
however, our leading order approximation of replacing (Vm)^ by its average misses 
the important interface pinning effects induced by the disorder, so this term has to 
be kept in full. A detailed discussi on of quenched disorder, using Renormahsation 
Group concepts, is given in section |8.3.2[ 
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5.3 HIGHER-ORDER CORRELATION FUNCTIONS 

Until now we have focussed exclusively on the pair correlation function C{r,t) and 
its the Fourier transform, the structure factor ^(k, t). These primarily probe the 
spatial correlations in the sign, or direction (for vector fields), of the order pa- 
rameter. However, one can also study the spatial correlations in the amplitude of 
the order parameter This is worthwhile for two reasons. In certain systems, 
such as superconductors and superfluids, the (complex scalar) order parameter ip 
does not directly couple to experimental probes. Rather, such probes couple to 
I -01^, and any scattering experiment, for example, measures the Fourier transform of 
(|'0(l)P|'?/'(2)p). The second reason to study these correlation functions is that the 
simultaneous calculation of two different correlation functions provides an exacting 
test of theory. This is because plotting one correlation function ag ainst another pro- 
vides an 'absolute' (i.e. free of adjustable parameters) prediction |107|| . Tested this 
way, the predictions of the gaussian theories of the 'OJK' and 'BPT' (or 'Mazenko') 
type are not quite as impressive as they at first seem. 

In this section we will be concerned specifically with the normalised correlation 
function 

_ ({i-0iin{i-y«> , 

(1-0(1)=) (1-0(2)=) 

where the '1' in each bracket represents the saturated (i.e. equilibrium) value of 

0^. The function C4(12) can be evaluated using any of the gaussian field methods 
discussed above ||55|. For definiteness, we adopt the 'systematic approach' of section 
|5l^ . The details of the calculation are qualitatively different for scalar and vector 
fields. 

5.3.1 Scalar Fields 

In terms of the gaussian auxiliary field m the numerator in ( |1 76| ) is given by 



= / dm{l) / dm{2) P(m(l), m(2)) {1 - 0(m(l))2}{l - 0(m(2))2} , (177) 



where P is the probability distribution ( 123 ). Since (1 — 0^(m)) approaches zero 



exponentially fast for scalar fields, the integrals are dominated by values of m(l) 
and m(2) close to zero (i.e. within an interfacial width of zero). The variation 
of P with m(l) and m(2), on the other hand, is set by the length scales r and 
L{t), which are both large in the scaling limit. Defining the interfacial width ^ by 
^ = / dm{l — 0(m)^) gives, in the scaling limit, 

C^ = eP{0,0)= , (178) 

/27r(l - 72) 



while the normalised correlator C4 is 

C4 = (1 - 72)-i/2 _ (179) 

Here we recall that 7 = 7(12) is the normalised correlator ( p.22| ) of the field m. In 
particular, 7(0) = 1 and 7(00) = 0. Using 7 = 1 — constr'^ /t for r <^ t^^'^, we see that 
Ci ^ L/r for r <^ L t^l"^ . This result will be derived using elementary arguments 



in section |6.4 Note that the 1/r dependence at small r implies a power- law tail, 
S'4(k) ~ Lk~^'''~^\ in the Fourier transform of C4. 

B y eli minating 7 between C4 and the pair correlation function C = {2/n) sin~^ 7 



[see ( |126|) ), we obtain the 'absolute' relation 

I/C4 = cos(7rC/2) (180) 
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betw een t he two correlation functions, with no adjustable parameters. We emphasize 
that ( |180| ) is a prediction of all gaussian theories, which differ only in the relation 

between 7 and the scaling variable x = r/t^/"^. Thus a test of (|180|) is a test of the 
gaussian assumption itself. 

In Figure 17 we show I/C4 plotted agains t C, where C4 and C wer e me asured 
simultaneously in 'cell dynamics' simulations ||108|| 'm d = 2 and = 3 ||105|| . Also 



shown is the prediction ( |180|) . It is clear that the agreement is much poorer than one 
obtains by fitting C alone (see Figure (12)). The agreement is significantly better, 
however, for c? = 3 than for d = 2, consistent with our claim in section ^]2| that the 
gaussian assumption becomes exact for c? — cx). 

5.3.2 Vector Fields 

The first step is a simple generalisation of ( |1 77| ) to vector fields: 

Cf = / dm{l) f dm{2) P(m(l), m(2)) {1 - ${m{l)f}{l - 4>{m{2)f} , (181) 



where P is a product of separate factors ( |1 23| ) for each component (since m is 
assumed to be gaussian). The subsequent analysis is different from the scalar case, 

however, because for vector fields ^(m)^ approaches its saturated value of unity for 

\m\ —>■ 00 only as a power law. To see this we recall that the function (f){m) is defined 
as the equilibrium profile function for a radially symmetric topological defect. The 

amplitude equation satisfies ( pSf ) with / — >• |0| and r —>■ |m|. From (^) we obtain 
directly 

1 -0(m)2 ^ ^Vl^r . |m| ^ 00 , (182) 

where = 2(?7. — l)/y(l). We will use this to define the 'core size' ^ for topological 

defects in vec tor f ields. 

Inserting ( |182| ) in ( |181| ), we see that for n > 2 the factors (1 — ^(m)^) do not, in 
contrast to scalar fields, converge the integral at small \m\ (i.e. at |m| ~ ^). Instead, 
the integrals are converged in this case by the probability distribution P, wh ich s ets 
a typical scale L(t) for |m|. This justifies the use of the asymptotic form ( |182| ) in 
the scaling limit: 



N 



(183) 

J |m(l)|2^ |m(2)|2 



It is now a straightforward matter to evaluate the m integrals Dividing by the 
large-distance limit (corresponding to 7 = 0), gives the normalised correlator ( |176D 

as 

C4 = P(l,l;n/2;72) , (184) 

where F is again the hypergeometric function 2P1. 

For 7 ^ 1, C4 has a short-distance singularity proportional to (1 — 72)("-4)/2 ^ 
(L/r)^^" (with logarithmic corrections for even n). It follows that the Fourier trans- 
form has the power-law tail S4 ~ j^i-nf^-id+n-i) n > 2. 

For the special case n = 2 one has to be more careful, as the integral ( |183| ) is 
formally logarithmically divergent at small |m(l)|, |m(2)|, and has to be cut off at 
^. A careful analysis |^ shows that exhibits logarithmic scaling violations 



\m 



in this case. However, in the scaling limit r 00, L(t) 00 with r/L{t) fixed, the 



extra logarithm cancels in the normalised correlator C4, and (|184|) is recovered, but 



with logarithmic corrections to scaling |105 



In fact ( |184D simplifies for physical (i.e. integer) values of n, giving (1 — 7^) 



for n = 2 and sin ^( 7)77 (1 — 7^)^^ for n = 3. As for the scalar theory, one can 
eliminate 7 between ( |184| ) and ( 131| ) to obtain a parameter- free relation between C4 
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and C that may be used as an absolute test of the gaussian assumption. Figure 18 



shows data for I/C4 plotted against 1 — C, from cell-dynamics simulations ||105 
and the corresponding predictions of the gaussian theory. It can be seen that the 
gaussian theory is again rather poor but, as for the scalar theory, it improves with 
increasing d, once more in accord with our argument that it becomes exact for large 
d. 

5.3.3 Defect-Defect Correlations 

As a final example, we consider the correlation functions of the defect density itself. 
In terms of the auxiliary field m, the defect density is p(x) = (5(m(x)) J, where 
J is the Jacobian between the field fh and the spatial coordinate x, for example 
J = |Vm| for scalar fields. A significant simplification is achieved by choosing 
Mazenko's definition of fh, near defects, as a coordinate normal to the defect. Then 
J = 1 holds identically at defects, giving simply p = 6{m). Making now the gaussian 
approximation for m, the one-point distribution function is 

P(m) = (27r^o)'"/^ exp(-mV25o) , (185) 

where Sq = [m?) is mean-square value of one component of m. Eq. (|158|) gives 
5*0 = 4t/(i. The mean defect density is, therefore, 

pTr' = {Km)) = P(0) = {d/S-ntr/^ . (186) 

The superscript 'gaussl' indicates that this is one way to calculate pdef within the 
gaus sian approximation. An alternative approach, pursued by Liu and Mazenko 
109|| , is to retain the Jacobian exphcitly. This gives a different result, because 



J = 1 at defects is true for the exact fh, but not for the gaussian approximation. 
With the Jacobian retained, the calculation has not been completed for general d 
and n. 

With J = 1, the pair correlation function is also trivially evaluated: 

(p(l)p(2)) = P(0, 0) = (d/Srrtr (1 - 7')""/' , (187) 
where we used ( |123| ) for P(0, 0). In the short-distance limit, r <^ t^^^, this becomes. 



-j , r«L(t). (188) 

Again, the short-distance behaviour can be evaluated exactly (see section [6^) , 
and the r~" short- distance behaviour recovered from simple geometrical arguments, 
which exclude, however, point defects (i.e. n = d). This failure to capture the 
correct short- distance behaviour for point defects is another weakness of the gaussian 
approximation. 

5.4 NEMATIC LIQUID CRYSTALS 

We have not succeeded in applying the systematic approach to the equation of 
motion (|74D for nematics. Application of the KYG method (see section |5.1.2| ), 



however, is relatively straightforward [B6||. For orientation purposes, we first recall 



the use of the KYG method for vector fields |5^. Recall that, in the scaling regime, 

the relation between the order parameter field (p and the auxiliary field fh can be 

simplified to = m, a unit vector, and that fh ma y be t aken to satisfy the diffusion 
equation dtfh = V^m. As was stressed in section 5.1.2|, this approach is somewhat 



ad hoc, and is not even guaranteed to yield the correct time-dependence for L{t). 
In practice, however, it gives good results for scaling functions since it builds in. 
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through the zeros of m, the correct topological defects. Therefore, we adopt this as 
a reasonable first attempt. It turns out that for nematics, we do in fact recover the 
correct growth, L ~ t^/^, as shown in section ^ 

The first step is to introduce the (traceless, symmetric) tensor auxiliary field m, 
satisfying the diffusion equation. The only tricky part is to determine the mapping 
Q{m), between the auxiliary field and the order parameter, analogous to (^{m) = 

sgn (m) for scalar fields and 0(m) = rfi for vector fields. The key observation is that 
these latter results simply represent the mapping from an initial value of m to the 
nearest minimum of the potential or, equivalently, they describe the attractors of 
the dynamics (|7^ for a spatially uniform initial state. It is easy to show that 



for a nematic, an equivalent procedure is the following. The director n at a given 
space-time point '1' is obtained as the eigenvector with largest eigenvalue of the 
tensor m(l) obtained by evolving the diffusion equation dtm = V^m forward from 
a random initial condition. The physical tensor Q{1) then has elements Qafe(l) = 
S [na{l)nb{l) — 5af)/3], where S is an arbitrary amplitude that has the value 3/2 
for the particular coefficients in the equation of motion (|7^. The pair correlation 
function is then obtained as 

C(12) = (2/3)(Trg(l)g(2)) , (189) 

where the factor 2/3 normalises (for = 3/2) the correlation function to unity when 
points '1' and '2' are the same. The average in (|189|) is over the (gaussian) joint 
probability distribution for m(l) and m(2), which can be deduced from the diffusion 
equation for m and the assumed gaussian initial conditions. 

The results for the pair correlation function and scaled structure factor are shown 
in Figure 19, along with the simulation data of Blundell and Bray ||67|, and the 
experimental structure factor data of Wong et al. [^]. The inset in Figure 19(a) 
shows that the real-space scaling function f{x) has the short- distance behaviour 
f{x) = 1 + ax'^lnx — hx"^ + ■ ■ ■. This is the same short-distance form as the 0(2) 
model and leads to the same tail in the structure factor, reflecting the presence 
of line defects (disclinations). The fit to the simulation data (with t he le ngth scale 
L{t) adjusted at each time) is good. Remarkably, the BPT function (|131|) for n = 2 



fits just as well, and indeed the simulation data for the two systems are essentially 
indistinguishable. This provides a dramatic illustration of the central role played 
by the topological defects: the nematic might naively be regarded as more like an 
n = 3 than an = 2 system. 

The theoretical curve in Figure 19(b) represents the 0(2) theory, as this was 
simpler to obtain, by numerical Fourier transform of the analytic result for the 
real-space scaling function, than the Fourier transform of the nematic correlation 
function, which had to be generated numerically ||6^ . Again the agreement is quite 



good. The data of Wong et al. can be shifted by an arbitrary amount, both hori- 
zontally and vertically, but we were unable to collapse it precisely on to the analytic 
result or simulation data. In addition, the experimental data have not yet reached 
the asymptotic regime expected on the basis of the string defects present. A 
line of slope -5 is included as a guide to the eye. 



5.5 CONSERVED FIELDS 

We have seen that the OJK scaling function ( |126| ) and its generalisation ( |131| ) to 
vector fields provide a very good description of the pair correlation function for 
nonconserved fields, subject to the caveat that the scale length L(t) is fitted when 
comparing with data. Furthermore, we have argued that the gaussian approximation 
for the auxiliary field fh is exact in certain limits, and provides a starting point for 
a systematic treatment. 

For conserved fields, the theory is less well developed. The most naive approach, 
for example, does not even give the correct growth law, L{t) ~ t^/^ (for scalar fields). 
One can still attempt to to make progress by introducing an auxiliary field m, but 
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in contrast to nonconserved fields, there is no evidence for any simple limit in which 
the theory for m becomes tractable. 

To put the difficulties into context, we start with scalar fields, and recall that the 
chemical potential /i satisfies the Laplace equation (|25|), with boundary conditions 
(p8|) imposed at the interfaces. In fact the interfaces act as sources of the field /z. To 
see this we integrate V^/i over a volume element dV enclosing an interface surface 
element dS. Using (^) gives a source density — 2f (r)(5(m(r))|Vm(r)|, where v is 
the interface velocity (measured in the direction of increasing (p) and m(r) is an 
auxiliary field whose zeros define the interfaces (so that 5(m(r))|Vm(r)| gives the 
volume density of interfacial area). Our usual choice for the field m, defined by 
(|132|), gives |Vm| = 1 at interfaces. This gives the Poisson equation 



V^fi = -2v{r)6{m{r)) , (190) 

with solution (for d = 3) 

1 r dr' 



Mr) = -/— ^^;(r')5(m(rO). (191) 
In J r — r 



The Gibbs- Thomson boundary condition gives /i = —aK/2 at an interface, where 
i^' = V ■ n = V^m i s the interface curvature. Using (|1 1 7|) (with |Vm| = 1) for the 
interface velocity in (|191|) , and the Gibbs-Thomson boundary condition for /i, gives 

(j„2 1 f dr' dm(r') ^ . . ... 

= i / I ^^-^ 

at interfaces. The same result could, in fact, be obtained directly from the Cahn- 
Hilliard equation, dt(f) = V^/i, by operating on both sides by the inverse Laplacian 
and setting /i = —aK/2 = — (cr/2)V^m on the interfaces. 

For nonconserved fields, the extension of the interface equation away from the 
interfaces (as in, for example, the OJK theory) is a mathematical convenience which 
does not change the underlying physics of interfaces moving under their local cur- 
vature. For conserved fields, however, the interface dynamics are nonlocal and the 
extension of ( |192| ) away from interfaces is non-trivial. Any extension should sat- 
isfy the following criteria: (i) the equation reduces to ( |192[ ) at interfaces, (ii) the 
chemical potential satifies V^/i = except at interfaces, (iii) V/i is discontinuous at 
interfaces, the normal component of the discontinuity generating the interface ve- 
locity as in (|29|) , (iv) the conservation of the order parameter should be preserved. 
Unfortunately, it is very difficult to construct an approximate theory that satisfies 
all these criteria, and I am not aware of any successful attempts. In addition, a 
good theory would ideally incorporate two further features: (v) the structure factor 
should vanish as k'^ for A; — > (see section ^.5.2|) , and (vi) the short-distance expan- 



sion of the real-space scaling function should cont ain (after the leading '1') only odd 
powers of r, the so-called 'Tomita sum rule' [|100|| , deriving from the smoothness of 
the interfaces. 

A number of approximate theories have been proposed which satisfy a su bset 
of these requi reme nts. The theories of Ohta and Nozaki ||110|| , Tomita ||111|| and 



Yeung et al. [|101|| all involve a gaussian approximation for the correlator of the 
auxiliary field m. The correct t^^'^ growth is obtained, and the scaling functions 
describe the real-space simulation data very well out to reasonable values of the 



scaling va riable, bu t violate the conservation law. A recent attem pt b y Kramer and 
Mazenko ||112| , |113|| corresponds to an off-interface extension of ([l92|) in which the 
left-hand side replaced by —fj, = ((j/2)(V^m -|- u(j)/L), with u a constant. The real- 
space scaling function is obtained by multipling through by at a different space 
point, and averaging both sides with the usual gaussian assumption for m. The 



45 



resulting scaling function has the desired fc^ small-fc form in Fourier space, although 
this be haviou r does not have the same origin as in the derivation of this form in 
section |5.5.2| . The real-space fit is not as good as earher theories. In addition, the 
chemical potenti al (ra ther than its gradient) is discontinuous at interfaces. 

Yeung et al. ||101|| have critically analysed approximation schemes based on a 
gaussian assumption for m. Using data of Shinozaki and Oono P3| for the pair 
correlation function C{r,t) in d = 3 to infer a value for the normalised correlator 
of m, namely 7 = sin(7rC/2) (see ( |126| )), they found that the Fourier transform 
7fc would have to be negative at small k to fit the data. However, this is impos- 
sible since 7fc > by definition. Yeung et al. concluded that no gaussian theory 
could adequately describe the data, at least for scalar fields in ci = 3. Nevertheless, 
we conclude this subsection by considering the approach of Mazenko for conserved 
fields, which is explicitly built on a gaussian assumption for m. This is especially 
interesting for vect or fi elds, since it allows us to make contact with the large-n cal- 
culation of section [4.3| . The reader may recall that the exact solution of Coniglio 
and Zannetti for conserved fields with n = 00 exhibits a novel 'multiscaling' 
behaviour. A natural question is whether this behaviour survives at finite n. Em- 
ploying the Mazenko approach, we find multiscaling behaviour for n strictly infinite, 
but conventional scaling for any finite n W^ . 



5.5.1 The Mazenko Method for Conserved Vector Fields 

A naive application of Mazenko's technique to the Cahn-Hilliard equ ation (|]), and 
its generalization to vector fields, yields, in complete analogy to (|138|) , 

(1/2) dC/dt = -V\V^C + a{t) -fdC/d-f) , (193) 

with a{t) still defined by ( |137| ). The only difference between ( |138| ) and ( |193D is the 
extra (— V^) on the righ t-han d side. The form ( |193| ) is obtained for any n, with the 
function C(7) given by ( |131|) . 



If we seek a scaling solution of ( |193| ), of the form C{r,t) = f{r/L(t)), it is 



immediately clear that consistency requires a{t) ~ 1/L(t)^ and L{t) ~ t^/^. We 
shall show in section |^ that this in fact is the correct growth law for n > 2, while 
there is a logarithmic correction, L(t) ~ (tlnt)^/^, for n = 2 (and d > 2). For 
n = 1, however, ( |193|) fails to give the correct t^^^ growth. The reason is clear. 



Taking to be a sigmoid function of a gaussian field from the outset overlooks 
the vital role of the bulk diffusion field in transferring material between interfaces. 
Recognizing this fact, Mazenko |^ writes as a superposition of 'ordering' and 

'diffusing' components, = ip{m) + 0, with m a gaussian field. It is then possible 

to construct a consistent th eory with t^^^ growth although the results do not 



agree well with simulations ||114|| 

Here we will concentrate on vector fields with n > 2, for which ( |193|) does give 

the correct t^^'^ growth. F or ge neral n this equation can be solved numerically for 
the scaling function f{x) [|115|| . Somewhat surprisingly, the solution for d = 3, 



n = 2 is very close to the simulation data of Siegert and Rao fTg], despite the 
(logarithmically) wrong growth law. For large n, however, we can make analytic 
progress, and co ntrib ute to the debate on the possibility o f mu ltiscaling for finite n. 

For large n, (|193|) is simplified as follows. Expanding ( p.31| ) to first order in 1/n 
gives C = 7-7(l-72)/2n + 0(n-2), and so -fdC/d'y = C + /n + 0{n-'^). Putting 
this in ( |193| ) gives 

(1/2) dC/dt = -V^C - a{t) V\C + C^/n) (194) 

correct to 0(l/?7,). The solution of this equation is very different for n strictly infinite 
than for n large but finite. 

For n = 00, the C'^/n term can be dropped and the resulting linear equation 
solved by Fourier transformation to give 

S{k, t) = S{k, 0) exp[-2kH + 2k%{t)] (195) 
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for the structure factor, where b(t) = j^dt' a{t'). This result is identical to the 
exact result obtained by Coniglio and Zannetti in the sam e lim it (compare for 
example, equation (|97|)), and leads to same multiscaling form ( |103| ). 

For n finite, we try scaling forms consistent with the expected t^^^ growth, namely 
^(k,^) = f^/^g^kt^/*), a{t) = qi/t^/\ and C{r,t) = f{r/t^^^), where f{x) is the 
Fourier transform of g{q) and is a constant. Using these in ( 194| ) gives 



Big) = i8qi/n)if\, (197) 



dg/dq = -id/q + 8q'-8qiq)g + qB{q) , (196) 

,2 /„\ I eZ\ 



where (/^)q indicates the Fourier transform of f{x)^. Note that g{0) = must hold 
for a conserved order parameter, otherwise ^(k, 0) would grow as L{tY, violating 
the conservation law. Integrating (|196| ) with initial condition g{0) = gives 

g{q) = q-" exp(-2g^ + Aqlq^) P dq' q"'+'B{q') exp(2g'^ - Aqlq'') . (198) 

JO 

The constant q^ is fixed by the condition /(O) = 1, i.e. I]qfl'(Q') = 1- For very 
large n we find a posteriori that is large. Then g{q) is strongly peaked near 
9 = Imi justifyi ng a steepest descent evaluation of the sum over q. For q near g^, 
the integral in ( |198D is dominated by q' values of order g~^, giving 



g{q) ^ 2-('^+3)r(l + d/2)B{Q)q;<'+'h~' exp(-2g^ + Aqlq') (199) 

for q near g^- Using this form in I]q5'(g) = 1, and evaluating the sum by steepest 
descents, gives 



1 = 2-('^+^) v^irrfr(l + d/2)qj''+''^ exp(2g^)5(0) , (200) 

where Ka = 2/{A7rY/'^T{d/2). Using this to eliminate B{0) from (|T^) gives the 
desired scaling solution, valid for g~^ ^ q ^ q-m, 

g{q) = {A/K,V2^) ql q-" exp[-2(g2 - q^f] . (201) 

In the limit qm —* oo, the width of the peak at g = g^ vanishes as g~^, so in this 
limit we can write 

giq) g^--^ 5(g - g^) , g^ ^ 00 . (202) 



The final step is to use ( |197| ) with g = to obtain a second relation (in addition to 



poq ) between B{0) and g^. From (m, B{0) = {8q^Jn) jd'^xf{xf. Using ( p02|) 



for g{q) gives f{x) = A{qmx), where A{y) = const. jy{y)/y'', with v = [d — 2)/2, 
and hence 

5(0) = (8g^/n) j d'^xA{q^xf = const, q^'^/n . (203) 

Putting this in ( pOQ ) gives 1 = const, q^^^^'^^^ 6xp(2g^), giving g^ — [(lnn)/2]^/'^ for 
large n. This in turn implies a characteristic length scale L(t) ~ (t/lnn)^/*^. 

To summarise, we have shown that, within the Mazenko approximation, scaling 
solutions are obtained for any finite n. Only for n strictly infinite is the multiscaling 
form (|103|) of Coniglio and Zannetti recovered. Note that the amplitude of the t^/^ 
growth depends in a singular way on n (as [Inn]"^/^) for n — 00, i.e. our result is 
non-perturbative in 1/n and could not be obtained by expanding around the large-n 
solution. The scaling solution is obtained when the limit t ^ 00 is taken at fixed n. 



47 



5.5.2 The Small- A; Behaviour of the Structure Factor 

One slightly unsatisfactory aspect of the above treatment is that is does not recover 
the correct small-g behaviour (indeed, the same shortcoming afflicts the scalar ver- 
sion of the calculation For g ^ 0, Eq. ([198D gives g{q) B{0)q^/{d+2). There 



are c omp elling arguments, however, for a behaviour at small q [|116| , [T7| , [Tg, 117 



omp 

1181, 11911 , strongly supported by numerical simulations , as well as experiment 
1201 1 , for a scalar order parameter. Here we discuss both scalar and vector fields. We 



begin by deriving an i nequality for the small-g beaviour, using an approach based 



on that of Yeung ||116 



We recall that the equation of motion for conserved fields takes the form = 
V^/i, where fi is the chemical potential. Multiplying through by at a different 
space point, averaging, and Fourier transforming, gives 



1 dSik) 



2 dt 

< e[S{k)]'/'{f,^fi.^y/\ (204) 

where the final line follows from the Cauchy-Schwartz inequality. Now impose scal- 
ing: ^(k) = L'^ g{kL). Since /i ~ 1/L for scalar fields, the analogous scaling form 
is (/ik/i-k) = L'^^'^g^(kL). Putting these into (|204|) and using L ~ t^^^ for scalar 
fields, gives 

dg{q) + qg'{q) < const, [g{q) g^{q)Y''' . (205) 

For g — i> one expects g{q) ~ q^ and g^{q) — >■ const.., because /i is not a conserved 
field. Then (|205|) gives the inequality 5 > 4. An approximate treatment which gives 
the expected q^ small-g behaviour has been proposed by Kramer and Mazenko ||112 | 



For vector fields, it is shown in sections ^ and ^ that L{t) ~ t^/^ provided n > 2 
(see section ^ for a discussion of n = 2). This suggests that the (vector) chemical 
potential scales as /2 ~ 1/L^ for n > 2, which gives (|205|) again and 5 > 4 as before. 



It is easy to show, using an argument of Furukawa |T7|, T8| , |1 1 7|| , that the lower- 
bound for 5 implied by Yeung's argument is in fact realized, i.e. 5 = 4. Integrating 
the equation of motion 9f0k = —k^l^k gives 0k(^) = 0k(O) — /o /^k(^i)- This 
yields 



^(k, t) = -S(k, 0) + 2(0k(t)0-k(O)) + k^ f dt^ I dt2 (/ik(ti)/i-k(t2)) . (206) 

Jo Jo 

It is clear that the first two terms on the right must be neglible in the scaling regime, 
otherwise Yeung's inequality would be violated! It is simple, however, to show this 
explicitly. In the scaling limit (A; — *• 0, L — >■ oo with kL fixed), S'(k, t) increases 
as L'^. Therefore the time-independent first term on the right of (|206D is certainly 
negligible. Now consider the second term. The autocorrelation function, A{t) = 

((/>(x, t)0(x, 0)), decreases as A{t) ~ L{t)'^^. This implies that (0k(^)0-k(O)) = 
L'^~^a{kL), where a{q) is a scaling function. Since this term grows less rapidly 
than L"', it also is a negligible contribution to S'(k, t) for large L. It follows that 
the structure-factor scaling function g{q) is obtained entirely from the final term in 
( |^06| ). It vanishes as q^ because /i is not a conserved field, so (/Xk/^-k) is non-zero at 
A; = 0, as discussed above. Furukawa has gone slightly further, and shown explicitly 
that this final term in indeed of order L'^. Inserting the two-time scaling form (for 
scalar fields) (/Xk(^i)yU-k(^2)) = Lf~'^g^{kLi, L2/ Li) (a natural generalization of the 
equal-time scaling form given above), using L ~ t^/^, and evaluating the double time 
integral by power counting, gives S'(k, t) ~ /;;4^2+(a!-2)/3 ^ j^ij^d+i ];gqui];g(i_ 

Finally we note that Tomita has given a rather general argument for the A;^ 



behaviour based on the isotropy of the scaling functions ||118 |. 
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5.6 BINARY LIQUIDS 

The equation of motion appropriate to binary hquids, Eq. (pll), was derived in 



section |2.7| for the case where the inertial terms in the Navier-Stokes equation can 
be neglected. Eq. (^) leads to an asymptotic linear growth, L{t) ~ t. Here we will 
discuss how one might attempt to calculate an approximate scaling function for pair 
correlations in this regime. 

In the regime where L{t) ~ t, the 'advective term' in ( ^l]) dominates the 'diffusive 
term', AV^/i, on the right hand side, so we will discard the latter. In the spirit of 
the 'syst ematic approach' of section |5.2|, we introduce an auxiliary field m defined 
by ( |132|) , with the additional choice ( |145|) , corresponding to a convenient choice of 
the potential ^(0). Then the chemical potential /x can be expressed as 

H = 5F/5(t) = y'(0)-VV 

= -0'(m){V^m + [1 - (Vm)2]m} . (207) 

Expressing the equation of motion (^) in term of m, and using (|207| ) for fi gives 



dm{r) /dt = J dr' [Vm(r) ■ T(r - r') ■ V'm(r')] {0'(m(r'))}^ {V'^m(r') 

+ [l-(Vm(r'))V(r')} (208) 

Now we recall that {0'(m)}^ acts very much like a delta function on the interfaces. 
In fact, the result (compare Eq. (|TT|)) a = J dm{d(j) / dmY for the surface tension 
leads to the identification {d(j)/dmY = (y5{m). Using this in ( ^08| ) gives 

dm{r) /dt = a J dr [Vm(r) ■ T(r - r') • V'm(r')] 6{m{r')) V'^m(r') . (209) 

With the identification 0(r) = sgnm(r), this equation can be rewritten as an 
equation for 0, by multiplying both sides by (5(m(r)): 



dt 



^ J dr' [V0(r) ■ r(r - r') ■ V^r')] \/'^m{r') . (210) 



This equation could serve as a convenient starting point for approximate treat- 
ments of the pair correlation function C{r,t). Note, however, that Eq. ( piQ ) is 
fundamentally nonlinear, and approximate scaling functions that capture the cor- 
rect physics are difficult to construct. The correct growth law is, nevertheless, built 
in to ( ^1U| ): simple power-counting (remembering that m scales as L{t)) gives im- 
mediately L(t) ~ (Jt/rj, as required. 

In the above discussion, the field m was taken to have zero mean, corresponding 
to a critical quench. An off-critical quench can be handled in similar fashion by 



allowing m to have a non-zero mean ||121|| . An outstanding challenge is to devise 



an approximate treatment which includes the 'switching off' of the linear growth 
at small volume fractions, when the minority phase is no longer continuous, and 
to properly incorporate thermal fluctuations in this regime (see the discussion in 
section pTTj) . 



6 SHORT-DISTANCE SINGULARITIES AND POROD 
TAILS 

In the previous section various approximate treatments of correlation functions were 
discussed. In this section we show that exact statements can be made about the 
short distance behaviour or, more precisely, the short- distan ce s ingularities, of these 
functions. In particular, the qualitative arguments of section can be made precise 
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56| , and the amplitude of the k ('^+") Porod tail obtained in terms of the density of 
defect core, pdef, which scales as L~". The basic result is a generalization of ([68|), in 



which the leading singular contribution to C(r,t) is a term of order Irl*^ for n odd 
(or n real, in a continuation of the theory to real n), and jrl'^ln |r| for n even. This 

in turn implies a power law tail in S(k,t). 

We first illustrate the method for the case of point defects {n = d). The exten- 
sion to the general case n < d is relatively straightforward. Next we discuss the 
short- distance singularities of son ie hi gher-order correlation functions, namely the 
function C4 defined by equation (|176|) , and the defect-defect correlation function. 
We also calculate the 'Porod tails' in the corresponding structure factors. Finally, 
we compute for scalar fields the joint probability distribution P(m(l), m(2)) of the 
auxiliary field m of section |^, in the limit where |m(l)|, |?t2,(2)| and the distance 
r = |ri — are all small compared to L{t). We find that the distribution is not 
gaussian, except perhaps for d = 00. 



6.1 POINT DEFECTS (n = d) 

Consider the field at points x and x + r in the presence of a point defect at the 
origin. We consider the case where |x|, |x + r| and |r| are all small compared to a 
typical inter-defect distance L, but large compared to the defect core size ^. Then 
the field at the points x and x -|- r is saturated in length (i.e. of unit length) and 
not significantly distorted by the presence of other defects. Moreover, the field can 
be taken, up to a global rotation, to be directed radially outward from the origin, 
as illustrated in Figure 10. Thus 

0(x).0(x + r) = i^f^^. (211) 
|x| |x + r| 

With r held fixed we average ( |211|) over all possible relative positions of the point 
defect, i.e. over all values of x within a volume of order around the pair of points, 
with the appropriate probability density pdef- Focussing on the singular part of the 
correlation function we obtain 

Csing{r,t) = Pdef J c/'^x ^^j^^^^ - analytic tcrms j . (212) 

The 'analytic terms' in ( pl2|) serve to converge the x-integral at large |x|, and allow 
us to ex tend the integral over all space. We include as many terms in the expansion 
of ( |211|) ) in powers of r as are necessary to ensure the convergence of the integral. 



When n is even, there is a residual logarithmic singularity. This case can be retrieved 
from the general n result by taking a suitable limit (see below). 

At this point two comments are in order. Firstly, by taking the field to be directed 
radially outward from the origin, we seem to be limiting ourselves to 'defects', and 
excluding 'antidefects'. The antidefect of a point defect, however, can be generated 
(up to arbitrary rotations), by 'inverting' (0j —(pi) an odd number of Cartesian 

components of the vector 0. Reference to Figure 10(a) shows immediately that, 
for = 2 = d for example, the 'antivortex' can be generated from the 'vortex' by 

this construction. Clearly, however, the scalar product 0(1) • 0(2) required for the 
evaluation of Cging is invariant under this operation. Secondly, we seemed to need 
the assumption that the field near a given defect is not significantly distorted by 
the presence of other defects. Actually, this is not strictly necessary. Any distortion 
generated by the other defects provides an analytic background fiel d tha t does not 
affect the contribution of the given defect to the singular part of C ||122| . 



The integral ( pl2| ) (extended over all space) is evaluated in The result is 



C.„, = n."/-^5(^,^) r(-^) pdef\r\\ (213) 
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where r(x) is the gamma function, and B{x,y) = T{x)T{y)/T{x + y) is the beta 
function. Note that the dependence on |r| can be extracted simply by a change of 
variable in (|2T2|) . 

The pole in the r(— n/2) factor for even values of n signals a contribution of the 
form |r|" ln(|r|/L) to Csing for those cases. We shall discuss these cases explicitly 
when we have the result for general n < d. 



6.2 THE GENERAL CASE n<d 

For n < d the defects are spatially extended, but the analysis is only slightly more 
complicated. The defect defines a surface, or subspace, of dimension d — n in the 
(i-dimensional space. On the length scales of interest (small compared to L{t)), the 

defect is effectively 'flat' (walls) or 'straight' (lines), etc., and the vector (p can be 
taken to lie in the n-dimensional subspace orthogonal to the defect (the 'orthogonal 
subspace'). The vector r can be resolved into components r_|_ and r|| lying in this 
subspace and in the d — n dimensional subspace of the defect respectively. Now 
consider the points x and x + r, where x lies in the orthogonal subspace, with the 
origin of x lyi ng on the defect (see the illustration in Figure 20 for the case d = 3, 



n = 2). Then ( 211 ) has the same form, but with r replaced by r_|_ on the right-hand 
side. Proceeding as for point defects, the integration over the n-dimensional vector 
X, with r fixed, gives 

C.., = nW-^5(^,^) r(-^) p.e/|rxr , (214) 

where pdef is, as usual, the density of defect core. 

The final step is to take the isotropic average of ( pi4| ) over the orientations of 
r, i.e. to compute (|r_|_|'^) where the brackets indicate an isotropic average. For 
generality, and because we will need it later, we in fact compute (|r_|_|") for general 
a. To do this we set up generalized polar coordinates with the first d — n polar axes 
in the subspace of the defect. Then 

d—n 



r±\ = rY[\ sin6'i| , 

Jo^^ dOi (sin^i)°+'^-i-* 

.=\~f^^^de~{Mif^ 
T{d/2)T{{a + n)/2) 



i=l 
d—n /-t/S 

{\r±n = r-n 



r(n/2)r((a + c/)/2) 



(215) 



Using ( |215| ) with a = n to perform the isotropic average of ( |214| ) gives the final 
result for the singular part of the correlation function, valid for all n < d: 

^ T{{d + n)/2)T{n/2) ^ ' ' ^ ' 



which reduces to ( |213|) for n = d. 



We remarked in the previous section that for even n the leading singularity is of 
the form r'^lnr. The precise result can be extracted by setting n = 2m + e, with m 
an integer, letting e go to zero, and picking up the term of order unity. The leading 
pole contribution, proportional to (r^)™, is analytic in |r| and therefore does not 
contribute to Csing- The 0(1) term (in the expansion in powers of e) generates the 
logarithmic correction from the expansion of |rp™'+'^. This gives, for even n, 

Cs^ng = - (4/n) tt"/^'^ (-1)"/^ ^{^{^^^^/f ^ p,, .r" lu r . (217) 
stng \ I ) \ ) Y{{d + n)/2)T^{n/2)'^ * ^ ' 
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It turns out that the Fourier transform S'(k, t) of Csinn{T'^-,t) has the same form for 
even, odd and real n, so we will not need to cons ider the even n case separately^ 

We can now compare the exact result ( pi6| ) with the equivalent result ( |164| ) 
obtained within the gaussian theory for nonconserved fields. Equating these provides 
another way of estimating the defect density within the gaussian theory, namely 

gaussZ _ 

^^^^ " (47rt)"/2 r{d/2) • ^ 

Comparing this with ([186D , we see that the two estimates differ for general d, but 
agree for oo. This is another indication that the gaussian approximation 

becomes exact in this limit. 



6.3 THE STRUCTURE FACTOR TAIL 

It remains to Fourier transform ( pi6|) to obtain the tail of the structure factor. 
Although the Fourier transform of ( P16|) by itself does not technically exist (be- 



cause the required integral does not converge), the following method gives the large- 
momentum tail correctly, as may be checked by back-transforming the result. Of 
course, the Fourier transform of the complete correlation function C(r, t) does exist, 
since C{r,t) vanishes at infinity. 



Simple power counting on ( pi6| ) gives immediately the power-law tail ^(k, t) ~ 
j^-{d+n) ^ To derive the amplitude we employ the integral representation 

/■oo 

r(-n/2)|r|" = / duu''^'^'^ {exp(-Mr^) - analytic terms} , (219) 

JO 

where 'analytic terms' indicates, once more, as many terms in the expansion of 
exp(— -ur^) as are necessary to converge the integral. These terms will not contribute 
to the tail of the Fourier transfor m, an d can be dropped once the transform has been 
taken. The Fourier transform of (El 91) is, therefore. 



j duu ^ J d'^r exp{-ur'^ - ik ■ r) 

POO 

7r"/2 / dnn-('^+")/2-%xp(-A;V4n) 





Inserting the remaining factors from (|216|) gives the final result. 

We note that this expression is smooth as n passes through the even integers. The 
generality of the result should be noted: it is independent of any details of the 
dynamics, e.g. whether the order parameter is conserved or non-conserved, and holds 
independently of whether the scaling hypothesis is valid. We note that, as well as 
providing an exact result against which to test approximate theories, Eq. ( ^21| ) can 
also be used to determine the defect density experimentally from scattering data. 

Measure men ts of the tail amplitude in num erica l simulations are in good agree- 
ment with ( p21| ) for both scalar ||5^ and vector [|105|] fields. As an example, we show 



in F igure 21 the results for the cases n = d = 2, n = d = 3 and n = 2, d = 3 
105 1 . It is interesting that in all cases the asymptotic behaviour is approa ched from 
above. For a scalar order parameter (n = 1), the leading correction to (|221|) has 

been calculated by Tomita ||10CI|| . It is a term of order k~^'^^^\ associated with the 



curvature of the interfaces and obtained from a short- distance expansion of the form 
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C(r = 1 — ar + br^ — ■ ■ ■. The absence of an term leads to the 'Tomita sum rule' 
for the structure factor ||10CI|| , dk [k'^^^S{k) — A] = 0, where A is the amplitude 
of the Porod tail. 

The discussion has so far been restricted to the cases n < d, where singular 
topological defects exist. What can be said about the structure factor tail for n > dl 
The case n = d + 1 may be complicated by the presence of topological textures 
|123|| . For n > d + 1, preliminary numerical studies for nonconserved dynamics ||124 



suggest that the structure-factor scaling function g{y) has a 'stretched-exponentiai 
tail, g{y) ~ exp{—y^), with 5 ~ 1 for d = l,n = 3 and 5 ~ 1/2 for d = 2,n = 4. 
Similar studies by Toyoki were analysed in terms of a power-law tail, g{y) ~ w 



with X > d + n, but are also consistent with stretched exponential decay [|125 



In contrast to n < d, the tail behaviour for n > d may be different for conserved 



and nonconserved dynamics. Recent results for conserved dynamics fT^j suggest a 
stretched exponential form but with 6 ~ 1.7 for both d = 2,n = 4 and d = l,n = 3, 
while 6 ~ 2.7 for d = l,n = 2. 

6.4 HIGHER-ORDER CORRELATION FUNCTIONS 



We consider first the correlation function C4, defined by ( 176 ), of the square of the 



order-parameter. We concentrate here on the numerator , defined by, 

= ({1 - - 0(2)^) • (222) 

For n < d, the presence of topological defects leads to a singular short- distance 
behaviour that can be evaluated in direct analogy to that of the usual pair correlation 
function C. As in t he a pproximate calculation of C4, using gaussian auxiliary field 
methods, in section |5.3|, we have to distinguish between scalar and vector fields. 



6.4.1 Scalar Fields 

For scalar fields, 1 — 0^ is sharply peaked near domain walls. It is convenient to 
introduce the auxiliary field m defined as in Mazenko's approximate theory, but not 
assumed to be gaussian! We recall that the function 0(m) represent the equilibrium 
domain wall profile, with m the coordinate normal to the wall. Since 0^ saturates to 
unity with a width of order ^ of the wall, we can use, for the calculation of scaling 
functions, 

1 - (p'^(m) = ^5{m) , (223) 

where we have used our usual definition of the 'wall width', ^ = dm (1 — 0^(m)). 
Putting (HI) in (|222D gives 



Cf(l,2) = e{5{mil))5im{2))) (224) 
= ^'^^(1,2) (225) 
= ^2p(i|2)p(2) . (226) 

Here we have used the fact that the wall area density (area per unit volume), 
5(m)|Vm|, can be simplified to 6{m) since |Vm| = 1 at interfaces from the def- 
inition of m as a coordinate normal to the interface. In ( p26| ), therefore, P(l,2) 
indicates the joint probability density to find both '1' and '2' in an interface. In the 
final equality, P(l|2) is the probability density to find '1' in an interface given that '2' 
is in an interface. Clearly P{2) = p, the average wall density. For r -C L{t), P(l|2) 
is dominated by cases where '1' and '2' lie in the same wall, which can be regarded 
as fiat on this scale, as illustrated in Figure 22. For general d, P(l|2) = Sd-i/Sar, 

where Sd = 2'k'^''^ /T{d/2) is the surface area of a d-dimensional sphere of unit radius. 
Assembling everything in ( p26| ). 



C^^ep^. i^r<^L{t). (227) 
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This result breaks down when r becomes comparable with ^, since it is n o longer 
adequate to neglect the thickness of the walls. The small-r behaviour ( |227| ) implies 
the power-law tail Si(k,t) ~ k~^'^~^^ for the corresponding structure factor, valid 
when fc^ -C 1 -C kL{t). 

To obtain the normalised correlation function C4, we divide (|227|) by its large-r 
limit (1-02)2. 

(1 - = ^(5(m)) = , (228) 

giving 

C,^^, e « r « Lit) , (229) 

from which the wall width ^ has dropped out. The exact result for the tail of the 
corresponding structure factor is obtained by Fourier transformation. Inserting the 
expression for Sd gives 

^4(k, t) (2/p) (47r)('^-2)/2 p(^/2) ^-('^-i) , fee < K kL{t) . (230) 

This could be measured by small-angle scattering in a situation where all the scat- 
tering was from the interfaces. The d = 3 result has be en de rived by Onuki, together 
with the leading correction, of relative order l/(/cL)2 |p.26|| . 

A heuristic derivation of the k~^^~^^ tail, based purely on scaling, proceeds as 
follows ra. Since (1 - 0^) ^ p ^ i/^, Cf has the scaling form Cf = L-^f{r/L), 
giving 5*4 = L'^~'^g(kL). But for kL ^ 1, the scattering intensity should scale as the 
defect density, i.e. as 1/L. This requires g{y) ~ y-('^-i) for ?/ 1. This argument 
can be generalised to vec tor f ields, as we shall see. 

It is interesting that ( |229| ) can be combined with ( pi6| ) (with ri = 1) to obtain 
an exact relation, valid at short distances, between the two correlation functions. 
For n = 1, ( |216| ) implies the short- distance behaviour 

r((rf + 1)/2) 

Eliminating pr between ( P29| ) and ( p31| ) yields 

, _ vr r((rf-l)/2)r((rf+l)/2) 
"2 (^-^)^ 

valid at 'short' distances (i.e. ^ ^ r ^ L{t)). Let us compare this exact result with 
( |18CI|) , obtained using gaussian auxiliary field methods. For short distances, when 
C i s clo se to unity, (|18CI|) becomes C^^ = (vr/2) (1 — C). This has the same form 



as (|232|) , but with a different coefficient of (1 — C*). However, the exact coefficient 
approaches the 'gaussian' value of 7r/2 in the limit d 00, consistent with our 
argument that the auxiliary field m is indeed gaussian in this limit. 

6.4.2 Vector Fields 

For vector fields, 0^ saturates to u nity only as an inverse power of the distance from 



a def ect, and the representation (|223|) is no longer appropriate. Instead, ( ^91) (see 
also (|182D ) implies 

1 - 0^ ^ ^V^^ (233) 
when the distance r from the core satisfies ^ ^ r ^ L{t). The calculation of the 
singular part of , due to the presence of defects, then follows that of the usual 
pair correlation function C: 



CLn, = eP.ef J d-X ( ^^ ) , (234) 
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where is the usual component of r in the plane perpendicular to the defect, 
and the angled brackets indicate an isotropic average over the orientations of r. 
Evaluating the x-integral first gives 

r^Af - n/2 r(2 - n/2)r^(n/2 - 1) 4 



Using ( |215|) , with a = n — A, gives (after some algebra) 

r(n/2- l)r(l-n/2)r(rf/2) 

r((rf + n-4)/2) 



casing = Pdef^ ' t^^^^ , ^ ^^7^:^^ • (236) 



This result gives the singular pait of for all n > 2. The poles in the numerator 
at even n signal additional factors of Inr (actually ln(r/L)), as in the calculation of 
C. The amplitude of the logarithm can be extracted by setting n = 2m + e, with m 
an integer, letting e go to zero, and picking up the term of order unity. This gives, 
for even n greater than 2, 

Q"... - -4?V.. (--)"^ ^„ ^2)ra.Ti-4)/2) -""'- ■ P^^) 

This singular short- distance behaviour implies, as usual, a power-law tail in 
Fourier space, S4(k,t) ~ j^-jd +n-i ) n > 2, in agreement with the prediction of 
the gaussian theories (section |5.3| and |^). Again, there is a heuristic argument 
for this [Q. The result ((1 — 0^)) ~ for n > 2 suggests the scaling form 

= L^^f{r/L), with Fourier transform L'^~^g{kL). Extracting the expected 
proportionality to pdef ~ L~'^ for kL ^ 1 generates the required ^^^jj^ 

The case n = 2 is more complicated, due to an additional pole in ( |236[ ). The 
same technique, however, can be used. We put n = 2 + e in ( ^3(j| ), let e go to zero, 
and pick up the term of order unity. The result is 

CLng = ^'Pdef rr{d - 2) In^ r/r' , d>n = 2. (238) 



For d = 2, we must set = 2 in ( |236| ) before taking e to zero. This gives 

~^ising = '^'^i'^Pdef 



CLn, = '2<'Pdef \nr/r\ d = n = 2. (239) 



The corresponding structure factor (the Fourier transform of C^^^^) has a large- 
momentum tail of the form ^'^Pdef ^.y? {k^) / k'^"'^ ■ These results are discussed in more 
detail in reference [ p.05|| . In particular, it is shown that the result for c? = n = 2 is 
inconsistent with conventional scaling form. 



6.4.3 Defect-Defect Correlations 

The short-distance behaviour of the defect-defect correlation functions, introduced 



in section p.3.3| , may also be determined exactly, at least for extended defects. The 
argument is a si mple extension of that used to calculate C4 for scalar fields. From 
the first part of ( |187| ), the correlator (p(l)p(2)) is just the joint probability density 
P(0,0) for points and '2' both lying on a defect. Clearly P(0, 0) = pde/P(2|l), 
in the usual notation. But for r <^ L, P(2|l) is dominated by cases where '1' and 
'2' are in the same defect (pr ovid ed the defects are spatially extended, i.e. n < d). 
An obvious generalisation of (^27]) gives 



(240) 



(p(l)p(2)) = pde/' 



Sdr^ 

Pdef T{d/2) 



W^r" T{{d-n)/2) 
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This result differs from ( |188| ), obtained by using the gaussian approximation, but 
approa ches it in the hmit d oo, as we have by now come to expect. In the exact 
result (^401) , the amplitude of the divergence vanishes for point defects (n = d ), 



an important physical feature that is missing from the gaussian approximation (|18 



6.5 THE PROBABILITY DISTRIBUTION P(m(l),m(2)) 

To conclude this section on short-distance behaviour, we compute the exact form of 
P(m(l),m(2)) for scalar fields, valid for |m(l)|, |^(2)| and r all much smaller than 
L(t), and show explicitly that it is not gaussian. Technically, the regime in which 
are working corresponds to taking the limit L{t) oo with |m(l)|, |m(2)| and r 
fixed but arbitrary. This situation is illustrated in Figure 23, where the domain 
wall can be regarded as flat in the limit of interest. The identity P(m(l),m(2)) = 
P(m(l)) P(m(2)|m(l)), where P{x\y) indicates a conditional probability, gives 



P(m(l),m(2)) = p(5(m(2) -m(l) -rcos^)) 



(241) 



where the angled brackets indicate an isotropic average over 9, and we have used 
P(m(l)) = p, the wall density, for <^ L. Carrying out the angular average 

(with weight proportional to sin''"^ 6) gives 



P(m(l),m(2)) 



B 



'I d-V 



2' 2 



P 

r 



1 - 



(m(2) -m(l))2y^-^)/' 



(242) 



for |m(l) — m(2)| < r, and P = otherwise. Clearly P(m(l), m(2)) is not gaussian 
for general d. Howe ver, it approaches a gaussian for large d, as we now show. 

It is clear from ( p42| ) that, in the limit d oo, P(m(l),m(2)) vanishes except 
when |m(l) — m(2)|/r is of order 1/ ^d. Therefore we define 



A = \fd{m{l) - m(2))/r . 



(243) 



Now the large-d limit can be taken at fixed A. Taking d large in the beta function 
too gives 



P(m(l),m(2))-. ("1- 



1/2 



- exp( 
r 



-AV2) 



oo 



(244) 



Let's compare this result with that of t he systematic approach, which we argued 
is exact for large d, by evaluating ( |123| ) in the same limit. Equation (|158|) gives 
5*0 = 4t/d and 7 = exp(— r^/8t). Inserting these in ( |123D , and taking the limit 
t — s> 00 with r, m(l) and m(2) held fixed, gives 



P(m(l),m(2)) 



d 



exp(-AV2) 



(245) 



Eqs. (HI) and ([245|) a gree if p = ((i/STrt)^/^ which is just Eq. (p^ . We conclude 
that exact result ( P42D is consistent with the gaussian approximation in the limit 
d — > 00, but not for any finite d. 



7 GROWTH LAWS 

The exact short-distance singularities derived in the previous section, together with 
the scaling hypothesis, provide a basis for deriving exact growth laws for all phase- 
ordering systems with purely dissipative dynamics. 
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Although the growth laws for both nonconserved and conserved scalar systems, 
and conserved fields in general, have been derived by a number of methods, there 
has up until now been no simple, general technique for obtaining L{t). In particular, 
the growth laws for non-conserved vector fields have, until recently, been somewhat 
problematical. Here we describe a very general approach, recently developed by 
Bray and Rutenberg (BR) |Q, to obtain L(t) consistently by comparing the global 
rate of energy change to the energy dissipation from the local evolution of the order 
parameter. This method allows the explicit derivation of growth laws for 0{n) 
models, but the results can be also be applied to other systems with similar defect 
structures. 

The BR approach is based on the dissipation of energy that occurs as the system 
relaxes towards its ground state. The energy dissipation is evaluated by considering 
the motion of topological defects, when they exist. The defect contribution either 
dominates the dissipation or gives a contribution that scales with time in the same 
way as the total dissipation. The global rate of energy change, computed from the 
time derivative of the total energy, is equated to the energy dissipation from the 
local evolution of the order parameter. For systems with a single characteristic scale 
L{t), this approach self-consistently determines the time-dependence of L{t). 

7.1 A USEFUL IDENTITY 

We begin by writing down the equation of motion for the Fourier components (p^. 

dt$^ = -k^{dF/d$,^), (246) 

The conventional non-conserved (model A) and conserved (model B) cases are = 
and /i = 2, respectively. (Elsewhere in this article, the symbol /i has been used for 
the chemical potential: the meaning should be clear from the context). 

Integrating the rate o f ene rgy dissipation from each Fourier mode, and then using 
the equation of motion ( p46|) , we find 

de/dt = j^{{dF/d$^)-dA) 

= - J^k-^ (5t4-5*0-k) , (247) 

where e = (F) /V is the mean energy density, and J^^ is the momentum integral 
/ d'^k / {2t[Y . We will relate the scaling behaviour of both sides of ( p47| ) to that of 
appropriate integrals over the structure factor, ^(k, t), and its two-time generali- 
sation. Either the integrals converge, and the dependence on the scale L{t) can 
be extracted using the scaling form (0) (or its two-time generalisation (|^)), or the 
integrals diverge in the ultraviolet (UVj and have to be cut off at kmax ~ l/'C? corre- 
sponding to a dominant contribution from the core scale. It is just this small-scale 
structure that is responsible for the generalised Porod law (|70D for the structure 
factor, and the time-dependence of any integrals controlled by the core scale can be 
extracted from a knowledge of the defect structure. 

7.1.1 The Energy Integral 

To see how this works, we first calculate the scaling behaviour of the energy density, 
e, which is captured by that of the gradient term in (|I|): 

e ~ ((V0)2) 

k'^L'^g{kL), (248) 

where we have used the scaling form (^) for the structure factor. For n > 2 the 
integral in UV convergent, and a change of variables yields e ~ L~^. For n < 2, 
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when the integral in UV divergent, we use Porod's law (^) and impose a UV cutoff 
at ~ 1/^, to obtain p8[] 



L-^ ln(L/0 



n <2 , 
n = 2 , 
n>2 . 



We see that the energy is dominated by the defect core density, pdef 
n < 2, by the defect field at all length scales between ^ and L for n - 
variations of the order parameter at scale L{t) for n > 2. 



(249) 

- for 
2, and by 



7.1.2 The Dissipation Integral 



We now attempt to evaluate the right side of ( p47| ) in a similar way. Using the 
scaling hypothesis for the two-time function. 



(250) 



which is the spatial Fourier transform of (P), we find 

92 



dtdf 



t=t' 
h{kL) 



(251) 



where L = dL/dt. 

When the momentum integral on the right of (p47D in UV convergent we obtain, 
using ( |251|) , de/dt ~ —U'L^~'^. If, however, the integral is UV divergent, it will 



be dominated by the behaviour of the integra nd n ear the upper limit, so we need 
to know the form of the scaling function h in (|251| ) for kL ^ 1. It turns out that, 
in general, the large-kL form is quite complicated, with many different cases to 
consider |2^. However, we only need the result for those cases where the dissipation 
integral requires a UV cut-off, otherwise simple power counting is sufficient. For 
those cases, one additional assumption, which can be verified a posteriori, yields a 
simple and rather general result (Eq. ( ^59| ) below). 



7.2 EVALUATING THE DISSIPATION INTEGRAL 

7.2.1 An Illustrative Example 

To see what difficulties arise, and how to circumvent them, it is instructive to con- 
sider a scalar field. We want to calculate {dt4>i(_dt4>-k) in the limit kL ^ 1. It 
is clear that dt4> is appreciably different from zero only near interfaces. In fact, 
since d(j)/dt = in a. frame comoving with the interface, we have, near an interface, 
dt(f) = — V.V0, where v is the interface velocity. In real space, therefore, 

(9^0(1) 9*0(2)) = (v(l) ■ V0(1) v(2) . V0(2)) . (252) 

The large kL behaviour in Fourier space is obtained from the short-distance (r <^ L) 
behaviour in real space. For r <^ L, the points '1' and '2' must be close to the same 
interface. For a typical interface, with radius of curvature of order 1/L, the speed 
V is slowly varying along the interface. Furthermore, the interface may be regarded 
as 'flat' for the calculation of the short-distance correlation, just as in the derivation 
of Porod's law. It follows that the averages over the interface velocity and position 
can be carried out independently, giving 

{d^l) dtm) = (l/rf) {v') (V0(1) ■ V0(2)) . (253) 
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Fourier transforming this result gives 



~ (v^) /Lk'^-^ , A;L > 1 , (254) 

where the Porod result ( |69D was used in the final line. We will see that Eq. (|254|) 
requires a careful interpretation. 

The next step is to evaluate (f^). Since the characteristic interface velocity is 

L, we expect (f^) ~ L^. This assumes, however, that the average is dominated by 
'typical' values. This, as we shall see, is the key question. Consider a small spherical 
domain of radius r in a non-conserved system. The interface velocity (see Eq. (0)) is 
f ~ 1/r. For a conserved system, Eq. ( ^I]) gives f ~ Thus the relation v{r) ~ 

l/r^~^ where z = 2 and 3 for nonconserved and conserved systems respectively 
(and the growth law is L ~ t^^^) covers both cases. The fact that v diverges 
at small r raises the possibility that (f^) is dominated by small domains. The 
domain-size distribution function n(r) has the scaling form n{r) = L^^'^^^^ f{r / L) 
(in order that / drn{r) ~ L~''', consistent with scaling). The important small- a; 
behaviour of the function f{x) can be determined as follows. Consider a small 
time interval At. The domains that will have disappeared after this time interval 
are those with radius smaller than ta* ~ (At)^/^. The number of such domains 
is of order L"'^'^"'"^) /J^"^* dr f{r/L). The requirement that this be linear in At forces 
f{x) ~ x^~^ for X ^ 0. Using v ~ r~^^~^) we can estimate the contribution to (f ^) 
from short scales: 

„ jl' drr'^~^ n{r) v'^{r) 

{v ) ~ — T 

drr'^~^n{r) 

~ drr'^-''/ drr^+'-^ . (255) 

The integral in the denominator converges at short scales, giving a result of order 
j^d+z-i ^ For nonconserved fields {z = 2), the numerator converges for all d > 1, 
giving L'^~^ for the numerator, and (f^) ~ Since L ~ 1/L for this case, we 

have (f^) ~ as expected. For conserved fields {z = 3), however, the numerator 
only converges at short scales for d > 2. For those cases, one again finds (w^) ~ L^. 
For d = 2, though, the numerator is of order ln(L/^). This gives (f ^) ~ ln(L/^), 

i.e. there are c ontributions from all scales, and (f ^) ~ no longer holds. Putting 
this into ( p54| ) gives 

(<9t0k5t0-k) ~ \n{L/i)/L^k , fcL > 1 , (rf = 2, conserved) , (256) 

We will now show tha t th is result is wrong! 

The factor 1/L in ( p54| ) represents the total interfacial area density: ( |254|) i mplie s 
that interfaces contribute additively to (9i0k9t0_k)- In the derivation of Eq. ( p54| ) , 
however, we explicitly assumed that only interfaces of typical curvature, of order 
contribute. For a piece of interface of local curvature 1/-R, the condition that 
the interface be regarded as locally fiat on the scale 1/k requires kR 3> 1, not 

simply kL ^ 1. For fixed k^ ^/L, sharply curved interfaces, with i? ~ do not 
contribute to the Porod tail. This means that, as far as the computation of the large 
kL behaviour is concerned, there is an effective short- distance cut-off ai 1/k: only 
interfaces with radius of curvature 1/k should be included. For the calculation 
of the usual Porod tail in {(j)]^ 4>--k) this makes no difference, because interfaces with 

R ^ 1/k make a negligible contribution to the total interfacial area as kL — > oo. 
For the calculation of {dt4>\t_dt(f)-]s), however, it can make a big difference, because 
of the extra factor of v"^ inside the average. This means that, in evaluating (f^), 1/k 
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rather than ^ is the appropriate short-distance cut-off. Applying this to conserved 
scalar fields in d = 2 gives (f^) ~ \n{kL), and 



(9i0k<9t0_k) ~ \n{kL)/V>k , A;L > 1, {d = 2, conserved) , (257) 



instead of (|256D. 



The final step is to insert Eq. (E^) into the dissipation integral (|247D, with 



d = 2 = n. One immediately sees that the integral is UV convergent: the L- 
dependence can be extracted trivially by a change of variable, de/dt ~ —1/L^. So 
we did not actually need the asymptotic form of (<9t0k f?t0-k) after all (except to 
show that the dissipation integral is UV convergent)! Note that the final result, 
de/dt ~ — is consistent with e ~ 1/L (Eq. ( |249| ) with n = 1) and the result 
L ~ for conserved (i.e. /z = 2) scalar fields. 



There is, however, one last complication. Eq. ( P57|) is still not quite correct! 
This is because the expression v{r) ~ for the velocity of a small drop (i.e. with 
r <^ L) breaks down for d = 2 due to the singular nature of the Greens function for 
the Laplacian. For this case one finds instead v{r) ~ l/r^ln(L/r). Using this 



gives finally {dt(f)]idt4>-k) ~ ln[ln(fcL)]/L^fc for kL ^ 1, instead of (P57| ). This does 
not alter, of course, the conclusion that the dissipation integral is UV convergent, 
and can therefore be evaluated simply by a change of variable. 



7.2.2 The Way Forward 

We have gone through this one case in some detail, because we can extract from 
it a general principle that avoids treating every case separately (although this can 
be done |^). The central point, given extra emphasis by the discussion above, is 

that we only need to know the asymptotics of (^dtcpk ■ dt$-\^ in those cases where 

the dissipation integral is UV divergent. The main result (with exceptions that can 
be enumerated) is that in all such cases the 'naive' estimates, obtained by using 

(f^) ~ L^, are correct. 

To make further progress we introduce the additional assumption, which can be 
checked a posteriori, that the dissipation is dominated by the motion of defect struc- 
tures of 'characteristic scale' L{t). By the 'characteristic scale' we mean the typical 
radius of curvature for extended defects (n < d), or the typical defect separation 
for point defects {n = d). That is, we are assuming that the dissipation is domi- 
nated by the motion of typical defect structures, and not by the disappearance of 
small domain bubbles, small vortex loops, or by the annihilation of defect-antidefect 
pairs. If the latter were true, the dissipation would be dominated by structure at 
the core scale, and the arguments given below would fail. We recall that for the 
case d = 2 = fi discussed above, the final k-integral for the dissipation was con- 
vergent, implying that the dominant /c- values are order 1/L, and the worries about 
the possible importance of small-scale structure were ultimately groundless. If the 

final k-integral were UV divergent, and the large kL limit of (dt$]i ■ dtcj)-]^ had 

important contributions from short scales, then the dissipation would be dominated 
by structure at the core scale, violating our assumption. Therefore, when our as- 
sumption holds, either (dt(j)\i ■ dt(j)--k^ is dominated by defect structures of scale L, 

or the final integral is UV convergent, or both. 

For t he re quired cases where the final integral is UV divergent, the large- /cL(t) 
limit of ( |251|) can be extracted from the physical/geometrical arguments used to 
obtain the generalised Porod law ([70|). According to our assumption, we can treat 
the defects as locally flat (or well separated, for point defects) for kL ^ 1. From 
( P51|) , we are interested in the behaviour of the two-time structure factor, ^(k, t, t') = 

(^(fk{t) ■ 0_k(i'))5 limit that the two times are close together. In the limit 

kL ^ 1, this will be proportional to the total density L^'"- of defect core. Introducing 
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L = {Lit) + Lit')) 12 and A = - L{t'))/2, we obtain 



Mt) ■ 0-k(t')) ~ jjTk^ ciikA) , A;L » 1 , (258) 

wher e consistency with Porod's law for t = t' requires a(0) = const. Using this in 
( P51|) gives 

■ ~ ^^^^^^ , fcL » 1 . (259) 

This reduces to ( p54| ) for n = 1 (with (f^) ~ L^). It should be stressed that 
we are not claiming that (|259|) is a gener al res ult, only that it is valid when we 



need it, i.e. when the dissipation integral ( 247 ) requires a UV cut-off. There are 
three possibilities: (i) The integral is UV convergent, its dependence on L{t) can be 

extracted by a change of variable, and the large-/cL behaviour of (dt(j)\<i ■ dt(j)--kj is not 

required, (ii) The integral is UV divergent, but the dissipation is still dominated by 
structures of scale L{t). Then we can use ( |259| ). (iii) The dissipation has significant 



contributions from structures with local curvature (or spacing) of order the core 
scale. Then one cannot treat the contributions from different defect core elements 
as independent, (|259|) no longer holds, and the present approach is not useful. For 
the moment we will proceed on the assumption that (i) or (ii) obtain. We will show 
that these possibilities cover nearly all cases. Examples of when (iii) holds will also 
be given. These include the physically interesting case d = n = 2. 



7.3 RESULTS 



Putting ( P59| ) into the dissipation integral (|247| ) shows that the integral is UV con- 



vergent for kL ^ 1 when n + fx > 2. Otherwise the integral is dominated by k near 
the upper cut-off 1/^. This gives 

J^k-f^ (dA-dt$-k) ~ Ln-^C-"'-^ n + /i<2, 

~ L'^L-'' ln(L/0 , n + /i = 2 , 

~ L^ Lf"-^ , n + fi> 2 . (260) 

The final step is to equate the dissip ation rate ( p60| ) to the time derivative of 
the energy density ( |249| ) , as required by (|247|) , and solve for L{t). The results are 



summarized in Figure 24, as a function of n and /x, for systems with purely short- 
ranged interactions. The two straight lines separating regimes of different behaviour 
are the lines n = 2 and n + /i = 2 at which the energy and dissipation integrals 
change their form. Note that conservation of the order parameter (which applies in 
a global sense for any /i > 0) is irrelevant to the growth law for fi < 2 — n, where n 
is treated here as a continuous variable. At the marginal values, logarithmic factors 
are introduced. The growth laws obtained are independent of the spatial dimension 
d of the system. 

For non-conserved fields (/i = 0), we find L ~ t^/^ for all systems (with d > n 
or n > 2). Lead ing corrections in the n = 2 case are interesting: the InL factors in 
( p49|) and ( p60| ) will in general have different effective cutoffs, of order the core size 
^. This leads to a logarithmic correction to scaling, L ~ t^/^(l -|- 0(1/ Int)), and 
may account for the smaller exponent (~ 0.45) seen in simulations of 0(2) systems 
| |127| , |128| , |105||. Note that for nonconserved scalar fields, the energy ( |249| ) and the 



dissipation (£6^) have the same dependence on the core size ^ (i.e. both contain a 

factor ^~^), so this dependence cancels from L(t). The fact that the correct t^^'^ 
growth is obtained from naive power counting on the linear terms in the equation 
of motion should therefore be regarded as fortuitous. For example, with long-range 
interactions, this 'cancellation of errors' no longer occu rs, a nd naive power counting 
gives an incorrect result for nonconserved scalar fields ||129| , 
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For conserved fields {fi > 0) our results agree with an earlier RG analysis 0, |T0 
with additional logarithmic factors for the marginal cases n = 2 and n + /i = 2. 
Note that the conservation law is only relevant for n + n > 2. Therefore for vector 
fields (n > 2), any /i > is sufficient to change the growth law, while for scalar 
fields {n = ij the conservation law is irrelevant for /i < 1, in agreement with the 
RG analysis [|I0|"| and earlier work of Onuki ||130| . 

Siegert anoRao |]76| have performed extensive simulations for n = 2, d = 3 and 
fi = 2. In their original paper they fitted L{t) to a simple power, and obtained 
a growth exponent slightly larger than 1/4. Recently, however, Siegert has shown 
that a very much better fit is obtained when the predicted logarithmic correction is 
included ||131|. 



7.3.1 Exceptional Cases: n = d <2 

In what cases is our key assumption, that dissipation is dominated by the motion 
of defect struct ures of characteristic scale L, correct? Certainly for any n > 2, the 
energy density ( p49| ) itself, and hence dissipation, is dominated by variations at scale 
L{t). Therefore, we limit the discussion to the case n < 2. 

For n < 2, the energy density is pro port ional to the defect core volume (with 
an extra factor ln(L/^) for n = 2, see (|249|) ) , but we will show that, in general, 
dissipation is still dominated by defect structures with length scales of order L. To 
see this, we investigate the contribution to the energy dissipation from small-scale 
structures (e.g. small domains, vortex loops, or defect-antidefect pairs): 

POD 

de/dt = dt dln{l,t)e{l) 
= - dlda{l,t)e{l) 

= j(Oe(0+ / dlj{l)de/dl, (261) 



where n{l,t) is the number density of defect features of scale /, e{l) ~ is the 
energy of a defect feature (with an extra ln(//^) factor for n = 2), and j{l,t) is the 
number fiux of defect featur es. W e have used the continuity equation, dtn + dij = 
to obtain the second line of ( |261|) , and the t-dependence has been suppressed in the 
final line. The total number of defect features, N, scales as ~ L~'^, and so N 
does not change significantly over times smaller than L/L. Since defects only vanish 
at the core scale, we have = j(0- It follows has a finite, non-zero, short- 
distance limit of order N ~ —L/L'^^^. We can use this to examine the convergence 
(at short-distance) of the final integral in ( |261| ). 

For d > n, the integral in ( |261j ) is well-behaved at small I, because e(/) ~ 
(x ln(//^) for n = 2), and the integral dominates the J (0^(0 term. The integral 
can be estimated by setting ~ j(^) and introducing a large-distance cut-off at 

/ ~ L. This gives de/dt ~ jXO^"^"" (x ln(L/0 for n = 2). This is 

just what one gets from differentiating ( [^49| ) (for the cases n < 2 considered here), 
verifying the consistency of the calculation. 

For d = n and n < 2, however, e{l) ~ constant, since the (point) defects only 
interact weakly through the tails of the defect profile. (The one physical example 
is the d = 1 scalar system). The leading contribution to the energy of a defect 
pair is just the core ener gy o f the individual defects, and dissipation is dominated 
by the J (0^(0 term in ( p61| ), which describes defect pairs annihilating. Since the 
dissipation occurs at separations Z ~ ^ <C L, the derivation of ( 259 ) no longer holds. 
In fact since, the energy of a defect pair depends only weakly on the separation for 
/ ^ L, the system will be disordered, with an equilibrium density of defects at any 
non-zero temperature. At T = 0, we expect slow growth that depends on the details 
of the potential V{(f)) ||132|| . These cases, including the d = 1 scalar system, are at 
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their lower critical d imen sion, and are beyond the scope of the simplified approach 
presented here (see [ [132|| for a fuller discussion). 

The 2d planar system (n = = 2) is a special case. The logarithm in the energy 
of a vortex pair, e ~ ln(//^), leads to a logarithmically divergent integral in ( p6 1] ) , 
i.e. vortex pairs with separations bet ween ^ and L contribute significantly to the 
energy dissipation. In this case Eq. ( |259| ), which depends on the kL ^ 1 limit 
being a single defect property, is again questionable. As a result the present method 
cannot address this case. Indeed, the contributions to the dissipation from all length 
scales suggest a possible breakdown of scaling. 

7.3.2 Systems Without Defects 

Since systems without topological defects {n > d + 1) will have convergent momen- 
tum integrals for kL ^ 1, we obtain L ~ t^/C^+t^) for these cases. We can also 
apply this result to systems with topological textures {n = d + 1), even though the 
appropriate Porod's law is not known. Since defects with n > d must be spatially 
extended and without a core, they will have a smaller large-/c tail to their structure 
factor S'(k, t) than any defects with cores. So for n > 2, w hen the energy dissi- 
pation clearly occurs at length scales of order L(t) (see (|249|) ) and the momentum 



integrals for defects with cores converge, our results should apply (|^ contains a 
fuller discussion of this point). Consequently the results in Figure 24 will apply for 
any system, apart from those systems explicitly excluded above. 

Of course, all this is subject to the caveat that the two-time structure factor 
exhibits the scaling form (|251|), on which the whole of this section is built. One 



explicit counterexample is the d = 1, n = 2 system discussed in section [4.5| , for 
which ( p51| ) explicitly fails. As a result, the growth law characterising equal-time 
correlations is not L ~ t ' , as suggested by Figure 24, but L ~ t^/^ In fact one use 



the two-time result ( p.l3| ) to calculate (^c?t0k ■ dt(j)-]ij explicitly for this system |]78 
and show that it is consistent with L ~ t^^'^ growth. 

7.3.3 Other Systems 

The strength of this approach is that it can be applied to systems with more com- 
plicated order parameters than ra-component vectors, provided they have purely 
dissipative dynamics. Then an equation of the form ( |247|) ca n be written down. 
The details of the energy functional (|1|) are unimportant ||13^ . The important in- 
gredients are the existence of an 'elastic energy', associated with spatial gradients 
of the order parameter, the conservation law (if any), characterised by fi, and the 
defect structure if any. The derivation is independent of the initial conditions, and 
so, e.g., applies equally to critical and off-critical quenches as long as the system 
scales at late times. We simply choose a Porod's law ( p^ ) to represent the dominant 
defect type, which is the one responsible for the asymptotic tails of the structure 
factor scaling function, i.e. the one with the smallest 'n' . Wh en the energy density 
is domin ated by defects, i.e. when the energy integral (|248|) is UV divergent, the 
relation ( |24S ) between the energy density and the structure factor, shows that the 
'dominant' defects will also be the ones which dominate the energy density. As 
examples, we consider nematic liquid crystals and Potts models. 

In bulk nematic liquid crystals, the 'dominant' defects (in the above sense) are 
strings, giving a Porod tail of (|69|) with n = 2, which with no conservation law 

implies L ~ t^/^, consistent with recent experiments [^2[ ^ and simulations 



The g-state Potts model has q equivalent equilibrium phases, which give rise 
to q{q — l)/2 different types of domain wall. These can indexed ajS, where a, 
(3 = 1, . . . g, are the phases separated by the wall. Three domain walls of type a/?, 
a7 and /57 can meet at a point {d = 2) or line (rf = 3), which represents a new type 
of defect. It is clear, however, that the Porod tail and energy density are dominated 
by the walls, so that the Potts model behaves as an n = 1 system. As a result, 
L{t) ~ t^/^ and t^^^ for nonconserved and conserved order parameter respectively. 
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Recent numerical results ||134| , [135|| support these predictions, after initial suggestions 

that the growth was slower. The t^^^ growth for conserved systems is also predicted 
by the Renormalisation Group approach of section 

It should be emphasized that the classification of nematic liquid crystals and 
Potts models as 'n = 2-like' and 'n = 1-like' respectively, pertain only to the Porod 
tails and the growth laws. As far as scaling functions (e.g. for pair correlations) are 
concerned, these systems belong to their own universality classes. Similarly, for off- 
critical quenches of conserved systems, the growth law is independent of the volume 
fractions of the phases, but the scaling functions are not. 



8 RENORMALIZATION GROUP RESULTS 

As with any other scaling phenomenon, it is tempting to try to apply Renormalisa- 
tion Group (RG) concepts to the late stages of phase ordering. The basic idea is to 
associate the scaling behaviour with a fixed point of the equation of motion under 
a RG procedure consisting of a coarse-graining step combined with a simultaneous 
rescaling of length and time. Such a procedure, if successful, would indeed provide 
a first principles derivation of the scaling behaviour itself, which has, up to now, 
been lacking (except for specific soluble models discussed in section^). 

Underlying such an idea is the schematic RG flow for the temperature, depicted 
in Figure 4. The critical point corresponds to a fixed point of the RG transfor- 
mation. At temperatures above (below) Tc, coarse graining the system leads to a 
system which is more disordered (ordered), corresponding to a system at a higher 
(lower) temperature. This schematic flow is indicated by the arrows in Figure 4. It 
follows from this that a quench from any T > to any T < should give the 
same asymptotic scaling behaviour. Any short-range correlations present at the ini- 
tial temperature will become irrelevant when L{t) ^ where is the correlation 
length for the initial condition. A different universality class is obtained, however, 
when the initial condition contains sufficiently long-range (power-law) spatial cor- 
relations, e.g. following a quench from Tc. Such cases will be discussed below. It 
follows from the previous section that the initial conditions cannot affect the growth 
law (provided that they still yield scaling behaviour). 

According to Figure 4, asymptotic scaling is controlled by the zero-temperature 
(or strong coupling) fixed point, justifying the neglect of thermal noise in the equa- 
tions of motion. We will see below how this works out in practice. A classific ation 
of systems according to the role of thermal noise has been given by Lai et al.||136 



In some systems, such as the Cahn-Hilliard (CH) systems considered in this section, 
thermal noise can simply be neglected. For kinetic Ising models (with conserved 
dynamics), where freezing occurs for T strictly zero, the temperature modifies the 
bare transport coefficient A, but in a scale-independent way that does not change the 
growth law. In systems with quenched disorder, however, there is a scale- dependent 
renormalisation of the kinetic coefficients that leads to logarithmic growth with T- 
dependent amplitudes. This case will be discussed in detail in section |8.3.2| 



The idea of using RG methods in this context is not new (see for example 



136|, 
y m 



13711 ), but in practical applications the RG framework has be en ex ploited most 
the numerical context via, for example, the Monte Carlo RG ||138|| . 

The difficulty with applying the RG to phase ordering is that, due to the absence 
of a convenient small parameter, analogous to e = 4 — d for critical phenomena, 
one cannot obtain explicit RG recursion relations. However, one can still make 
some progress. For conserved fields, a very simple and general result taken over 
from critical phenomena can be used to determine growth exponents exactly 10 



without the need to construct explicit RG recursion relations for the entire set 
of parameters specifying the equation of motion. Without such explicit recursion 
relations, of course, the very existence of a fixed point has to be taken on trust. 
This is tantamount to assuming the validity of the scaling hypothesis ab initio, and 
inferring the existence of an underlying RG fixed point. This is the approach we 
will adopt. It will take us surprisingly far. 
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8.1 THE RG PROCEDURE 



8.1.1 Equation of Motion 

We start by recalling the Cahn-Hilliard equation for a conserved order parameter, 
generalised to vector fields. Introducing the transport coefficient A explicitly on the 
right-hand side (we have previously absorbed A into the timescale) gives 

d$/dt = A V'^{6F/6(f) . (262) 

Next we Fourier transform, and divide through by Afc^. For generality, and antici- 
pating future requirements, we will write the equation in the form 

+ il^ = -^+aW. (263) 



Xk^' Tj dt 

Here we have replaced fc^ by the more general k^^, as in section 0, and included on the 

left-hand side a term {l/r)d(j)k/dt appropriate to non-conserved dynamics, which 
is recovered in the limit A —>■ oo . For any finite A, however, the order parameter is 
conserved by the dynamics (|263| ). We include the extra term because it will in any 
case be generated (along with many other terms) after one step of the RG procedure. 

A gaussian white noise term C,]i{t), representing thermal noise, has also be in- 
cluded in (|263|) . We require that the canonical distribution be recovered in equilib- 
rium, i.e. P[(j)] oc exp(— F[0]/T). The usual fiuctuation-dissipation relation fixes the 
noise correlator, 

mh)e^{h))=2T6,,6{h-t,) + , (264) 

where i, j = 1, . . . ,n indicate Cartesian components in the internal space. We have 
argued previously that thermal noise is irrelevant. The RG approach shows this 
explicitly. 



8.1.2 The Coarse-Graining Step 

One RG step consists of the following four stages: (i) The Fourier comp onents 0k (^) 
for the 'hard' modes with A/b < k < A are eliminated by solving (p63|) for the time 
evolution of these modes, and substituting the solution into the equation of motion 
for the remaining 'soft' modes with k < A/b. Here A ~ 1/^ is a UV cut-off, and b 
is the RG rescaling factor. 

(ii) A scale change is made, via the change of variable k = k'/b, in order to 
reinstate the UV cut-off for the soft modes to its original value A. Additionally, time 
is rescaled via t = Vt' . The requirement, imposed by the scaling hypothesis, that 
the domain morphology is invariant under this procedure, fixes z as the reciprocal 

of the growth exponent, i.e. L{t) ~ t^/^. Finally, the field 0k(^) for k < A/b is 
rewritten as 

k{t)=k'/b{hH')=b^$'^,{t') . (265) 

The scaling form (|^ for the structure factor becomes ^(k,t) = t'^'^g{kt^/^). From 
the definition of S and equation ( |265|) , 

^(k,t) = 62^(0;,,(t') ■ 0U,(t')) = b''^t"''-''g{k't'^'') = b^^-H""/' g{kt^/') , (266) 

from which we identify ( = d/2. 

(iii) The new equation of motion for the soft modes is interpreted in terms of a 
rescaled transport coefficient A' and free energy F'. In addition, terms not originally 
present in ( |263[ ) will be generated, and must be included in subsequent RG steps. 
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Similarly, one must allow for a more general structure for the thermal noise than 

( |^64|) ■ Finally, the distribution Po({0k(O)} of initial conditions will also be modified 
by the coarse graining. 

(iv) Scaling behaviour is associated with a fixed point in which both the equation 
of motion and Pq are invariant under the RG procedure. In particular the fixed- 
point free energy is the that appropriate to the 'strong-coupling' fixed point, which 
is attractive for systems below Tc- Note also that the fixed distribution Pq contains 
the scaling morphology. 

8.1.3 RG Recursion Relations 

Unfortunately, the above procedure cannot be carried out explicitly, due to the 
absence of a small parameter [|139|| , and remains largely a 'gedanken RG'. Neverthe- 
less, on the assumption that a fixed point exists (equivalent to assuming scaling), 
the recursion relations for the transport coefficient A and the temperature T can 
be written down exactly. This is sufficient to determine z and to test the stability 
of the nonconserved fixed point against the conservation constraint. These results 
agree with those of section 0. In addition, however, we can also identify universal- 
ity classes, clarify the role of the initial conditions in determining the large-scale 
structure, and make strong predictions about the effects of quenched disorder. All 
of these are beyond the scope of the methods of section |^. In this sense, the two 
approaches are complementary. 

The observation that enables further progress is that the term in the 

equation of motion ( p63|) is singular at /c = 0. Since no new large-distance singular- 
ities can be introduced by the elimination of small length-scale degrees of freedom, 
it follows that the coarse- grain ing step (i) of the RG does not contribute to the 
renormalisation of 1/A in (By contrast it can, and does , contribute to the 

renormalisation of l/F UTof, which is a non-singular term in ( p63| )). As a result, 1/A 
is change d on ly by the rescaling step (ii). Exactly the same argument at the critical 
point [§, |140|| leads to the identity z = A — i] between the dynamic critical exponent 



z and the static critical exponent rj. It is important to recognize that the latter 
result is non-perturbative, and is not restricted to the conventional Wilson-Fisher 
fixed point. 

Since the strong-coupling fixed point is attractive (see Figure 4), the free-energy 
functional scales up at this fixed point, -F[{0k7b}] = F[{4>k'}]: where the exponent 
y can be determined by elementary arguments. Using this and ( |265| ) in ( |263| ) gives 
the coarse-grained equation of motion 

+ ''-'f II + -I + ••■)# + • - £ + . (267) 

where . . . indicates additional terms generated by the coarse graining step. Dividing 
through by to restore the right-hand side to its previous form, gives 

("-'--^ + f n + ■ • -1 + ■ • •) # ^ ■ - ^ «^'(''' • (^^«' 

where the new noise term is 

^Ut')=b^''^^'/t{bH') , (269) 

with correlator 

(d^(t'i)a'(i2)) = b''^-''-' 2T6,^ - Q (6^3^ + ^ [1 + ■•■]) . (270) 
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The absence of contributions to the terms, either in the equation of motion 

or the noise correlator, from the coarse graining step, means that the recursion 
relations for 1/ A and T can be written down exactly: 

JJ2(+^.-y-z ^ (271) 



XJ VA. 
r = b-^T . (272) 

The T-equation is just what one would expect. Since the free energy functional 
scales up as b^ at the strong-coupling fixed point, rewriting the equation of motion, 
as in ( ^68| ), in a form in which the free energy functional is unchanged is equivalent to 
scaling temperature down by a factor b~y. At the same time, the transport coefficient 
A renormalises as in ( ^71| ). This last equation determines the growth exponent for 
all cases in which the conservation constraint is relevant, in a sense to be clarified 
below. 



8.2 FIXED POINTS AND EXPONENTS 

In the strong-coupling phase (i.e. for T < T^), T fl ows to zero under repeated 
iteration of the RG procedure, implying y > in (|272|) . If the dynamical flxed point 
controlling the late-stage scaling regime is described by a non-zero value of 1/A (i.e. 
when the conservation law is relevant), then ( |271|) implies 

z = 2C + fi-y = d + fi-y , (273) 



where we inserted the value ( = d/2 at the last step. Eq. ( 273 ) is exact, given the 
scaling assumption underlying the RG treatment. 

It is interesting to consider the same argument at the critical flxed point. Then 
T' = T = Tc implies y = and z = 2( + fi. The structure factor scaling relation 
reads, for this case, 5'(k, t) = L?'~'^g{kL) = t^'^~^^^^g{kt^/^), so the analogue of Eq. 
( ^66| ) flxes C = (2 — r))/2, and z = 2 + fi — rj. For /i = 2, this is the familiar result 
z = 4 — 7 ? for model B |14U| | , which we stress is an exact, nonperturbative result. 
Eq. ( |273| ) is just the generalisation of this result to the strong-coupling flxed point. 

To determine y for the strong-coupling flxed point, we coarse-grain the system 
on the scale L{t). At this scale the system looks completely disordered, so the excess 
energy per degree of freedom (i.e. per volume L{tY) is of the order of the local excess 
energy density at that scale, i.e. of order L(ty. The excess energy density on the 
origi nal scale therefore decreases as e ~ L(ty^'^. Comparing this with the result 
( p49|) for e obtained in section |^ gives 

y = d — n , n <2 , 

= d-2 , n> 2 . (274) 

Note, however, the extra logarithm in (p49|) for n = 2. 



For the usual scalar (n = 1) and vector (n > 2) flelds, ( 274 ) gives the usual 
results y = d — 1 and y = d — 2 respectively, familiar from statics: b'^~^ is just 
the energy cost of a domain wall of linear dimension b, while b'^~^ is the energy 
cost of impos ing a slow twist of the vector fleld over a region of size b"^. The extra 
logarithm in (|249|) for n = 2 is due to the vortices, which dominate over slow 'spin- 
wave' variations for this case. As an amusing aside we note that (|274| ) gives the 



'lower critical dimension' di, below which long-range order is not possible for T > 0, 
for the continuation of the theory to real n. Since the existence of an ordered phase 
requires y > 0, we have di = n for n < 2 and di = 2 for n > 2. The result for n < 2 
recovers the known results for n = 1 (the scalar theory) and n = (the self-avoiding 

walk) . 

Inserting ( |274| ) into ( p73| ) gives the flnal result for z: 

z = n + , n <2 , 

= 2 + fi , n> 2 . (275) 
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These results agree with those derived in section ^ which are summarised in Figure 
24. For scalar model B (n = 1, /i = 2) we recover the usual t^^^ Lifshitz-Slyozov 
growth, while for vector model B with n > 2 we obtain t^^^ growth. At the crossover 
value n = 2, there is an extra logarithm that the RG method does not see (since 
it determines only the growth exponent), but which is captured by our previous 

approach (section |7 ). 

Comparison of ( 273| ) with Figure 24 shows that ( p73| ) is not valid below the line 
n + fi = 2. How do we see this within the RG approach? Recall that to derive 
( |^73| ) we have to assume that (1/A) is non-zero at the fixed point, i.e. that the 
conservation constraint is relevant, in the RG sense. Consider now t he fix ed point 
of the non-conserved system, with A = oo in the equation of motion ( p^63| ). Let the 
corresponding value of z be Znc- Now in trod uce the conservation law through an 
infinitesimal 1/A. The recursion relation (|271D then gives 



VA 



(276) 



wher e Zc = 2( + ^ — y is the value of z (Eq. (|273|) ) at the conserved fixed point. 
Eq. ( p76| ) shows immediately that (1/A) iterates to zero for Zc < Znc, i-e. the con- 



servation law is irrelevant when < Znc- Since L{t) ~ t^^^, this means that the 
conserved system cannot exhibit faster growth than the nonconserved system. This 
is intuitively reasonable: an additional constraint cannot speed up the dynamics. 
There is an interchange of stability of RG fixed points when Zc = z^c, the fixed 
point with th e lar ger z being the stabl e on e. It follows that z = max {z^, Znc), with 
Zc given by (p73|) in general and by (|275|) for the 0(n) model. Since z = 2 for 
the non-conserved 0(n) model (see 0), this interchange of stability accounts for the 
crossover line n + = 2 in Figure 24. (We note that the same reasoning implies a 
similar interchange of stability, and the result z = max [zc, Znc), at the critical fixed 
point [|141I] ). 

The generality of Eq. ( |273|) deserves emphasis. For any system with purely 
dissipative conserved dynamics, one only needs to insert the value of y, which can 



be determined from the energetics as in the derivation of (|274 ) for the 0{n) model. 
As an example, consider again the g-state Potts model. The energy density is 
dominated by domain walls, so the energy density scales as e ~ 1/L ~ L^"*^, giving 
y = d — 1 just as for the Ising model. Therefore, the usual Lifshitz-Slyozov t^^^ 
growth is obtained for /i = 2 . Of course, we already obtained this result in section 
0. Recent numerical studies [|135|| confirm this prediction. 

As a second example we note that, in agreement with our findings in section |^, 
the growth law in independent of the nature of the initial conditions (which played 
no role in the derivation), provided scaling is satisfied. A case of experimental 
interest is a conserved scalar field - the t^^^ Lifshitz-Slyozov growth is obtained for 
all volume fractions of the two phases. By contrast, the scaling functions can depend 
on the form of the initial conditions. This should not be too surprising, since the 
fixed point distribution for the initial conditions contains the scaling morphology. 



This will be discussed in detail in section B.3.1 



8.3 UNIVERSALITY CLASSES 

The present RG approach cannot, unfortunately, determine Znc, since it rests on 
(1/A) being non-zero at the fixed point. Neither can it explicitly pick up the loga- 
rithms on the boundary lines n = 2 and n + fi = 2 oi Figure 24. So does it have any 
advantages over the seemingly more powerful energy scaling approach of section 0? 
The answer is an unequivocal yes. The reason is that the RG identifies universality 
classes as well as exponents. As an example, consider a scalar n = 1 system with 
< 1. The energy scaling method tells us that L{t) ~ t^^^, as for nonconserved 
dynamics, but tells us nothing about correlation functions. The RG, by contrast, 
tells us that when the conservation is irrelevant not only the exponents but also all 
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correlation scaling functions are the same as those of the nonconserved system, i.e. 
for /i < 1 the scalar system is in the nonconserved universality class. 

At first this result seems paradoxical: in the scaling form for the structure factor, 
5'(k, t) = L'^gikL), the scaling function g{x) has a non-zero value at x = for 
nonconserved dynamics, whereas for conserved dynamics g{x) must vanish at x = 0. 
So how can a system with conserved dynamics be in the nonconserved universality 
class? To understand this one needs to remember that the scaling limit is defined 
by /c — >■ 0, L — >■ cxo, with kL fixed but arbitrary. Onuki has argued that, for /i < 1, 
5'(k, t) ~ k'^^L'^^'^ for — ^ ||130|| . If we imagine plotting S'(k, t) in scaling form, i.e. 
g{x) = L~'^5'(k, t) against x = kL, then Onuki's small-fc form gives a g{x) of order 
unity when x ~ L^^^^^, which is vanishingly small as L ^ oo for /i < 1. In other 
words, on a scaling plot the region of kL where the nonconserved scaling function 
is inaccurate shrinks to zero as L — cxd. 



8.3.1 The Role of the Initial Conditions 

To what extent do the scaling functions depend on the probability distribution 
for the initial conditions? The RG answers this question [|4|. New universality 
classes are obtained when sufficiently long-ranged (power-law) spatial correlations 
are present immediately after the quench. These could either arise 'physically', as in 
a quench from Tc, or be put in 'by hand' as initial conditions on the T = dynamics. 
Consider initial conditions with a gaussian probability distribution of variance 

{<Pm<P\m = A{k) 6., , (277) 



We recall the definitions, introduced in section |4.2| , of the response to and correlation 
with the initial condition: 

- (HI) 

C(k,t) = (0LW0-k(O)) (279) 

respectively, where C(k, t) is a shorthand for the two-time structure factor S(k, t, 0). 
The gaussian property of {0k (0)} means that these two functions are related by 

C{k,t) = A{k)G{k,t) , (280) 

a trivial generalisation of (^) that can be pr oved e asily using integration by parts. 
The RG treatment proceeds as in section 8.1.2[ The only additional feature is 



that the scaling form (p4[) for G(k, t) implies that the initial condition (/>k(0) acquires 
an anomalous scaling dimension related to the exponent A. (The exponent A should 
not be confused with the transport coefficient: the meaning should be clear from 
the context). Therefore we write, analogous to Eq. ( P65| ) for the rescaling of the 
field at late times, 

(O) = 4'/6(O)=&^0k'(O) (281) 



5k 



for the rescaling of the initial condition. This gives, analogous to ( 266 ) 



G{k,t) = b^-^ {^^) = b^-H"/' gnik't'^^n = h^-^-H"^ gn{kt''') , (282) 

from which we identify X = C ~ ^- The scaling of the equal time structure factor 
gives, as before, C, = d/2, so x = d/2 — A. 

Under the RG transformation, the correlator A{k) of the initial condition be- 
comes &^'^(0k'(O)0-k'(O)) = A{k'/b) + • • ■, where the dots indicate the contribu- 
tion from the coarse-graining step of the RG. So the new correlator, A'{k') = 
(0k'(O)0lk'(O)), is given by 

A'{k') = b^^"^ [A{k'/b) + ••■]. (283) 
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Consider now the case where A (A;) has a piece corresponding to long-range 
(power-law) correlations: 



A{k) = AsR + ALRk-'', (284) 

with < a < d. Then the real-space correlations decay as r^i^~'^). From (|283|) 
we can deduce the recursion relations for the short- and long-range parts of the 
correlator: 

^'sR = b''-'[^SR + ---], (285) 
Kr = b^^-''^'^ Alr . (286) 

Note that the long-range part A^r, being the coefficient of a si ngular (as k ^ 0) 
term, picks up no contributions from coarse-graining - Eq. (P86|) is exact. 



At the fixed point, both the equation of motion and the initial condition distri- 
bution must be invariant under the RG transformation. It follows from ( |286|) that 



the long-range correlations are irrelevant at the 'short-range fixed point' (i.e. A^r 
iterates to zero) if cr < cTc, where 

a, = d- 2XsR , (287) 

and XsR is the value of A for purely short-range correlations. When a > a^-, the 
invariance of Alr at the 'long-range fixed point' fixes A = Xlr = {d — a)/2, an 
exact result. Thus there is an exchange of stability of fixed points when a = ac- 
The determination of Xsr itself is non-tr ivial , since it requires explicit computation 
of the terms represented by the dots in (|285|) . 



For cr > cTg, the scaling behaviour belongs to a new universality class, in which 
the growth exponent is unchanged but the scaling functions, e.g. Qr^x), depend 
explicitly on a. Note that the function C(k, t), the correlation with the initial 
co nditio n, depends on a for any a > through the k~'^ term in the prefactor A{k) 
in (|280|) . Thus, in the scaling region, 

C{k,t)=ALRk-'^L^gRikL) . (288) 

Summing this over k gives the autocorrelation function: 

A{t) = (0(r, t) ■ 0(r, 0)) ~ L-'^d-a-x) _ ^289) 

In the 'long-range' regime, where X = {d — ct)/2, this gives A(t) ~ Con- 
sider, as an example, the 2d Ising model quenched from the equilibrium state at 
Tc. Then a = 2 — rj = 7/4. Measurements of A for the same model quenched from 
T = oo, with nonconserved dynamics, give Xsr — 0.75 [|19], as do experiments on 
twisted nematic liquid crystals, which are in the same universality class ||2T|. There- 
fore a > d — 2XsR and this system is in the 'long-range' universality classTlt follows 
that A(t) ~ ~ t"^/^^, which has been confirmed by numerical simulations 

mm- 

The scaling function g(x) for the equal-time structure-factor also has a dif- 
ferent form in the long-range regime. For any a > 0, ^(k, t) has a 'long-range' 
contribution S'ii?(k, t) varying as at small k. It is given by the diagram of 
Figure 25, where the circle represents A(A;) and the lines are exact response func- 
tions. Thus SLR{k,t) = ALRk-''G(k,t)G{-k,t). Using the scaling form (H) for 
G gives SLR(k,t) = AiRk~'^ L'^^[gR{kLyf' . Comparing this to the general scaling 
form ^(k, t) = L'^g{kL), we see that when a < ac (i.e. in the 'short-range' regime), 
Slr is negligible in the scaling limit (/c ^ 0, L — >• c>o, with kL fixed), and so does 
not contribute to the scaling function. In fact, since the long-range correlations are 
irrelevant in this case, the scaling function is identical to that for purely short-range 
correlations. 
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For a > ac, the contribution Sir survives in the scahng hmit, and the full 
scaling function is long-ranged, g{x) ~ for x — >■ 0. In real space, this means 
that the equal-time correlation function decays with the same power-law as the 
initial-condition correlator, i.e. C{r,t) ~ [L/rY'" for r ^ L. For the 2d Ising 
model quenched from T^, the predicted {L/rY^^ decay has been seen in simulations 



8.3.2 Systems with Quenched Disorder 

The influence of quenched disorder on the motion of interfaces and other defects is 
of considerable current interest in a variety of contexts. The new ingredient when 
quenched disorder is present is that the defects can become pinned in energetically 
favourable configurations. At T = this leads to a complete cessation of growth. 
For T > 0, thermal fluctuations can release the pins, but in general growth is much 
slower than in 'pure' systems, typically logarithmic in time. 

To see how logarithms arise, consider a single domain wall in a system with 
quenched random bonds. The typical transverse displacement of the wall over a 
length /, due to disorder roughening, is of order while the typical fluctuation of 
the wall energy around its mean value is of order These exponents are related by 



the scaling law ||142|| % = 2^-1-^ — 3, which can be obtained by estimating the elastic 
energy of the deformed wall as and noting that the pinning and elastic 

energies sh ould be compa r able . The barrier to domain motion can be estimated 
by arguing ||142| , |143| , |144| , |145|| that the walls move in sections of length /, where 
/ is the length scale at which the walls 'notice' their curvature, i.e. the disorder 
roughening should be comparable to the distortion, of order P/L(t), due to the 
curvature of walls with typical radius of curvature L(t). This gives / ~ and 
an activation barrier of order ~ /,x/(2-C) (assuming that the energy barriers scale 
in the same way as the energy fluctuations between local equilibrium positions of 
the wall). Equating this barrier to T gives a growth law 

L(t) ~ (Tint) (2-0/^ . (290) 



For d = 2, the exponents ( and x exactly known ||146|| : C = 2/3 and x = 
1/3, giving L{t) ~ (Tlnt)^. A number of a ttempts to measure L{t) in computer 
simulations have been made ||147| , |14(^ , |149|| , but it is difficult to obtain a large 
enough range of (Int)^ for a convincing test of the theoretical prediction. Recent 
experimental studies of the two-dimensional random-exchange Ising ferromagnet 
Rb2Cuo,s9Coo,nF4, however, suggest L{t) ~ (Int)i/^ with ip = 0.20 ± 0.05 ppj , 
consistent with the theoretical prediction ip = 1/A. 

Perhaps of greater interest than the growth law itself is the universality class for 
the scaling functions. It can be argued [l^ that, since L ^ I /,i/(2-C) for L ^ cxd 



(note that C < 1 for a system above its lower critical dimension, otherwise disorder- 
induced roughening would destroy the long-range order), on length scales of order 
L the driving force for domain growth is still the interface curvature: the pinning 
at smaller scales serves merely to provide the (scale-dependent) renormalization of 
the kinetic coefficient responsible for the logarithmic growth. This leads to the 
conclusion [|I^ that the scaling functions should be identical to those of the pure 
system, a prediction that is supported by numerical studies ||148| , |149|| . The same 



prediction can be made for systems with random-field (i.e. local symmetry-breaking) 
disorder |19[, and is supported by recent simulat ions [ 151| 



It is interesting that the argument leading to ( p90| ) makes no reference to whether 
the order parameter is conserved or not. The time taken to surmount the pin- 
ning barriers dominates all other timescales in the problem. The argument out- 
lined above suggests a scale-dependent kinetic coefficient T{L) ~ exp(— L^/'^^~^Y^)• 
Putting this into the usual nonconserved growth law L ~ [r(L)t]^/^ gives L ~ 
[Tln(t/L2)]x/(2-C)^ which reduces to ( p90| ) asymptotically, since InL <C Int for 
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t oo. For conserved dynamics, the same argument just gives t/L^ instead of 
t/L'^ inside the logarithm, and (|290D is again recovered asymptotically. 

While this physically based argument is certainly plausible, the RG makes a 
more powerful prediction: not only are the growth laws the same for conserved and 
nonconserved dynamics, but they belong to the same universality classl This means, 
inter alia, that they have the same scaling functionsl To see this, we simply note 
that since the fluctuations in the free energy, SF ~ L^, are asymptotically negligible 
compared to the mean, (F) ~ L'^~^ (provided the system supports an ordered phase 
at infinitesimal T) , the strong-coupling exponent y is given by the same expression, 
y = d — 1, as in the pure system. (Alternatively, and equivalently, the extra length 
of domain wall due to disorder roughening of the interfaces in a volume L'^ scales 
as ~ ];^d-3+2C ^ L'^-^y It follows from ( p73| ) that, provided the conservation law is 
relevant, the growth law is L{t) ~ t^^^ (for fi = 2) as in the pure system. Since, 

however, ( ^901 ) shows that L{t) grows more slowly than t^^^ for the nonconserved 
system, our previous arguments show that the conservation law is irrelevant for 
systems with quenched disorder. Therefore conserved and nonconserved systems 
are in the same universalit y cla ss. 

Numerical simulations ||152 , 153 | allow us in principal to test this prediction. 
They certainly s how l ogarithmic growth, but with an insufficient range of L for a 
definitive test of ( |290| ). The most striking conclusion of the RG is that the scaling 
functions are those of the nonconserved system. For example, a scaling plot for the 
structure factor, i.e. a plot of L'^S{k, t) against kL should give a non-zero intercept 
at kL = 0. For any fixed L, of course, the conservation law requires that S vanish 
at = 0, but in the scaling limit {k ^ 0, L ^ oo, with kL fixed) the region of 
small k where the conservation law is effective should shrink to zero faster than 1/L 
as L —>■ oo. There are indeed indications of this in the small-A; data of Iwai and 
Hayakawa | |153| |, but the range of L explored in not large enough to reach the true 
scaling limit. Indeed, this will always be difficult with growth as slow as ( ^9U| ). 



8.4 THE RG FOR BINARY LIQUIDS 

As a final application of the RG approach, we return to phase separation in binary 
liquids. The new element here is that the temperature T, though formally irrelevant, 
can enter scaling functions in cases where the minority phase consists of disconnected 
droplets, when the nominally dominant linear growth, due to hydrodynamic fiow, is 
absent. 

The analysis parallels that of model B, but including the extra hydrodynamic 
term of Eq. (pTl). In order to discuss the role of temperature, we have to include the 
thermal noise explicitly. The presence of the hydrodynamic term in (|5ll) implies, 
via the fiuctuation-dissipation theorem, that the noise correlator takes the form 



(e(r,t)e(r',t')) = -2\TV'Sir-r')6it-t') 

+2T V0(r) ■ r(r - r') ■ V'0(r') 6{t - t') . (291) 

Carrying out the RG step as before, the recurrence relations for A and the viscosity 
rj (implicit in the Oseen tensor (0)) become 

(292) 
(293) 

is the same as ( |271| ) with ( = d/2 and y = d — 1 inserted explicitly. The 




Eq. (|92 
renorma^ 
(with y = d — 1] 

T' = b^''^T. (294) 



isation of the noise again gives ( 272 ), which we display again for convenience 
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It is clear that the conventional conserved fixed point, A = A*, which has z = 3, 
is unstable against the introduction of hydrodynamics, since 1/r] (which measures 
the strength of the hydrodynamic interaction) scales up as 6^ at this fixed point. 
Rather, the 'hydrodynamic fixed point', with rj =J]\ must have z = 1, recovering 
the dimension-analytic result obtained in section At this fixed point A scales 
to zero, i.e. bulk diffusion of the order parameter is irrelevant at the largest scales. 
Temperature is also irrelevant, as expected. 

The physical arguments of section EO, |38[, however, show that the linear 
growth (i.e. z = 1), which is a consequence oThydrodynamic flow along interfaces, is 
possible only when the minority phase is continuous. What happens if the minority 
phase consists of isolated droplets? Then z > 1, and the relevant flxed point must 
be 7]* = 0. Let's reconsider the usual conserved ('model B') flxed point in this light. 
This flxed point, with A = A* non-zero and flnite, has z = 3, i.e. t^^^ growth. At 
this flxed point the recursion relations for rj and T are rj' = b'^rj and T' = b^^'^T. 
Therefore, rj and T both flow to zero, but their ratio remains flxed (for d = 3). Note 
that the ra tio T /r] is exactly what appears in the hydrodynamic part of the noise 
correlator (|291|) . This means that, while the temperature is technically irrelevant, 
the hydrodynamic part of the noise cannot be discarded - in fact it is just this part 
which drives the brownian motion of the droplets that is responsible for coarsening by 
droplet coalescence. The ratio is a marg'ma/ variable, so in principle we expect 
scaling functions to depend on it, reflecting the relative importance of evaporation- 
condensation and droplet coalescence to the coarsening (see, however, the discussion 
below) . 

To be more precise, we can use dimensional arguments to construct the important 
variables. The effect of thermal fluctuations on scales smaller than the correlation 
length ^ can be incorporated through the surface tension a (which scales F) and 
the equilibrium order parameter M (which scales (p). The length-scale associated 
with the Lifshitz-Slyozov mechanism is then p5| , L{t) = [Xat/M'^Y^^, while 
the dimensionless marginal variable is ksTM'^ / aXrj. The general form for L(t) is 
therefore L{t) = {Xot / M^f^ fiksTM'^ / aXi]), where f{x) is a crossover function 
with /(O) = constant. For large x, one must have f{x) ~ x^^^ so that L{t) is 

independent of A, giving L{t) ~ (kRTt/vY/^ in this regime, in agreement with the 
brownian motion argument of section |2.7| . Note that the function / also depends 
implicitly on the volume fraction v of the minority phase. 

This marginal behaviour is speciflc to d = 3. For general d, (|293|) and (|294|) 
give [T/t])' = {T /t]). So for d > 3 the ratio T /rj is irrelevant at the Lifshitz- 
Slyozov (LS) flxed point, and the evaporation-condensation mechanism dominates 
asymptotically for all T < T^, giving unique scaling functions for a given volume 
fraction. For < 3, on the other hand, T/rj is relevant at the LS flxed point. 
The dynamics is therefore controlled by a 'coalescence flxed point', with T /rj flxed, 
implying z = d, and A is an irrelevant variable. This agrees with the d = 2 result 
of San Miguel et al. [154]. To summarise, for d > 3 the LS mechanism dominates 
and L ~ t^^'^, for d < 3 coalescence dominates and z = d, i.e. L ~ t^^^, while for the 
physically relevant case d = 3 both mechanisms operate and marginal behaviour is 
expected such that, even for a given volume fraction, scaling functions will depend 
continuously on the dimensionless ratio ksTM"^ /aXrj. 

It is important to note, however, that in the above discussion the transport 
coefficient A in the equation of motion (pi]), and the viscosity 77 appearing in the 
Oseen tensor (0), have been treated as independent variables. While they can 
certainly be treated as independent in numerical simulations, in real binary liquids 
they are related [Q. Linearising (^) around one of the bulk phases gives equation 



nj ) as in model B dynamics (the hydrodynamic term in (|5T|) drops out at linear 
order). Inserting the transport coefficient on the right-handside (it was absorbed 
into the timescale in (^)) gives a bare diffusion constant Dq = XV" (1) ~ A/^^, 
where ^ is the interface thickness. The diffusion constant for a drop of size L 
is D{L) ~ Dq^/L ~ fiksT ~ ksT/riL, using the Einstein relation, and the usual 
relation fi{L) ~ ^/i]L for the mobility. This gives Xrj ~ ksT^. Using also a ~ M^/^, 
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which follows from (pT]), for the surface tension gives the crossover variable x as 
ksTM'^ /aXr] ~ 1. A more careful calculation |^0| shows that the evaporation- 



condensation and droplet coalescence mechanisms both lead to a mean droplet size 
R{t) given by = k{kBT/57iri)t, with k = 0.053 for the LS mechanism and k = 12v 
for the droplet coalescence mechanism. 

Effectively two-dimensional binary liquid systems can be achieved using the Hele- 
Shaw geometry, where the fluid is confined between parallel plates. The no-slip 
boundary conditions mean that the Navier-Stokes equation ( ^7|) simplifies, and the 
results are different from those obtained by simply putting a = 2 in the previous 
paragraph, which correspond to using 'free' boundaries. In the Hele-Shaw geometry, 
V^v is dominated by the term d'^\/dz'^, due to the rapid variation of v perpendicular 
to the plates (the z-direction). This leads to the 'Darcy's Law' form of the Navier- 
Stokes equation (the inertial terms can be neglected due to the frictional effect of 
the boundaries): 

V = {dyi2r]) {-Vp - (pVfi) , (295) 

where d is the plate spacing and \{x,y) is now the velocity averaged over the z- 
direction. Using the incompressibility condition V ■ v = to eliminate p yields an 
equation of the form (|50D, but with the Oseen tensor replaced by its Hele-Shaw 
equivalent 



Our starting point for the RG analysis is therefore equations (P7|) and ( P91| ), with 
the Oseen tensor T replaced by T^^ . 

Coarse-graining in the xy-p\ane, holding the plate spacing d fixed, gives the 
recurrence relations 

A' = b'-^X 

T' = b-^T , (297) 



which have the same form as (|292|) - ( p94| ), with the extra factor 6^ in the rj equation 



corresponding to the extra factor of k"^ in T^^ . Eqs. (|297|) still have the LS fixed 
point A = A*, with z = 3 and T irrelevant, but now t] (or the product Xt]) is marginal, 
suggesting again a continuous family of universality classes reflecting the relative 
importance of bulk diffusion and hydrodynamic flow. In their numerical studies of 
critical quenches in the Hele-Shaw geometry, Shinozaki and Oono ||155|| verify the 



^1/3 gj^Q-y^^i;}^ a^j^d find that the scaling functions do indeed depend systematically on 
the value of A?; (and propose essentially the same explanation). 

For any off-critical quench, fluid flow along the interfaces terminates when the 
droplets of minority phase become circular. Eventually, the LS mechanism domi- 
nates the coarsening, with L{t) ~ g^iu^ 

unique scaling functions for a given 
volume fraction. In this geometry the ratio T/t], representing the hydrodynamic 
noise, flows to zero at the LS flxed point, {T/rj)' = b^~'^{T/ri), so droplet coales- 
cence is subdominant asymptotically in time. If the transport coefficient A is small 
enough, however, the coalescence fixed point, with z = A so that T/rj is fixed, will 
dominate the coarsening for a range of times, giving L(t) ~ {Tt/viY^^. 

This result may be derived heuristically by extending to the Hele-Shaw geometry 
the argument given in the final paragraph of section |2.7| . For droplets of size R with 
areal number density n ~ vjR?^ where v is the volume fraction, the 'coalescence 
time' is given by the same expression, ~ R^/vD, as in the bulk case. In the 
Hele-Shaw geometry, however, the mobility of a droplet of size i? is /i ~ d/rjR^, so 
the Einstein relation D = ksT^ for the diffusion constant gives tc ~ rjR'^/vdksT. 
This implies a time-dependence R ~ [vdkBTt/T])^/'^ for the typical radius of a drop. 
Comparing this with the LS growth, R ~ (Xat/M'^Y^^, shows that the crossover from 
coalescence dominated to LS dominated regimes occurs when R ~ vdknTM^ /r]Xa. 
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This length is just vd times the dimensionless crossover variable x we identified in 
the discussion of bulk binary liquids. For real binary liquids (as opposed to computer 
simulations) we have seen that x is of order unity, so the crossover length is set by 
the product vd of the volume fraction and the plate spacing. Since this product 
is obviously less than d, it follows that a t^^^ coalescence regime (which requires 
d) should be unobservable in real binary liquids. 



9 SUMMARY 



In this article I have reviewed our current understanding of the dynamics of phase or- 
dering, and discussed some recent developments. The concept of topological defects 
provides a unifying framework for discussing the growth laws for the characteristic 
scale, and motivates approximate treatments of the pair correlation function. 

The most important consequence of the presence of topological defects in the 
system is the 'generalized Porod law', equation (ffOD, for the large kL tail of the 
structure factor. This power-law tail, whose existence has long been known for 
scalar sy stems, has rece ntly been observed in computer simulations of various vector 
systems ||105| , |128| , |156|| . It should be stressed that the form of the tail depends only 
on the nature of the dominant topological defects. In nematic liquid crystals, for 
example, the presence of disclinations or '1/2-strings', implies a structure factor 



67, 66 



6[, i.e. a. k tail for bulk systems. This tail 



, and is not inconsistent |^ with experimental 



tail described by ([70|) with n 
has been seen in simulations 
results ||6^, |63| . 

The Porod law (^), together with the scaling hypothesis, leads to a powerful 
and general technique for deriving growth laws [Q. The results are summarized in 
Figure 24. Again, the technique is more general than the simple 0{n) models to 
which it has been applied here. Nematic liquid crystals, for example, are described 
by the nonconserved dynamics of a traceless, symmetric, tensor field. However, the 
presence of dominant string defects implies the same growth law as for the 0(2) 
model, namely L{t) ~ t^/^, consistent with the simulations l |67|| (a llowing for the 
predicted logarithmic corrections to scaling) and experiment [p2|, |63[|. 

The dominant role of topological defects also motivates approximate treatments 
of the pair correlation scaling function /(x) 0, |5^, |5^, 60|, |85|, and the systematic 
treatmen t [pOl] discussed in section ^. All of t hese theories lead to the same scaling 
function ( |131|) , with the O JK scaling function (|126|) corresponding to the special case 

1. The form ( |131| ) is a direct consequence of the non-linear mapping 0(m), with 

— i> m for |m| oo, and the gaussian distribution assumed for the field m. The 
'OJK-type' theories |5|, |8|, g and the 'Mazenko-type' theories p|, 0, |8|, |8g 



n 



differ only in the equation for 7, the normalized pair correlation function for rh. 

These approximate scaling functions all give good fits to experiment and simu- 
lation data (see, e.g.. Figures 14 and 15). However, there is one important caveat. 
When fitting data to theoretical scaling functions, it is conventional to adjust the 
scale length L{t) for the best fit. An absolute test can, however, be obtain ed by 

m3 



For example, the normalized correlator 
calculated within 'gaussian' theories of t 



calculating two different scaling function s an d plotting one against the other 

of the square of the field can also be 

m 



leOJK or Mazenko type 



The result 



depends only on 7, the normalised correlator of the gaussian auxiliary field. Elim- 
inating 7 between C(12) and 6*4(12) gives an absolute prediction for the function 
C4(C). When this prediction is compared to simulation results, however, the agree- 
ment is found to be rather poor (Figure 17): C and C4 can be fitted separately, as 
functions of r/L(t), by choosing the scale length L(t) independently for each fit, but 
not simultaneously. However, the agreement improves with increasing d, in agree- 
ment with the idea that these theories based on a gaussian auxiliary field become 
exact at large d |Q. Including the 1/A^ correction in the systematic approach of 
section |5]^ will presumably improve the fit at fixed d. Mazenko has r ecently in- 
troduced an alternative way of including non-gaussian corrections [|102|| , and finds 
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improved agreement with the simulation results. 

The calculation of scaling functions for conserved fields is a significantly greater 
challenge, especially for scalar fields, where even obtaining the correct t^^^ growth 
law, within an approximate theory for the pair correlation function, is not straight- 
forward. Mazenko has extended his approximate theory to conserved scalar fields 
89| , but the agree ment with high quality simulation data is not as good as for 



nonconserved fields [|114|| . There is an additional complication that a naive applica- 



tion of Mazenko's method gives t^^^ growth, which Mazenko argues corresponds to 
surface diffusion only. In order to recover the t^^^ growth, he has to add an addi- 
tional term to incorporate the effect of bulk diffusion. For conserved vector fields, 
the naive Mazenko approach gives the expected t^^^ growth (see Figure 24), but 
without the logarithmic correction expec ted fo r n = 2. The approximate analytic 
treatment [|79l |115|| (presented in section |5.5.1| ) of the equation for C(12), valid for 



5.2, 



n ^ 1, gives good agreement with scaling functions extracted from simulations 

A systematic approach for conserved fields, generalizing the treatment of section 

would be very welcome, although it is far from straightforward. An even greater 
challenge is to develop good approximate scaling functions for binary liquids. 

To summarize, we have focussed on the role of topological defects as a general 
way of deriving, through the Porod law (|7D|) and the scaling hypothesis (represented 
by equations (|^) and (W)) the forms of the growth laws for phase ordering in various 
systems. The study ofsuch defects also motivates, through the mapping to an aux- 
iliary field that varies smoothly through the defect, approximate theories of scaling 
functions. For nonconserved fields, such methods are, in principle, systematically 
improvable (section One of the challenges for the future is to try to develop 
comparable methods for conserved fields. 

From a wider perspective, phase ordering dynamics is, perhaps, the simplest 
example of a scaling phenomenon controlled by a 'strong-coupling' RG fixed point 
(Figure 4). It may not be too much to hope that techniques developed here will find 
useful applications in other branches of physics. 
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Figure Captions 

Figure 1 Magnetization of the Ising model in zero applied field as a function of tem- 
perature (schematic), showing spontaneous symmetry-breaking at Tc- The arrow 
indicates a temperature quench, at time t = 0, from Tj to Tp. 

Figure 2 Monte Carlo simulation of domain growth in the d = 2 Ising model at 
T = (taken from J. G. Kissner, Ph.D. thesis. University of Manchester, 1992). 
The system size is 256 x 256, and the snapshots correspond to 5, 15, 60 and 200 
Monte Carlo steps per spin after a quench from T = oo. 

Figure 3 Typical form of the symmetric double- well potential V{(j)) in equation (Q). 
The detailed functional form of V{(f)) is not important. 

Figure 4 Schematic Renormalization Group flow diagram, with fixed points at T = 0, 
Tc and oo. All T > Tc are equivalent to T = oo and all T < Tc to T = 0, as far as 
large length-scale properties are concerned. 

Figure 5 Scaling function f{x) for the pair correlation function of the d = 2 Ising 



model with nonconserved order parameter (from reference [97]). The time t is the 
number of Monte Carlo steps per spin. 

Figure 6 Domain- wall profile function (j){g) (schematic). 

Figure 7 Asymmetric potential V{(j)) for a conserved order parameter, showing the 
common-tangent construction that determines the compositions of the separated 
phases. 

Figure 8 Sketch of the function g{x), given by Eq. (H^, for different 7, where 70 = 



4/27. 

Figure 9 The 'Mexican hat' potential V{(j)) for the 0{n) model with n = 2. 

Figure 10 Types of topological defect in the 0{n) model: (a) domain wall (n = 1) 
(b) vortex {n = 2 = d) (c) string {n = 2, d = 3) (d) monopole, or 'hedgehog', 
(n = 3 = d) (e) antivortex. 

Figure 11 Cross section of ±l/2-string configurations for a nematic liquid crystal. 

Figure 12 Scaling function f{x) for the nonconserved d = 2 Ising model, showing 

Monte Carlo data (MC) from Figure 5, and the approximations of OJK and 
Mazenko [B^. The scaling lengths L{t) for the theoretical curves were chosen to give 
the same slope as the data in the linear 'Porod' regime at small x (from reference 

0). 



Figure 13 Spatial variation (schematic) of the order parameter and the auxiliary 
field m, defined by Eq. ( [T32|) . 

Figure 14 Scaling plots for the pair correlation function of nonconserved systems 
with an 0(ri)-symmetric vector order parameter, plotted against rp^/" where p is 
the defect density (proportional to (1 — 0^) for n = 1). The data are taken from 
reference ||105|| . In (d), the length scale L{t) was chosen independently at each time 
to give the best collapse. The continuous curves are 'best fits by eye' of the BPT 
prediction [|^, |58|] . 

Figure 15 Log-log scaling plots of the structure factor for nonconserved systems 
with an 0(n)-symmetric vector order parameter. The data are taken from reference 



105|| . The continuous curves are the Fourier transforms of the corresponding curves 
in Figure 14. The data exhibit the expected k~'^'^^^^ tails for large kL{t). 

Figure 16 Schematic forms of (a) the potential Vi(0)and (b) the total potential V{(f)) 

used to incorporate external fields into the systematic approach (section |5.2.4| ). The 
dashed lines indicate parts of the potential that are irrelevant to the dynamics 
described by Eq. ( |170| ). 
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Figure 17 An 'absolute test' for theories of nonconserved dynamics based on an 
assumed gaussian auxiliary field. Here C and C4 are the pair correlation functions 
for th e ord er parameter and it square, the latter normalised by its large- distance limit 
fEq. (|176 )). The data are for a scalar order parameter in dimension (a) d = 2, and 



(b) d = 3. The continuous curve (independent of d) is the prediction of gaussian 
theories based on the OJK or Mazenk o ap proaches. The broken lines give the 
predicted short- distance behaviour (see [|105|| and section ||). 

Figure 18 Same as Figure 17, but for vector fields: (a) d = 2 = n, {h) d = 2, 
n = 2, and (c) d = 3 = n. The continuous curves are the predictions of the gaussian 
theories. 

Figure 19(a) Real-space pair correlation function for a nematic liquid crystal within 
the 'equal- const ant approximation', calculated using the 'KYG' approach as de- 
scribed in the text. The scaling variable x is r/y/si. The data are the Monte Carlo 
simulations from reference ||6^, with L{t) fixed from the best fit to the theory. Inset: 
short- distance behaviour of the theory, showing a leading In x singularity. 

Figure 19(b) Log- log plot for the scaled structure factor of a nematic liquid crystal. 

Continuous curve: the 0(2) theory; data points: simulation data from reference 
|67| , with L{t) chosen as in Figure 19(a). Experimental data from reference []63 



are 



shown on the left, arbitrarily positioned: they can be moved left-right and up-down. 
The straight line is a guide to the eye, with slope -5. 

Figure 20 Coordinate system employed for the calculation of the amplitude of the 
Porod tail for d = 3, n = 2. 

Figure 21 Simulation data from Figure 15 replotted to reveal the amplitude A{n, d) 
of the Porod tail, defined by 5(k,t) A{n,d)pdef/k'^^'' for kL{t) > 1. The hori- 
zontal dashed lines are the prediction of Eq. ( p21| ) for the asymptotic limit. 

Figure 22 Relative positions of points '1' and '2' for the small-r limit of C4(12). 

Figure 23 Calculation of the short-distance limit of P (m(l),m(2)). 

Figure 24 Time-dependence of the characteristic scale L{t) for systems with purely 
dissipative dynamics. Exceptional cases are discussed in the text. 

Figure 25 'Long-range' contribution to the structure factor, for long-range correla- 
tions in the initial conditions. The circle represents the initial con ditio n correlator, 
Eq. ( p77| ), while the external legs represent the response function (|278|) . 
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